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1.0 INTRODUCTION 

This report includes data for the 1972 (October, 1971 to June, 1972) 
and 1973 (October, 1972 to June, 1973) growing season for winter wheat in 
Kansas. Obviously, there was no ERTS— 1 imagery for the 1972 wheat crop 
since the satellite was launched July 23, 1972. However, under the PEIS 
(Pre-ERTS Investigator Support), we obtained three U-2 flights over our 
test areas (March 21, April 26 and June 6, 1972). The imagery was taken 
during the NASA— Ames to Wallop Island ferry trips. U— 2 flight lines for 
each date differed and made chronological comparison of a given test area 
impossible; however, we were able to detect virus infected wheat fields on 
U-2 color infrared film. A systematic analysis of U— 2 data was not possible 
because ground observations and U-2 flight lines did not always coincide. 
Since the actual U-2 film was not received until 3 months after the harvest, 
we could not generate the necessary ground truth. 

In addition to the U-2 flights, the University of Kansas-CRINC (Center 
for Research, Inc.) provided three low altitude flights with a Cessna 182 
equipped with four Hasselblad cameras, Cessna flight coverage was selected 
from ground observations of wheat fields. Diseased and healthy fields in 
the same frame were compared as to the date of planting, variety, disease 
severity, fertility, topography and crop rotation. In comparing many of 
healthy and diseased fields, a common denominator was crop rotation. For 
example, soil borne mosaic virus was more severe on fields that were cropped 
with wheat the previous year. 

The 1973 Kansas wheat crop produced a record 381,000,000 bu with an 
average yield of 37.0 bu per acre which is attributed to a 10% increase 
over 1972 in both harvested acres and yield. Part of the increase in yield 
can be attributed to abundant soil moisture because of ample precipitation 
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and the absence of severe disease epidemics over the state# The abnormally 
high precipitation adversely affected this project by (a) reducing water 
stress , (b) producing cool weather which affects disease infestation, and 
(c) association with clouds which affect quality of the imagery. 

Efforts were concentrated in two of the five test areas: Finney County 

(Garden City) and Riley County (Manhattan), (Fig. 1.1). The Finney county 
site had a large area infected with wheat streak mosaic virus (WSMV). Riley 
county was the location of the data collection systems (DCS) and the plots 
for the energy balance determinations. We obtained 5 ERTS— 1 observations 
for the WSMV area and 7 observations for the irrigated and nonirrigated 
wheat fields in Finney county. In addition, 8 observations were obtained 
for Riley county. Because of cloud cover during U-2 flights, no U-2 
imagery was available for the entire 1973 wheat crop. Two sets of U-2 
imagery were obtained after the wheat harvest. 

Objectives of the investigation were: (a) to evaluate the effect of 

water stress, disease, and leaf area on the reflectance characteristics of 
wheat, (b) to evaluate disease losses in terms of yield and water use, and 
(c) to predict disease severity and economic loss. 


In this report, the water use of irrigated and nonirrigated wheat and 
of healthy and diseased wheat are discussed in Chapter 2 (objective a and 
b). Data handling and reduction of ERTS-1 imagery is reported in a separate 
chapter (Chapter 3). The use of ERTS-1 imagery to determine leaf area, _ 
crop growth and disease severity is given in Chapter 4 and 5 (objective a 


and b) . Prediction of disease severity and economic loss (objective c) by 
the use of ERTS-1 data collection systems is reported in Appendix A. Appendix 


C contains a detailed electronic description of the data collection interface 
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Fig. 1.1 Map of test areas in Kansas. 
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system. The results of ground measurements of the spectral reflectance of 
sorghum, soybean and wheat are provided In Appendix B. As a result of the 
ERTS-1 project, a Master's Thesis was written and is provided as Appendix 
D. Appendices E, F, G and H are related to data handling and reduction 
phases of the project. 

The investigators wish to express appreciation to the National Aero- 
nautics and Space Administration for supporting this project. 



6 


Soils on the two fields are classified as a Ulyssess-Rich field silt 
loam with 1.5 percent organic matter and a pH of 6.9. A particle size 
analysis indicated the soils contained about 50 percent silt and 20 percent 

clay. 

Field A, which had been summer fallowed since July 1971, was planted 
to Scout winter wheat on September 15, 1972. Seeding rate was 29 kg per 
hectare in rows spaced 25.4 cm apart. The wheat was completely headed by 
May 24, 1973 and a yield of 2689 kg per hectare harvested on July 5, 1973. 

As field B had been in wheat the previous season, it was preirrigated. 
Ninety kg per hectare of nitrogen were applied as anhydrous ammonia prior 
to planting and 50 kg per hectare of Eagle winter wheat were seeded on 
September 22, 1972, in rows 30.48 cm apart. Two irrigations of 3.05 cm 
each were applied through a center pivot sprinkler system on May 23 and 
June 2, 1973. Wheat harvest was completed on July 5 and the yield was 
3496 kg per hectare. 

According to Variety Tests with Fall-Planted Small Grains (1971), 
Eagle wheat is a selection of Scout with nearly identical vegetative 
characteristics. Reflectance measured by the MSS system aboard ERTS-1 
should be independent of winter wheat variety. Scout or Eagle. 

Data Collection 

Fields A and B were both divided into four equal sized square plots 
with a sampling area in the center of each plot. Two additional sampling 
areas were established in field A where the corners of the field were 
double drilled. The sampling areas were sub-divided into one meter square 
plots. Leaf area index and soil moisture were measured within one day 
of each pass of ERTS-1 on one plot of each sampling area selected at 


random • 



7 


Leaf area was determined by measuring the length and breadth of each 
leaf from randomly selected plants. Area of each individual leaf was 
calculated from the equation (Teare and Peterson, 1971): 

LA - Q.813X - 0.64 <2 * 1) 

where : 

„ 2 
LA = Leaf area, cm 

2 

X = Product of length times breadth of leaf, cm 
Total leaf area on each one meter square plot was calculated by multiplying 
number of plants times average plant leaf area. Leaf area index was taken 
as the ratio of total leaf area to the land surface area. 

Soil samples were taken from the surface and at the following incre- 
ments of depth: 0 to 15, 15 to 30, 30 to 60, 60 to 91, 91 to 121, 121 to 
152, and 152 to 182 cm. The samples were weighed, dried in an oven at 105 C 
until they reached a constant weight and reweighed. Soil moisture was cal 
culated and expressed as percent on a dry weight basis. 

Bulk density, field capacity and the permanent wilting point for 
Ulyss ess -Richfield silt loam were obtained from the Garden City Experiment 

Station. 

Meteorological data from the area were collected for use in the ET model. 
Maximum and minimum temperatures, dew point temperatures and wind run were 
obtained from the Garden City Experiment Station. Solar radiation was 
obtained from the Dodge City Weather Service and rainfall was measured near 
the field site. 

Data Analysis 

Digital counts in each MSS band and various combinations of bands were 
compared by multiple regression techniques with leaf area index and available 


soil moisture. 
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Using an IBM 360/50 digital computer, estimates of the water use by 
wheat from the ET model developed by Jensen et al, (1971) were compared 
with changes in measured soil moisture. 

RESULTS 

Although ERTS-1 passed over the field site every 18 days, clear atmos- 
pheric conditions were encountered on only 6 flight days during the period 
from wheat seeding to wheat harvest. Table 2.1 gives weather conditions on 
each flight date and the schedule of data collection. Data for July 7 were 
excluded from the analyses because the wheat crop had been harvested. 

Prediction of Leaf Area Index 

Tables 2.2 and 2.3 give digital counts of MSS data for bands 4, 5, 6 
and 7 and ratios of 4 to 5, 4 to 7 and 5 to 7 along with measured LAI for 
fields A and B. The regression equations for predicting LAI are: 


LAI = -0.15MSS4 + 5.41 

* 

R 2 

= 0.80 

(2.2) 

LAI - -0.065MSS5 +2.66 

> 

R 2 

= 0.86 

(2.3) 

LAI = 1.94MSS6 - 9.37 

> 

R 2 

= 0.20 

(2.4) 

LAI = 0.15MSS7 - 3.53 

f 

R 2 

= 0.53 

(2.5) 

LAI = 2.92MSS4/5 - 2.63 

9 

R 2 

» 0.95 

(2.6) 

LAI = -1.22MSS5/7 + 2.08 

> 

R 2 

« 0.85 

(2.7) 


where 

LAI “ Leaf area index 

MSS « Digital counts for numbered band or ratio 
2 

R «= Regression coefficient 

For each equation there is some minimum or maximum value of the digital 
counts or ratios beyond which LAI goes negative and the results are meaningless. 



Table 2.1 Weather conditions during satellite pass over 
test fields. 


Date 


Weather 

Condition 

Data 

Acquired 

September 


1972 

Cloudy 


September 22, 

1972 

Clear 

X 

October 

*s 

o 

f-l 

1972 

Partly Cloudy 


October 

28, 

1972 

Cloudy 


November 

15, 

1972 

Cloudy 


December 

3, 

1972 

Partly Cloudy 


December 

21, 

1972 

Partly Cloudy 


January 

8, 

1973 

Cloudy 


January 

26, 

1973 

Cloudy 


February 

13, 

1973 

Rain 


March 

3, 

1973 

Foggy 


March 

21, 

1973 

Clear 

X 

April 

8, 

1973 

Heavy Snow 


April 

26, 

1973 

Rain 


May 

14, 

1973 

Clear 

X 

June 

1, 

1973 

Clear 

X 

June 

19, 

1973 

Clear 

X 

July 

7, 

1973 

Clear 

X 


indicates both ERTS-1 and field data taken. 



Table 2.2. Digital counts of MSS data and LAI for field A 


Date 


MSS4 

MSS 5 

MSS6 

MSS 7 

MSS4/5 

MSS4/7 

MSS5/7 

LAI b 

9/22/72 

Mean 

34.75 

37.89 

38.64 

19.55 

0.918 

1.779 

1.939 

0.00 


S.D. a 

1.41 

1.90 

2,18 

0.86 

0.040 

0.068 

0.080 

0.00 

3/21/73 

Mean 

33.26 

32.29 

45.87 

25.25 

1.031 

1.318 

1.280 

0.37 


S.D. 

1.28 

1.58 

1.74 

0.69 

0.040 

0.055 

0.069 

0.10 

5/14/73 

Mean 

29.74 

24.50 

48.11 

28.08 

1.218 

1.064 

0.877 

0.97 


S.D. 

1.69 

2.12 

1.79 

1.66 

0.066 

0.101 

0.104 

0.26 

6/1/73 

Mean 

33.43 

29.48 

52.32 

29.87 

1.138 

1.121 

0.990 

0.89 


S.D. 

1.72 

2.42 

1.84 

1.04 

0.062 

0.083 

0.104 

0.25 

6/19/73 

Mean 

41.14 

49.33 

55.26 

28.70 

0.835 

1.436 

1.722 

0.00 


S.D. 

1.62 

2.07 

1.49 

0.92 

0.033 

0.074 

0.090 

0.00 

7/7/73 

Mean 

59.46 

78.53 

77.68 

36.36 

0.758 

1.636 

2.161 

0.00 


S.D. 

2.14 

4.25 

2.72 

1.49 

0.030 

0.061 

o.iis 

0.00 


^Standard deviation. 
u Average of six sampling points 



Table 2.3, Digital counts of MSS data and LAI for field B 


Dace 


MSS4 

MSS5 

MSS6 

MSS7 

MSS4/5 

MSS4/7 

MSS5/7 

LAI b 

9/22/73 

Mean 

37.05 

40.41 

40.96 

20.78 

0.919 

1.786 

1.947 

0.00 


S.D. 

1.62 

2.54 

2.37 

1.02 

0.038 

0.094 

0.128 

0.00 

3/21/73 

Mean 

33.54 

32.99 

41.47 

22.57 

1.019 

1.488 

1.463 

0.44 


S.D. 

1.09 

1.91 

2.15 

0.96 

0.049 

0.073 

0.088 

0.07 

5/14/73 

Mean 

27.63 

19.22 

56.66 

36.78 

1.454 

0.760 

0.532 

1.53 


S.D. 

1.60 

2.68 

3.56 

3.18 

0.132 

0.109 

0.129 

0.39 

6/1/73 

Mean 

26.93 

20.03 

48.66 

31.61 

1.355 

0.858 

0.638 

1.23 


S.D. 

1.32 

2.23 

3.43 

2.54 

0.111 

0.083 

0.094 

0.36 

6/19/73 

Mean 

36.68 

37.94 

52.00 

29.97 

0.971 

1.227 

1.270 

0.00 


S.D. 

1.21 

3.11 

2.05 

1.56 

0.060 

0.079 

0.131 

0.00 

7/7/73 

Mean 

54.46 

73.87 

77.48 

38.24 

0.739 

1.425 

1.932 

0.00 


S.D. 

2.30 

4.37 

3.39 

1.31 

0.033 

0.060 

0.100 

0.00 


a 

k Standard deviation. 

Average of four sampling points. 
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The regression coefficients for all the LAI prediction equations except 
equation (2.4) were statistically significant at the 0.05 level. For the 
individual bands, 4 and 5 were better predictors of LAI. While ratios of 
band 4 to band 5 and band 5 to band 7 were good predictors of LAI, the ratio 
of band 4 to band 5 had the higher regression coefficient. Figure 2.1 shows 
measured LAI as a function of the ratio of digital counts in band 4 and band 
5. 

Prediction of Soil Moisture 

Available soil moisture for plant use was taken as the difference between 
soil moisture at sampling and at the permanent wilting point. Table 2.4 gives 
available soil moisture for fields A and B at depths of 0 to 15, 0 to 30, 0 to 
60, and 0 to 91 cm. The negative values are due to errors in the assumption 
of soil moisture percentage at permanent wilting point as the soil profile 
probably did not become that dry. Rather than adjust soil moisture percentage 
by changing the permanent wilting point percentages, the available soil mois- 
ture values were left negative. The only effects of this action are that 
the constants in the regression equations may be in error but the form of 
the equations will remain the same and the regression coefficients are not 
affected. 

Linear regression equations to predict available soil moisture from 
digital counts in the individual MSS bands and various band ratios were 
developed from digital counts in Tables 2.2 and 2.3 and available soil moisture 
in Table 2.4. Data for field B on March 21 were excluded from the analyses 
because of rain between flight of ERTS-1 and soil measurements. Only band 6 
predicted available soil moisture at the 0.10 level of significance. The 
equations for available soil moisture in the 0 to 15 cm zone, AV^, and in 
the 0 to 91 cm zone, AVgp 


were: 



MEASURED LEAF AREA INDEX 




Table 2.4. 


Available soil moistures (cm of 
water) . 


Date 


Soli Depth- an 


0-15 

0-30 

0-60 

0-91 


FIELD A (DRYLAND) 


9/22/72 

0.66 

1.37 

2.59 

3.29 

3/21/73 

0.66 

1.50 

3.13 

4.33 

5/14/73 

0.19 

0.40 

1.24 

2.38 

6/1/73 

-0. 26 a 

-0.34 

-0.14 

0.23 

6/19/73 

-0.63 

-0.99 

-1.46 

-2.02 


FIELD 

B (IRRIGATED) 


9/22/72 

0.61 

1.00 

1.58 

1.70 

3/21/73 

0.73 

1.71 

3.38 

4.48 

5/14/73 

0.52 

0.67 

1.35 

2.16 

6/1/73 

0.85 

1.20 

1.44 

1.51 

6/19/73 

-0.37 

-0.78 

-1.12 

-1.57 


Negative available soil moisture because soil 
moisture at permanent wilting point was lower 
than that measured at the Garden City Branch 
Experiment Station. 
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AV 15 = 7.52 - 0.15MSS6 , R 2 = 0.35 (2.8) 

AV 91 = 30.31 - 0.60MSS6 , R 2 - 0.33 (2.9) 

Measured LAI was added as an independent variable and multiple regres- 
sion equations developed to predict soil moisture. The resulting equations 
for available soil moisture in the 0 to 15 cm zone, AV^,-, and the 0 to 91 cm 
zone, AV^, which had a significant regression coefficient at the 0.05 level. 


were: 

AV 15 = 11.73 - 1.88 LAX - 0.31MSS4 , R 2 = 0.93 (2.10) 
AV 15 = 5.14 + 0.76 LAI - 0.11MSS6 , R 2 - 0.80 (2.11) 
AV 15 = 4.85 + 1.56 LAI - 0.20MSS7 , R 2 = 0.79 (2.12) 
AV 91 = 21.99 + 3.03 LAI - 0.46MSS6 , R 2 = 0.73 (2.13) 
AV 91 = 18.1 + 5.60 LAI - 0.71MSS7 , R 2 = 0.84 (2.14) 


As LAI was highly correlated with the ratio of counts in band 4 to 
band 5, equation (2.6) was substituted into equations (2.10) and (2.14) for 
predicting available soil moixture, The resulting equations were: 

AV 15 * 16 - 67 " 5.49MSS4/5 - 0.31MSS4 (2.15) 

and AV 91 = 3.37 + 16.35MSS4/5 - 0.71MSS7 (2.16) 

Figure 2.2 shows the comparison of available soil moisture predicted 
by equations (2.15) and (2.16) with measured available soil moisture. 

Predicting Water Use 

As no usable data were collected by ERTS-1 from the time wheat was 
seeded until March 21, 1973, water use calculations were made from that 
date until just prior to wheat harvest. Table 2.5 gives precipitation 



PREDICTED AVAILABLE SOIL MOISTURE (cm) 
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MEASURED AVAILABLE SOIL MOISTURE (cm) 

Figure 2.2. Comparison of measured and predicted available 
soil moisture. 


Table 2.5. Precipitation and irrigation at the test 
sites. 


Date 

Precipitation 

cm 

Irrigation (Field B) 
cm 

3/23/73 

1.10 


3/27/73 

0.70 


3/30/73 

1.00 


4/7/73 

0.25 


4/24/73 

0.80 


5/7/73 

1.25 


5/23/73 


3.05 

6/2/73 


3.05 

6/28/73 

0.90 
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measured by a local fanner near the test site and irrigation water applied 
to field B during the latter period. The irrigations were carried out over 
about a 3 day period and the total amount credited to the center day of the 
period. 

The ET model as developed by Jensen et al. (1971) was used to predict 
soil moisture depletion on both wheat fields using their suggested crop 
coefficients. Measured soil moisture depletion was entered into the 
computer program for March 21, 1973. The wheat crop coefficients were: 

Y = 0.233 - 0.0114X + 0.000484X 2 - 0.00000289X 3 (2.17) 

and Y = 1.022 + 0.00853D - 0.000726D 2 + 0.00000444D 3 (2.18) 

where: 

Y = Wheat crop coefficient 

X = Percent of period between seeding and 100 percent crop cover 

D = Days after 100 percent crop cover 

Table 2.6 gives the measured and predicted soil moisture depletion and water 
use for fields A and B where measured soil moisture depletion is the difference 
between field capacity for the 182 cm soil profile sampled and the measured 
soil moisture level. The ET model over-predicted soil moisture depletion 
and water use for the period March 21 to June 1 and greatly underpredicted 
soil moisture depletion and water use for the period June 1 to June 19. 

Soil moisture depletion and water use for the entire period were under- 
predicted by the ET model. 

A new crop coefficient for the ET model was computed from a multiple 
regression analysis using LAI measurements from field A. The new crop 
coefficient equations were: 



Table 2.6. Measured and predicted soil moisture depletion and water use 


Date 

Soil Moisture Depletion 

- cm 

Water Use 

for Period - cm 

Measured a 

Jensen* 5 

Revised 

c d 3 . 

1 Revised 2 Measured 

Jensen^ 

c d 

Revised 1 Revised 2 





FIELD A (DRYLAND) 



*3/01 /7Q 







J/ LLj / 5 

1/ • to 






5/14/73 

19.84 

25.07 

19.35 

7.29 

12.52 

6.80 

6/1/73 

27.74 

31.24 

27.86 

7.90 

6.17 

8.51 

6/19/73 

37.52 

32.16 

34.65 

9.78 

0.92 

6.79 

Total 

__ 

_ — 


24.97 

19.61 

22.10 


3/21/73 

19.28 



FIELD B (IRRIGATED) 












5/14/73 

23.44 

26.14 

20.80 

26.52 

9.26 

11.96 

6.62 

12.34 

6/1/73 

27.15 

29.24 

26.56 

32.66 

6.76 

6.15 

8.81 

9.19 

6/19/73 

37.77 

27.61 

31.52 

37.90 

13.67 

4.68 

8.01 

8.29 

Total 





29.69 

22.79 

23.44 

29.82 


Field measurements of soil moisture depletion. 
^Wheat crop coefficient suggested by Jensen et al . 
C Wheat crop coefficient from LAI of field A, 

^Wheat crop coefficient from LAI of field B. 


(1971). 


M 

lO 
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Y = 0.005 + 0.0165X - 0.000467X 2 + 0.00000402X 3 (2.19) 

and Y = 0.998 - 0.00297D - 0.000747D 2 (2.20) 

Predicted soil moisture depletion and water use are given in Table 2.6 
for fields A and B using the revised crop coefficients in the ET model. Pre- 
dicted water use for the period was within 2.87 cm or about 10 percent of 
measured water use for field A and is within the accepted accuracy for the 
ET model. Predicted soil moisture depletion for field A compared very closely 
with measured values on May 14 and June 1 but was low on June 19 when wheat 
was nearing maturity. The ET model with a revised crop coefficient based 
upon LAI for field A under-predicted water use for field B. 

As the ET model utilizing LAI as the crop coefficient from dryland wheat 
successfully predicted water use on field A, LAI from irrigated wheat was 
used as the crop coefficient to predict water use on field B. The crop 
coefficient equations were: 

Y = 0.0109X - 0.000288X 2 + 0.00000333X 3 (2.21) 

and Y = 1.52 - 0.000834D 2 (2.22) 

Table 2.6 gives the predicted soil moisture depletion and water use for 
field B using the revised crop coefficient in the ET model. Predicted soil 
moisture depletion and water use were about equal to measured values on 
June 19. However, there was considerable variation between measured and 
predicted values on May 14 and June 1. 

DISCUSSION 

Limited data were collected during this study because of cloud cover 
during many passes of ERTS-1. However, it appears that the MSS system has 
the potential for predicting water use of growing crops. 
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One possible method for predicting water use is from available soil 
moisture predicted from reflectance measurements. Water use would be the 
difference between available soil moisture on succeeding days. As the MSS 
system only sees the earth’s surface, soil moisture would be predicted from 
its effects on the soil surface and growing vegetation. 

Kondrat’yev (1965) reported that albedo varies between soils. The 
variability was attributed to different soil color, soil moisture content, 
organic matter content and soil particle size with soil moisture content 
the most important factor. Bowers (1971) found that reflectance increases 
as soil moisture decreases and concluded that reflectance techniques are 
precise enough to measure surface moisture. However, due to the effects of 
other soil factors on reflectance, a calibration will be necessary for each 
soil type. 

There are several factors which influence reflectance from growing 
vegetation. According to David (1969), a water deficit in the soil will 
result in increased reflectance. Severe nitrogen deficiencies also increase 
reflectance (Remote Sensing, 1970). Leaf reflectance is affected by variety 
and relative maturity of the crop (Remote Multispectral Sensing in Agriculture, 
1970). There are other factors including soil salinity, plant diseases and 
mineral deficiencies which may affect reflectance from vegetation. Whether 
soil moisture can be accurately predicted by the MSS system depends upon the 
relationships between soil moisture, vegetative growth and factors affecting 
reflectance. 

Only band 6 showed potential for predicting available soil moisture by 
a linear relationship with digital counts. The addition of a second band or 
band ratio which correlated with LAI improved the accuracy of available soil 
moisture prediction. 
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Water use of growing vegetation can be predicted from ET models. The 
MSS system, through prediction of LAI, has the potential to supply the 
numerical values of the crop coefficient equations for winter wheat developed 
from LAI measurements. Although LAI is greater for irrigated wheat than 
for dryland wheat, the curves have the same general shape as shown in 
Figure 2.3. Additional research is needed to determine the correct relation- 
ship between LAI and crop coefficient. When this relationship can be ex- 
pressed mathematically, the MSS system will be capable of supplying the 
crop coefficient for ET models. 

Practical use of data collected by ERTS-1 or similar vehicles is depen- 
dent upon timely acquisition and processing of the data. To be useful in 
water resources management, the data should be available within 24 hours of 
flight time. One potential use of the data is in irrigation scheduling. 

SUMMARY 

To effectively manage water on agricultural lands, daily water use of 
crops must be known. We hypothesized that the MSS system aboard ERTS-1 
could provide data for predicting water use of winter wheat. 

A linear relationship was found between digital counts in band 6 and 
available soil moisture at 0 to 15 and 0 to 91 cm depths. Prediction of 
available soil moisture was improved by adding the ratio of band 4 to band 
5, which predicted LAI, as a second independent variable. Daily water use 
is the change in available soil moisture on successive days. 

Crop coefficient equations, based upon LAI, were developed for use in 
an ET model to predict daily water use of dryland and irrigated wheat. 
Predicted water use for the period March 21 to June 19, 1973, was within 
accepted accuracy for ET models. 



Y = 1.52-0.000834D 
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As only limited MSS data were collected because of excessive cloud 
cover, additional research is needed to verify and extrapolate the results 
of this study to fields with different soil types. 

2.2 Water Loss from Diseased and Healthy Wheat 

The patterns of water use in healthy and diseased wheat are important 
to both basic and applied research in terms of a clearer understanding of 
the epidemiology and of management practices that will reduce economic 
effects due to the disease. By accurately determining the water use of 
healthy and diseased crops and predicting yield reductions due to disease, 
we can provide the grower with information so that he can intelligently 
decide whether he should allow the diseased wheat to attain maturity or 
plow under the diseased wheat and replant a following crop. In western 
Kansas, wheat and other crops (i.e. soybean and sorghum) yields are largely 
determined by the availability of water; therefore, the earlier the grower 
decides to replant, the greater the conservation in soil moisture and a more 
likelihood of a successful crop. 

Procedure 

Five fields, Erichsen (commercial, healthy), Hartner (commercial, 
healthy), rust-inf ected, wheat streak mosaic virus-infected and control 
(healthy), were planted to Tritlcum aestivum L. cv. Scout in late September 
1972. Cropping descriptions for the five fields are given in Table 2.7. 

From the amount of water stored in the soil profile and the precipitation, 
the seasonal water use can be estimated. Soil moisture samples were deter- 
mined weekly (in some cases, inclement weather prevented weekly measurements) 
at 15-cm increments in a 150-cm profile. 



Table 2.7. Cropping descriptions of five wheat fields 


Field 

Soil Texture 

Seeding Rate 
Kg /ha 

Row Spacing 
cm 

Harvested 
plants per 
m 2 

Erichs en 

Silty clay loam 

96 

20.3 

856 

Hartner 

Silt loam 

101 

20.3 

872 

Rust- infected 

Silty clay loam 

84 

17.8 

959 

Wheat streak mosaic 

Silty clay loam 

84 

17.8 

1125 

Control 

Silty clay loam 

84 

17.8 

1144 
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Hourly estimates of the evapo transpiration rate (ET) were determined 

by the surface energy balance, 

ET “ (R - G) / (1 + 3) (2.23) 

n 

where R and G are the flux densities of net radiation and soil heat, res- 
n 

pectively. The Bowen ratio (JO is determined by y AT/Ae where y is the 
psychrometric constant and AT and Ae are the gradients of temperatures and 
vapor pressure above the canopy. Temperatures and water vapor pressures 
were determined with wet and dry bulb thermocouple psychrometers where the 
wick of the wet bulb was composed of a porous ceramic tube. Net radiation 
was determined with a hemispherical shielded radiometer located 3.0 meters 
above ground. Soil heat flux were determined with heat flux plates (5-cm 
depth) and calorimetrically. The sensors were scanned every 10 minutes and 
recorded on punch tape. The instruments and data acquisition system are 
described by Brun et al. 0-972). 

Stomatal resistance, R g , was determined with a diffusion porometer 
described by Kanemasu et al. (1969) and given by 

1/R s ■ 1/E ab + 1/R ad <2 ‘ 24) 

where R , and R , are the resistances of the abaxial and adaxial surfaces 
ab ad 

of the leaves. Leaf-water potential was estimated by pressure bomb technique 
(Barr, 1968). In both the stomatal resistance and leaf -water potential 
determinations, three upper canopy leaves were measured at near midday. 

Results and Discussion 

Figure 2.4 shows the water loss (determined by soil moisture) and the 
precipitation pattern. The amount of water loss from all five fields were 
quite similar; the water loss for the Hartner field being the largest. 
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This could be attributed to the lighter textured soil (silt loam) and the 
heavy rains (72 cm compared to the 30 year normal of 51 cm) which would 
be conducive to large drainage rates below the root zone. The average 
seasonal water loss for the five fields was 78.83 ± 5.89 cm. Nearly one- 
half of this amount being lost after dormancy (late March). By predicting 
the yield reductions 30 days in advance (e.g. predicting the economic loss 
from April 23 to 29 meteorological data'*’), a grower could conserve up to 
15 cm of water (assuming about 7-cm loss through evaporation) by plowing 
under his diseased crop; this soil water would then be available for the 
following crop. 

Figures 2.5, 2.6, and 2.7 show the hourly energy balance terms for 
May 4, May 19, and June 15, respectively. The rust infection and the 
control fields have nearly identical fluxes of net radiation and evapotrans- 
piration. However, the soil heat flux was approximately 50% greater for the 
control than for the rust- infected. On May 19 net radiation was greater 
for the control but the evapotranspiration was slightly greater for the rust- 
infected (Fig. 2.6). On June 15 (hard-dough stage. Table 2.8) the evapo- 
transpiration rates were similar for the two fields but the net radiation 
was greater for the control than for the rust-infected. 

Although there was a greater heat load on the control field, as indicated 
by the net radiation, evapotransporation was greater in the rust-infected 
field at time when rust spores were strongly evident (May 19) on the leaves. 

The greater ET rate would suggest a lower stomatal resistance on rust- 
infected leaves which was confirmed by data presented in Fig. 2.8 (c.f. May 19) 

1 

Appendix A. ERTS-1 Data Collection System Used to Predict Wheat Disease. 
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Figure 2.5 Trends in the hourly energy balance of rust and control (healthy) 
wheat for May 4, 1973. ET = evapotranspiration, Rn = net radiation, 
and G = soil heat flux. 
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Figure 2 0 6 Trends in the hourly energy balance of rust and control (healthy) 

wheat for May 19, 1973. et = evapo transpiration, Rn = net radiation 
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Figure ;2«7 Trends in the hourly energy balance of rust and control (healthy) 
wheat for June 15, 1973, et = evapotranspiration, Rn = net 
radiation, and G = soil heat flux. 
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Table 2.8. Physiological calendar for wheat 


Date 

Identifying 

Characteristics 

September 23 

planted 

April 21 

early joint 

May 1 

boot 

May 8 

heading 

May 18 

flowering 

May 30 

milk stage 

June 6 

soft dough 

June 15 

hard dough 


Table 2.9. Yield data for test 

wheat fields. 


Field 

Yield 

(bu/Acre) kg/ha 

1000 Kernal 
weight 
gm 

Erichsen 

(36. 3) 2566 

3 3.0 

Hartner 

(41.7) 2892 

31.1 

Rusted 

(28.4) 2616 

33.4 

Streak (WSMV) 

(36.7) 2621 

33.7 

Control (healthy) 

(40.0) 2729 

33.3 
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Lower stomatal resistances enhance diffusion of both water vapor (ET) and 
carbon dioxide (photosynthesis). The larger stomatal resistances of the 
wheat streak mosaic virus (WSMV) -infected plants indicate a loss of turgor 
pressure by the epidermal leaf tissue (decrease in leaf-water potential) « 
Fig* 2*9 shows the lower leaf— water potential of WSMV -infected compared to 
the control. Rust-infected leaves had a lower leaf-water potentials later 
in the growing season (May 29 to June 5), The unfavorable water balance of 
the diseased plants was reflected in the yields (Table 2.9), 



STOMATAL RESISTANCE (sec cm-') 



Figure 2.8 Daily trends in stomatal diffusion 

resistance of control, rust-infected and 
wheat streak mosaic virus (WSMV)- infected 
wheat leaves. 
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DAYS 1973 


Figure 2.9 Dally trends In the leaf-water potential 
of control, rust-infected and (WSMV)- 
infected wheat leaves. 
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3.0 DATA REDUCTION 
3.1 General Data Handling Method 

Initially we requested a standing order for 9" by 9" positive trans- 
parencies of MSS4 and MSS5 of each of the two fixed test areas (Riley 
county and the irrigated and nonirrigated fields in Finney county) and 
the high probability wheat streak mosaic virus (WSMV) areas of Kansas, 

Upon receipt of the transparencies they were examined for lack of clouds 
over the test areas and general suitability for further data reduction. 
Initially* several dates were examined and cataloged before determining 
which transparency would be used for further data processing. Later in 
the project when more frequent data were required and the anticipated 
arrival of the computer compatible tapes (CGT’s) was approximately three 
months after date of the overflight, it was decided to order tapes as soon 
as possible in order to minimize the delays. 

Upon receipt of these tapes, the transparencies were used to determine 
the desired test area and the tapes were then sent to the Remote Sensing 
Laboratory at the University of Kansas so that the desired area could be 
stripped off onto another tape and a gray-scale map was generated (Fig. 3.1). 
Since the second tape was organized in a band by band structure, a third 
tape was generated organized on a point basis. The computer maps were then 
put together for the location of the test areas by line and- cell numbers. 
Numerous methods were attempted for locating the. test fields (including the 
use of U-2 data); however, the best method was to use a clear piece of 
acetate over the entire area. This was laid over the computer maps and all 



I t M t t to D 0 

J 9 . * « 3 I I 

i , ; ; . ; . 


Reproduced from §|]j||| 

Figure 3.1 Computer Generated Gray Scale Map [best available copy, w 
8 33*9 H 5 | i'H!l| itsasf i^|3*afs| + + ♦ + } { t::; + ssfts« = x 


*♦*[?*+» .. 

* . . . $ =* 


. » .>***»« 

V e o o 

A 6 p 0 *, 0 ff . e o e — 8 o 

» • • *• ► , » ■ 

4 • • 

* -■ M » - • i , ■ r y 9 

i+5 : ■ ; . * n 

3*. « l'f*i p ? : j 9 t+to 


.99 ? * f *$tf<JMXX 8 X!tt()f 4 ! 9 s ! 9 8 5 * = ==33 = 333 + *** M** 3 5;|>.* ++ + *=^t = + =>! 

' 9 f 9 f 9 6 44 4 j 9 *3** J *4+=4$f 4*44 0 *?*=*?+ f l * 9 I * ? + 4 3 ss + +3 

{■*n v ! ! j9995*=4^^{ijn?n;9i*=sss+5+6s+^^" a t++nt5*+**3i=***? 
H + I v 9 9 tss:?*i 8 j ; 8 ; j j { s **+} 8 0 » . J . | « ! + s + l * U^+l 1 * 3 ! * M ♦*.» » !*==?= 
*4*+tf! | *<.3 53*1 } * + * }■*■** 5 | { ; f * | ft } 5 js^ssp HUM ‘ f l*ssttas4t* + * + 
? 3f"5 + + = f^MK4+ § liiiil'll |=^=3*8 S J+? 3 tW$XXXaX$ii + + + * ? 
; ; J 9 33SJ t » ? 3 ^m$w$**u 5 ft m * s+-?if4 a ! i +4W9WW *= ++4*+*+= 
i !?^s=*l ? UBHIOHBHIVidf'l 5?8i | +*34X*¥+ ? S = «3?XW$X1U = * M , H + + = = 
9N«^=* + *1 1 H8Sa@8a = pfMt + *s ===- = f^? t = « 9 1 t I 

9^^088091^5 f 9 8 !H*aW 

8i88$&==+*s ’xxiiixxfxftt+xxxf^jiHf^^t » • » t**avs&s 

9B88&Ktfc&:fc44?i*s5+ + ****$¥*4 !5]<B=sai48»fc*ij4 + 5%8$4*i l 1 8 «“8iS8$8 

in5llM*t»*+**s»**+*+i f ? 5 |*tigtiif * t4^i »?*#*+ 1 *=?=“?$ H 1 n+sXWHBQB 


8888SBB8S88 

$$@3011188®® 


9888808 

1180880 


***? s u*=*+5 = +H*4$W55?P88W*Ht t + 

^^3**j**=ti43| 53^« S Hlpp)!|g¥4e4?n » »M¥*8B!0BM|VVM4 
^£3*? 8 \ J 4 = 4XX4* = Xp4$+P®X$ 3 *» 5 *? in 5 5- 34»t&RVilM|8P9ltf 
£i s 3 j 9 i j + 4m44W+i+888f^ + + * +++ + S4*' t, *’*?^8BaB8ii8®i8i 
|j|**l*4s384^8^M4«t j&VWttfttKgff t*?9M*t*499i|BRI89!lSRB 
3**i&pSfSV888^f8iiSWli|8XXBlilBJ*lffH«»?X*5> + ?*)!l®»» M5* a *gB 

54i»lXXXXXXX?4 3 ^**5 5 t*W = ¥P = = S=l * t- 
Bl84;=ljia^X08B^iilM4+i|$881l|^lil^B8s#?45^l' HSH = S = 4W = **+ ? it 
@8fcS 33£=4XXH!il'590^l^«>?B0P* ; !.*l* U4+5 J4s+?«XBii|4f «'X!j!K 

B8X*t s + =?34XM8fc0^^XW§ + M » » ' H^ + 5 HXXS888B0s*,>HW 

B8|=t 1**8 isXiXltaXtf+S !*^4 + *9> 1 t 8 U?t?¥H¥VBVyXX¥9VKXXt 

@0|st*s*u |(tlii:KMI 5*£ssis»55?»XX«8X|18»*XXkHi«8X^s 
BBBp* t s u t i 1 1 ! » ? ? s !+? l i*+f 8 : l 5 f , ??RffR9iB9BRIR8RI|ii3Xlf4 


t0B8BgmR$B8^8 

taRRRRHHBBltiggRilB 


:ssfi9 9 s 


X* 9 9 t t j | j + = ** 
^4X^3 | *S j » I 9 f 


*X?4i 
I + *M« 




+ 4* g *4 j o t y* j 

!4$++V» ? * t I i « 
5^ a +];;;;»i? 

tu ; ^ ^ ^ 


!®^81itiH«8«S»U 


J- $ 'W »■ I* « ♦ ^ ^ w 21 m *¥ 9 * I 

4X$P?*li*«B9*|tH. 


3 4^SPi 


3^s*l 1 5 9 \ t4M . ;.H » » | = |ssf^s*3 = s + *X^S 

1X+ J ' t89«?«?»** ; » 9 I o 3 H I*l*t t f$ 

|f?n ,, 8 ¥ » ; i f »■», I M 5+41*,?* 

i q o , , ►,*».? f • • .• i^i [ i M 4 f * • | H 

4't o 5 j o.^sb ?#4 ( | i 9 9 I ^ V | t ft 

i.*- . ., > .*9 i Ki»'. ? 4 , » • f I • * * 5S4+1 4 = 

3jp«0 ?opo^6 pOO0e>« 55 

!?» J . , f f 1 ? ♦ , 99 J •* » M +, +?¥ 

POO ? PO e ?^P 5 Q .» * 4 b # 

fl . t 8 ? , M f fl ; f? ^ 8 

^ 9 |p p qo« 

p p 9 ^ ^ 0 00 0 r/ I | 

V 9 Jt J ^ Q 0 

■ 6 ^ H gpo a 

, S«Ui ?Ai d , t ni5?¥M***s¥H$B 

M ; ? M 

o 9 ; o J m ; . 

m 9 , S » * < u » 

9+434*1 5 n 1 + 8 H s =*t=!S8!j)¥ 3 + s ¥8«ia8 . 

* 9 d 4o**4*f?. t*8 , iS1? + + : '^9^^* = 3 408i8m§@ 

. ;.,;is+-»'i'f: u <r?4X@X+s*3«lS8llli^0 
4t4t*H£8?*tV8ii9 | 99?><$38 


+ * r*4^fpp 

j 1 1 1 9 T 9 n * 

t 9-4 9 *3 8 9 9 S 
* ! 9 

9MM |4QBS 




.“*5*^88888885 


! i?M * 


*. ■*■ v ^ 

§ 9 * 9 


< . .n. 'i j. 


* ■ ? 15 
• f I ^ 

” j o f t 


B J 3 9 

. P P §- ? k , 

• * tr *■ 1 i> V - 

Ap40 ?oOO0 


? . ? 9 


M i j j 9 ,, j , 

rf* (? » tf ’ 0 

1^9 0 5 M 0 ^ 

i ^ * 0 0 o Q ' 

* * ***** 


r ? n : : s 5 *♦» ?«X 
*» t o /: ? »» 1 1 ‘ # |4 
f 1 1 ? » • ; i r 3 1 s *| 

J 4 *} 444 * 4 * 4 | 4 ^ 
1 + 4} 44 9 M t f ? 9 | 


* 9. a G • 4 5 0 


m o t 

an n ■ 


R J Q 0 Q P 


M,Pf 


? 0 J , 

* 9 F 
:• f 0 

9 f . 


r/* j 

1 § • i 


0 



38 


pertinent identifiable points sketched on the acetate with a grease pencil* 
One of the most difficult tasks was the consistent identification of a test 
field on the maps. 

For the Riley county test fields the U-2 transparencies were projected 
over the computer maps in order to approximate the location of the fields* 
There were several identifying points that were initially thought to be 
identifiable on each date but were not* After the initial location had 
been verified by several people familiar with the test sites, the acetate 
sheet was used to identify the area on succeeding dates* 

For the Finney county test fields the most useful identifying marks for 
a large area were the county mile lines (roads every mile). Knowing their 
approximate spacing, a best fit to these mile lines was drawn on the computer 
map and again, by using several people familiar with the test area, an 
original layout of the test fields was made. It must be emphasized that 
one major check in this process must be the discussion of more than one 
trained person as to the location of a specific field and even then errors 
occur. The acetate overlay reduced the time in locating the fields on 
succeeding dates, especially in the case of the WSMV test area of Finney 
county with its over fifty fields. It should be noted that the N-S mile 
lines are much easier to locate than the E-W ones , 

After locating the fields, the comers of a straight line approximation 
to the field were recorded by line and cell number and key -punched. These 
were used to pick the data for a specific field from the point-oriented 
computer tapes* In order to check the locating of the field, two methods 
were used. First the points taken from a specific file were displayed as 
a n l H on a printout with all other points as a "O" which allowed for the 
detection of many small errors* This was made easier by the fact that each 
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file contained the data for a 64 by 64 block or area. However, the most 
useful method turned out to be the use of a Calcomp plot of the boundary 
of all fields on a specific data. These boundaries were generated using 
the identical data cards. Different dates could be overlaid and differences 
in size and location easily detected. An error of a hundred lines was 
detected in one case, 64 cells in another, and one line or cell in several 
cases. At this stage the data for each of the test areas were handled 
differently; however for each area one common program was run which indicated 
the mean, minimum, maximum, standard deviation, correlation coefficient, and 
histograms of each of the MSS bands, the difference of bands 6 and 4, 7 and 4, 
6 and 5, 7 and 5 (where band 7 had been multiplied by 2 to make the ranges of 
all band 128) and the ratios of bands 4 and 5, 4 and 6, 4 and 7, 5 and 6, and 
5 and 7. An example of this printout is shown in Figure 3.2. This particular 
printout will be discussed later in this chapter. 

3.2 Special Computer Programs Generated 

General Program 

A number of special computer programs were generated to assist in the 
location of the fields and in the data reduction for each field. The most 
frequently used program found the means, mlnimums, maximums, standard 
deviation, correlation coefficients and histogram. The input to this 
program consisted of the tape to be used, the inset or points to be dis- 
carded around each field, and the boundary data for each field (line and 
cell number for each comer). A program listing is given in Appendix E. 

A sample of the output is shown in Figure 3.2. This sample is for one 
entire file of data for tape number 1294-16521 which covered the WSMV 
region of Finney county on May 13, 1973. A total number of 114,688 points 
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are covered. It is interesting to note that only 3 quantities have a 
normalized cross correlation of less than 0*2 and these are MSS band 4 
and MSS band 6 which have a correlation of 0.20, MSS band 5 and MSS 
band 6 with a value of 0.04, and the ratios of MSS 4/5 and MSS 4/6 with a 
correlation of -0.12 which is probably due to the low correlation between 
bands 5 and 6. 

Calcomp CQNTUR Subroutine 

A Calcomp, contour-generating program was developed which is a very 
efficient program running in from 1/3 to 1/10 the time of other programs 
tested. The algorithm utilized 4 adjacent points requiring data which is 
equally spaced in the X direction and the Y direction but not necessarily 
the same spacing. This program was used to attempt to define the location 
of fields and specific targets for location; however this method did not 
meet with much success. This program can be easily modified to cross hatch 
a region between two or more given contours by modification of program to 
intercept the branch out when it is determined that the contour does not 
cross the rectangle. The computer subroutine, brief description and flow 
chart are given in Appendix F. 

Enhancing Category Variation 

A series of computer programs were written to implement a method of 
taking vector data and reducing it to a single number for ease of plotting. 
Data from unwanted categories would be near zero, and the data from the 
desired category would have maximized variation about some number. The proof 
of the algorithm is given in Appendix G. The computer programs are given 
in Appendix H. The first computes the eigenvectors and eigenvalues using a 
set of sample data. The second program computes the variation of data when 
the test eigenvectors are inputed. 
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3.3 CCT Data Used 

Table 3.1 gives a listing of the observation numbers of each of the 
tapes used for each of the test areas. The Riley county test field did not 
receive coverage from Hay 10, 1973 until July 4, 1973 due to cloud coverage 
even though it was in the overlap and was covered on two consecutive days 
each period. An almost three month period was missed in late fall for the 
irrigated and non— irrigated fields of Finney county due to cloud coverage. 

3.4 Data Reduction Recommendations for Other Investigators 

The following are a summary of our recommendations to investigators 
utilizing ERTS CCT's. 

1. Print a computer gray-scale map of the area. Be sure to include 
at least a 3 mile buffer all the way around a small test area to 
aid in location. If possible include large man-made objects such 
as airports, interstate highways, water bodies, etc. 

2. On the first attempt to locate the field have as many people as 
possible that are familiar with the test area go over the tenta- 
tive location. If possible, have them try to locate it without 
any information as to where others have placed it and then collect- 
ively discuss the location. 

3. After the initial location is defined, make an acetate overlay 
including the pertinent characteristics such as large dark areas, 
sharp boundaries, etc. MSS5 seems the best suited for this. Use 
this overlay to place the location of the test areas on later 
dates. Even then it is useful to have others duplicate your 
overlay. There will be a few point variation from date to date 
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Table 3*1* Computer compatible tapes used in the project . 

Satellite Pass 


Observation Number 

Date 


Comments 

Riley County Test Area 





1022-16391 

14 

Aug 

72 

Before Planting 

1058-16392 

19 

Sep 

72 


1237-16345 

17 

Mar 

73 


1256-16403 

5 

Apr 

73 


1274-16403 

23 

Apr 

73 


1291-16344 

10 

May 

73 


1328-16400 

15 

Jun 

73 

Clouds obscured field 

1346-16395 

4 

Jul 

73 

After Harvest 


Irrigated and Non-lrrigated Test Area Finney County 


1061-16564 

22 Sep 72 


1132-16514 

5 Dec 72 

On edge of tape 

1240-16523 

22 Mar 73 

1 

1295-16573 

14 May 73 


1312-16520 

1 Jun 73 


1330-16515 

19 Jun 73 

1 

1348-16514 

7 Jul 73 

After Harvest 1 


Wheat Streak Mosaic 

Virus Test Area Finney County 


1240-16523 

22 Mar 73 

1 

1294-16521 

13 May 73 


1313-16520 

1 Jun 73 


1330-16515 

19 Jun 73 

1 

1348-16514 

7 Jul 73 

After Harvest 1 


1-Used for both Finney County Test Areas 
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due to magnification and other optical errors hut a local region of 
possibly 4 to 10 square miles does not change appreciably. 

4. Provide cross checks on your locations whenever possible, such as 
test field boundaries overlaid on reflectance contours. If more 
than one test area is in the view, their relative positions should 


remain constant between dates. 
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4.0 EFFECT OF CROP GROWTH ON GRTS-1 MSS RESPONSE 

Determination of crop growth from spacecraft has received considerable 
attention by agriculturists , Such efforts have been brought about by the 
increasing awareness of shortages in food and fiber production. In a given 
area, agronomic crops develop and mature at somewhat predictable rates and 
abnormal growth patterns are exhibited when the photosynthetic process is 
interupted or reduced (e.g. disease, insects, nutrition, drought or flooding). 
Therefore, the monitoring of crop growth can provide valuable information 
on the prediction of production. 

Plants appear green to the human eye because the relatively large 
reflection in the green wavelength (500 nm) . This relatively high reflec- 
tance in the visible wavelengths can be attributed to the strong absorptance 
of blue (450 nm) and red (500 nm) wavelengths by plant pigments, namely, 
chlorophyll. More characteristic of healthy plant leaves than the low 
reflectance in the visible wavelengths (400-700 nm) is the high reflectance 
(about 50 percent) and transmittance (40 percent) to near infrared (700-1300 
nm) radiation. Because of this optical characteristic, leaves stacked on 
top of one another exhibit greater near infrared reflectance as the leaf 
layer increases; therefore, near infrared reflectance has important conse- 
quence in indicating differences in vegetation density (Allen and Richardson, 
1968). 

The reflected radiation stream from a crop canopy is composed of rays re- 
flected from the vegetation and from the soil surface. The surface, which domi- 
nates the scene reflectance depends upon percentage of crop cover. For many 
crops, especially row crops, there is a good correlation between percentage 

2 • 

Symposium on significant results obtained from ERTS-1. NASA SP-327. 

March 5-9, 1973. NASA-GSFC. 
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cover and leaf area index (leaf area to ground area) ; however as the percentage 
approaches 100 (about 85 %) large increases in leaf area index can occur with 
slight changes in percent cover. Since the spectral reflectance of soil 
differs from that of chlorophyll-containing tissue, vegetation density 
can be deduced from the signal strength of reflected rays in the visible 
and near infrared wavelengths. Kanemasu (Appendix B) suggests that ratio of 
the reflectances in wavelengths of MSS4 and MSS 5 is less than unity when 
soil exposure dominates the scene and tends to follow the leaf area index 
while Wiegand et al. (1973) states "the photosynthetic potential of green 
plants cannot be deduced directly from the photosynthetically active wave- 
lengths" (MSS4 and 5). However, the apparent inconsistency in the two 
studies may be due to the MSS4:MSS5 ratio being sensitive at low LAI (<2) 
while the MSS6 and 7 bands being more sensitive at high LAI (>2). At low 
LAI, the soil reflective properties in the visible wavelengths dominate while 
at high LAI the leaf reflective properties in the near infrared wavelengths 
dominate. 

METHODS AMD MATERIALS 

Four commercial wheat fields were used in study. Table 4.1 shows their 
location and cropping description. 

Plant samples for leaf area determination were collected at frequent 
intervals (usually 10 days) throughout the growing season. Samples in two 
of the fields (Hartner and Erichsen) were measured with an optical planimeter 
while the leaf area of samples from the other two fields were determined 
with an empirical equation using leaf width and length (equation 2.1). 
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Table 4.1. Location and cropping description of the Riley county fields 
(Hartner and Erichsen) and the irrigated and non-irrigated 
wheat fields in Finney county. 


Location 

Hartner - 39° 08' N, 96° 37* W 
Erichsen - 39° 07' N, 96° 35' W 
Irrigated - 38° 8.5’ N, 101° 4.9' W 
Nonirrigated - 38° 9.6' N } 101° 5.9' W 



25.4 
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Data reduction 

The greater the digital count in each band, the greater the radiance, 

-2 -1 

Maximum digital counts for MSS4, 5 and 6 is 127, The radiance (raw cm sr 
ym b is given by 

E =■ (digital counts) (maximum radiance) /127 (4.1) 

where the maximum radiance for 4, 5 and 6 are 24.8, 20.0 and 17.6, respec- 
tively. 

The radiance for WSS7 is 

E^ = (digital counts) (15.3)/63 (4.2) 

Hence, conversion factors (mw cm sr^ yirsT counts ) for digital counts in 
bands MSS4, 5, 6, and 7 to radiance are 0.19528, 0.15748, 0.13858, and 0.24286. 
The radiance, E, measured above earth T s atmosphere, has been affected by the 
optical properties of the atmosphere. 

Satellite reflectance is the ratio of the radiance, E ? to the irradiance 
(incoming radiation flux above atmosphere). In order to compare satellite 
reflectance with ground reflectance, the spectral transmission of the atmost- 
phere must be determined. The transmission of atmosphere becomes increasingly 
important to surfaces with low reflectances where errors can overwhelm low 
signal strength. 

Assuming equal atmospheric transmission for MSS4 and 5, we can describe 
the reflectance ratio of MSS4 and 5 as 



where and are the incident radiation fluxes in the wavelengths of MSS4 
and MSS5. Assuming (I,-/!^) equals 0.8696 (Smithsonian Meteorological Tables), 
equation (4.3) can be rewritten in terms of the conversion factors 



49 


^4 _ (Counts 4) (0.19528) 
Rj. (Counts 5) (0.15748) 


(0.8696) 


Counts 4 
Counts 5 


(1.078) 


(4.4) 


Similar relationships can be derived for the other combination of the 
reflectance ratios; however, caution should be used since the assumptions 
are not completely valid. The reflectance ratio of MSS4;MSS5 is slightly 
larger than the ratio of their digital counts. In this report, we will use 
digital count ratios and not attempt to estimate reflectance ratios. 


RESULTS 

Table 4.2 shows (a) the digital counts for MSS4, 5, 6, and 7, (b) the 
difference in digital counts, (c) the ratio of the digital counts, and 
(d) their standard deviations for the eight ERTS-1 observations on two Riley 
county wheat fields. Linear regression equations were calculated for the 
correlation of leaf area index (LAI) and percent cover (P.C.) with the 13 
digital parameters listed in Table 4.2. The regression equations for the 
pooled data (both fields) are summarized in Table 4.3 for the September to 
May observations. 

The lowest correlation coefficients were obtained directly from the 
MSS bands while the highest correlation coefficients were the difference 
between two bands, MSS6-5 and MSS (2 x 7)—5 ; however, band differences have 
high standard deviations. In general, the linear correlations were high 
whenever two bands were combined either by ratio or difference. Of the 
band ratios used, the poorest correlations were MSS4/6, and 4/7. It is 
significant that the MSS4/5, which consists of bands confined to the 
visible wavelengths, is one of the better indicators of crop growth. 

Table 4.4 shows the pixel correlation coefficients between the MSS 
bands for the Hartner field, A positive correlation existed between bands 



Table 4,2a. Mean digital counts and standard deviations for ERTS-l observations of Hartner field 


Observation No. 

NSS4 

MSS5 

MSS6 

MSS 7 

6-4 

(2x7)-4 

6-5 

(2x7)-5 

4/5 

4/5 

4/7 

5/6 

5/7 

1022-16391 

41.46 

44.82 

49.91 

24.14 

8.46 

6.82 

5.09 

3.46 

.95 

.83 

1.73 

.90 

1.88 

August 14, 1972 

(4.86) a 

(8.92) 

(2.64) 

(1,70) 

(3.62) 

(7.02) 

(7.57) 

(11.13) 

(.11) 

(.08) 

(.26) 

(.16) 

(.43) 

1058-16392 

34.74 

37.48 

37.39 

17.65 

2.65 

.57 

-.09 

-2.17 

.93 

.93 

1.98 

1.00 

2.14 

Sept. 19, 1972 

(3.83) 

(4.20) 

(2.43) 

(1.43) 

(3.59) 

(4.11) 

(4.14) 

(5.31) 

(.12) 

(.10) 

(.23) 

(.11) 

(.29) 

1076-16393 

.35.04 

36.38 

38.17 

19.17 

3.13 

3.29 

1.79 

1.96 

.97 

.92 

1.83 

.95 

1.90 

Oct. 7, 1972 

(3.30) 

(5.26) 

(2.99) 

(1.24) 

(2.25) 

(3.45) 

(3.30) 

(4.73) 

(.07) 

(.06) 

(.17) 

(.09) 

(.24) 

1237-16345 

30.41 

28.46 

36.55 

19.32 

6.14 

8.23 

8.09 

10.18 

1.07 

.83 

1.58 

.78 

1.48 

March 17, 1973 

(1.74) 

(2.30) 

(2.06) 

(1.39) 

(2.49) 

(3.34) 

(3.16) 

(3.92) 

(.07) 

(.06) 

(.15) 

(.08) 

(.18) 

1256-16403 

29.96 

27.18 

40.77 

22.14 

10.82 

14.32 

13.59 

17.09 

1.12 

.74 

1.37 

.67 

1.26 

April 5, 1973 

(2.75) 

(4.89) 

(2.73) 

(2.10) 

(4.28) 

(6.05) 

(6.48) 

(8.30) 

(0.11) 

(.10) 

(.225) 

(.15) 

(.35) 

1274-16403 

27.90 

22.80 

44.45 

25.10 

16.55 

22.30 

21.65 

27.4 

1.25 

.63 

1.12 

.51 

.92 

April 23, 1973 

(3.16) 

(5.04) 

(1.40) 

(1.59) 

(2.80) 

(5.70) 

(4.66) 

(7.53) 

(.12) 

(.06) 

(.21) 

(.11) 

(.28) 

1291-16344 

30.17 

22.11 

48.39 

27.89 

18.72 

25.61 

26.78 

33.67 

1.37 

.62 

1.08 

.45 

.79 

May 10, 1973 

(1.69) 

(1.91) 

(2.19) 

(1.08) 

(2.11) 

(2.73) 

(2.29) 

(2.79) 

(.06) 

(.03) 

(.07) 

(.04) 

(.07) 

1346-16395 

37.0 

39.75 

46.25 

22.90 

9.25 

8.80 

6.50 

6.05 

.93 

.80 

1.62 

.87 

1.76 

July 4, 1973 

b 

(3.46) 

(4.87) 

(2.47) 

— 

— 

(5.69) 

(6.30) 

— 

— 

— 

(.11) 

(.24) 


a (Standard deviation) 
^Discontinuity in telemetry 



Table 4* 2b. Mean digital counts and standard deviations for ERTS-1 observations of Erichsen field 


Observation No. 

MSS4 

MSS5 

MSS6 

MSS7 

6-4 

(2x7)-4 

6-5 

(2x7)-5 

4/5 

4/6 

4/7 

5/6 

5/7 

1022-16391 

26.06 

21.13 

27.88 

13.75 

1.81 

1.44 

6,75 

6.38 

1.24 

1.01 

2,13 

.82 

1.75 

August 14 , 1972 

(1.98) a 

(1.36) 

(8.98) 

(5.70) 

(7.20) 

(9.65) 

(8.47) 

(10,91) 

(.08) 

(.24) 

(.62) 

(.22) 

(-57) 

1058-16392 

30.54 

28.21 

29.33 

13.83 

-1.21 

-2.88 

1.125 

- .54 

1.09 

1.05 

2.24 

.97 

2.07 

Sept* 19, 1972 

(2,11) 

(3.08) 

(2.85) 

(1.86) 

(3.04) 

(3.79) 

(3.51) 

(4.23) 

(.07) 

(.11) 

(.31) 

(.12) 

(.30) 

1076-16393 

22.40 

19.18 

18.69 

9.41 

-3.69 

-3.56 

- .49 

- .36 

1.18 

1.23 

2.48 

1.04 

2.10 

Oct. 7, 1972 

(1.58) 

(2.51) 

(3.54) 

(2.38) 

(2.48) 

(3.80) 

(2.09) 

(3.27) 

(.10) 

(.17) 

(.46) 

(.11) 

(.30) 

1256-16403 

23.57 

17.93 

26.14 

14.64 

2.57 

5.71 

8.21 

11.36 

1.33 

.92 

1.65 

.70 

1.27 

April 5, 1973 

(2.21) 

(3.03) 

(2.96) 

(2.34) 

(4.33) 

(5.80) 

(5.24) 

(6.69) 

(.12) 

(.16) 

(.33) 

(.17) 

(.33) 

1274-16403 

25.77 

17.82 

41.59 

24.59 

15.82 

23.41 

23.77 

31.35 

1.48 

.64 

1.10 

.45 

.78 

April 23, 1973 

(2.31) 

(3.97) 

(6.28) 

(4.78) 

(7.82) 

(10.92) 

(9.68) 

(12.77) 

(.18) 

(.15) 

(.34) 

(.18) 

(.35) 

1291-16344 

27.82 

18.12 

50.44 

29.32 

22.62 

30.82 

32.32 

40.53 

1.55 

.57 

.99 

.37 

.66 

May 10, 1973 

(1.09) 

(1.74) 

(7.78) 

(5.74) 

(7.57) 

(11.28) 

(8.95) 

(12.68) 

(.15) 

(.10) 

(.23) ( 

(.10) 

(.21) 

1346-16395 

49.0 

59.47 

62.00 

30.33 

13.00 

11.66 

2.53 

1.20 

.82 

.79 

1.62 

.96 

1.96 

July 4, 1972 

b 

(2.45) 

(2.73) 

(1.84) 

— 

■ 

(1.55) 

(2.70) 

— 

— 

— 

(.02) 

(.09) 


3 

(Standard deviation) 
b 

Discontinuity in telemetry 



Table 4.3a. Linear regression equations of leaf area 
index (LAI) and digital counts from MSS 
4, 5, 6, and 7 taken from 6 ERTS-1 ob- 
servations of Erichsen and Hartner fields. 


Correlation 

Linear Regression Equation Coefficient 

LAI = -0.081 (MSS 4) + 3.333 «.308 

LAI - -0.101 (MSS 5) + 3.414 -.704 

LAI - 0.072 (MSS 6) - 1.824 .734 

LAI = 0.126 (MSS 7) - 1.717 .838 

LAI - 0.095 (MSS 6-4) - .001 .910 

LAI - 0.071 [MSS (2 x 7) -4] - .028 .937 

LAI = 0.074 (MSS 6-5) - .128 .958 

LAI - 0.058 [MSS (2 x 7)-5) - .107 .962 

LAI = 4.148 (MSS 4/5) - 4.195 .918 

LAI = -3.946 (MSS 4/6) + 4.095 -.815 

LAI = -1,696 (MSS 4/7) + 3.472 r^.864 

LAI = -3.665 (MSS 5/6) + 3.400 -.945 

LAI - -1.565 (MSS 5/7) + 2.952 -.920 
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Table 4.3b. Linear regression equations of percent 
cover (PC) and digital counts from MSS 
4, 5, 6, and 7 taken from 6 ERTS-1 
observations of Erichsen and Hartner 
fields. 


Correlation 

Linear Regression Equation Coefficient 

PC * -3.457 (MSS 4) + 137.492 -.410 

PC - -3.142 (MSS 5) + 115.858 -.729 

PC « 2.579 (MSS 6) - 63.474 .717 

PC - 4.573 (MSS 7) - 59.997 .843 

PC - 3.364 (MSS 6-4) + 4.214 .915 

PC - 2.465 I MSS (2 x 7)-4^ + 4.403 .942 

PC = 2.536 (MSS 6-5) + 2.100 .967 

PC - 1.953 IMSS (2 x 7)-5] + 3.318 .968 

PC = 131.342 (MSS 4/5) - 123.603 .892 

PC = -146.141 (MSS 4/6) + 152.736 -.869 

PC = -61.194 (MSS 4/7) + 128.859 -.900 

PC = -123.929 (MSS 5/6) + 122.025 -.963 

PC = -52.937 (MSS 5/7) + 107.027 -.952 



Table 4.4. Correlation 
field. 

coefficients 

between the various MSS 

bands for the 

Hartner 

Date 

4 vs 

5 4 vs 6 

4 vs 7 

5 vs 6 

5 va 7 

6 vs 7 

August 14 

.90 

.65 

-.41 

.59 

-.51 

.01 

September 19 

.41 

.39 

.26 

.30 

-.10 

, 66 

October 7 

.88 

.72 

.30 

.78 

.42 

.69 

March 17 

.60 

.14 

-.03 

-.04 

-.17 

.73 

April 5 

.87 

-.21 

-.47 

-.38 

-.64 

.75 

April 23 

.91 

.44 

-.59 

-.38 

-.63 

.23 

May 10 

.79 

.41 

.01 

.36 

.06 

.68 

July 4 

.13 

.09 

.09 

.09 

-.09 

.86 
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4 and 5 and between 6 and 7. The lowest positive correlation existed between 
a visible band (4 or 5) and a near- infrared band (6 or 7). This relationship 
results because vegetation absorbs strongly in the visible and reflects 
strongly in the near- infrared wavelengths. 

Radiometric response . Fig, 4.1 shows the scene radiance from a single wheat 
field at various growth stages. In general, as the plants develop, radiance 
in bands 4 and 5 decrease while bands 6 and 7 increase. The higher radiance 
at LAI ■ 0 than at LAI = .14 may be due to soil moisture differences (Appendix 
B). The stubble field (post harvest) shows a high reflectance in the visible 
(MSS 4 and 5) and in the near infrared wavelengths (MSS 6 and 7). 

Finney County Fields . In the Finney county area, leaf area indices were 
determined periodically on an irrigated and a nonirrigated wheat field 
(c.f. section 2.1). The Riley county (Hartner and Erichsen fields) and 
Finney county (irrigated and non-irrigated fields) data were pooled and 
linear regression equations were determined for the various band ratios 
(Table 4.5). The MSS4 to MSS5 ratio gave the highest correlation coefficient 
with the MSS 5 to MSS6 ratio having a slightly lower correlation coefficient. 
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Fig. 4.1. Radiometric response of ERTS-1 bands for 

Hartner wheat field during 1973 growing 
season. 



Table 4.5. Linear regression equations of pooled 
LAI and MSS digital count data from 
Hartner and Erichsen fields (Riley 
county) and irrigated and non-irrigated 
fields (Finney county) . 


Linear Regression Equation 

Correlation 

Coefficient 

LAI - 

3.305 (MSS4/5) - 3.089 

.92 

LAI = 

-3.395 (MSS4/6) + 3.275 

-.72 

LAI = 

-1.384 (MSS4/7) + 2.650 

-.73 

LAI = 

-3.034 (MSS5/6) + 2.804 

-.90 

LAI = 

-1.307 (MSS5/7) + 2.402 

-.87 
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5.0 DETECTION OF DISEASE SEVERITY AND ECONOMIC LOSS 

ERTS-1 imagery was processed to determine the satellite’s usefulness 
in the early detection and estimation of wheat disease severity and economic 
loss. The experimental site was a 450 square mile area of Finney and Gray 
counties which contained healthy and wheat streak mosaic virus (WSMV)- 
infected wheat fields. A detailed map of this area was prepared and the 
crop and its condition was determined in every field. The following severity 
ratings were used: 

0 = Healthy, No WSMV infection (very green) 

1 = Trace WSMV infection (few yellow plants) 

2 = Moderate WSMV infection (whole field slightly yellow) 

3 = Severe WSMV infection (whole field very yellow) 

With farmers permission, four random samples (a sample = 16 square 
foot plot) were harvested from each of 54 fields. Samples were threshed 
for yield determination, and four yield groups were assigned: 

A = 13.8 - 20.00 Bushels/Acre 
B = 20.1 - 25.0 Bushels /Acre 
C = 25.1 - 30.0 Bushels/Acre 
D = 30.1 - 41.6 Bushels /Acre 

Grey scale maps of the site were prepared for MSS Band 5 for each clear 
date (March 20, May 13, May 31 and June 18, 1973). Test fields were readily 
identified using center-pivot irrigators, airports and peculiar field 
shapes for registration. Digital data for the four MSS bands on each data 
were processed manually or by machine for each field. To minimize border 
effects a boundary of one data point around the edge of each field was 


discarded. 
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Means were determined for digital counts on all MSS bands and ratios 
of the bands. These means were tabulated relative to yield (Table 5.1) 
and disease severity (Table 5.2). Correlation coefficients for MSS digital 
counts versus yield and MSS digital counts versus severity were determined. 
None of these relationships was significant at the 1% level, but 11 of 72 
relationships were significant at the 5 % level. Of these, eight occurred 
in the 5/31/73 data. On this date a negative correlation existed between 
yield and MSS bands 4 and 5, with a positive correlation between yield and 
the MSS 4/5 band ratio (Table 5.3, Fig. 5.1-C). A significant negative 
correlation existed (barely) on 5/31/73 between yield and the MSS 4/6 and 
4/7 band ratios (Table 5.3). MSS band 6 was correlated with yield on 
3/20/73 (Table 5.3). MSS band 7 never showed a correlation with yield at 
the 5% level (Table 5.3). 

Similar correlations occurred between MSS digital counts and disease 
severity. For 5/31/73 a positive correlation existed hetween severity and 
MSS band 4 and for the MSS 4/6 and 4/7 band ratios (Table 5.4, Fig. 5.2-C). 
A higher positive correlation existed on 6/18/73 between MSS hands 4 and 5 
and severity (Table 5.4, Fig. 5,2-D). There was no significant correlation 
between MSS bands 6 or 7 and severity (Table 5.4). 

Significant correlations between MSS digital counts and wheat yields 
or disease severity were demonstrated only for the 5/31/73 and 6/18/73 
data. The negative correlation between yield and digital counts for MSS 
bands 4 and 5 on 5/31/73 (Fig. 5.1-C) is reasonable and probably resulted 
from premature coloration of the crop, greater reflectance by the soil due 
to thinning of the crop or both. WSMV prematurely colors and thins the 
crop. Positive correlations between severity and digital counts for MSS 
bands 4 and 5 on 5/31/73 and 5/18/73 (Fig. 5.2-C,D) are also plausible for 


the same reasons. 



60 


We have correlated ERTS-1 imagery with ground truth for both wheat 
yield and disease severity with significant correlation obtained at the 5% 
level. Also, in both cases the effects of the disease are being detected 
near the end of the crop season rather than the disease per se being detected 
early in the season. It is reasonable to assume that a disease sufficiently 
severe to reduce yields by over 50% and readily detectable by eye should 
have been more readily detected and quantified by ERTS-1. However, WSMV 
was most obvious from very late March to May 1 when no ERTS~1 imagery was 
available due to inclement weather. Therefore, although we report some 
positive results on the quantitative effects of the disease, we were unable 
to adequately test the ERTS"1 system for early detection and estimation 
of severity and yield reduction in wheat due to wheat streak mosaic virus. 
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Table 5.1, Means of MSS digital counts in relation to wheat yields. 


MSS Band MSS Band Ratio 


Yield 


Date 

Group 

Yield 

4 

5 

6 

7 

4/5 

4/6 

4/7 

5/6 

5/7 

3/20/73 

A 

17.2 

25*8 

22.4 

29.7 

16.2 

1.2 

0*9 

1.6 

0.8 

1.4 


Std. Dev. 

2.4 

0.7 

1.1 

2*8 

1.9 

0.0 

0*1 

0.2 

0.1 

0.2 


B ’ 

22.8 

27.0 

24.3 

29.8 

15.8 

1*1 

0.9 

1.7 

0*8 

1*6 


Std. Dev. 

1.4 

1.7 

2.9 

1.5 

1.2 

0.1 

0.1 

0*2 

0*1 

0.3 


C 

28.3 

26.4 

22.8 

31.9 

17.3 

1*2 

0*8 

1*5 

0*7 

1.3 


Std. Dev. 

1.3 

1.5 

3.1 

3.6 

2.6 

0.1 

0.1 

0*2 

0*1 

0.3 


D 

35.7 

26.8 

23.6 

31.9 

17.5 

1.2 

0.8 

1.6 

0*8 

1.4 


Std. Dev. 

3.6 

2.6 

4.4 

3*3 

2.1 

0.1 

0.1 

0*2 

0.1 

0.3 

5/13/73 

A 

17.5 

34.4 

29.7 

49,7 

27.3 

1.2 

0*7 

1.3 

0.6 

1.1 


Std. Dev. 

2.5 

2.7 

3.6 

6.8 

4*6 


0.1 

0.2 

0*1 

0.2 


B 

23.3 

35.5 

30.6 

54.8 

30.7 

1.2 

0.7 

1.2 

0*6 

1*1 


Std. Dev. 

1.0 

3.8 

5.7 

7.4 

5*2 

0*1 

0.1 


0*1 

0.3 


C 

28.0 

34.7 

29.3 

54.0 

31*1 

1.2 

0*7 

1*2 

0.6 

1*0 


Std. Dev. 

1*3 

2.8 

4*6 

6.4 

5*1 


0.1 

0*2 

0*1 

0.3 


D 

36.3 

34.7 

29.7 

52.0 

29.4 

1.2 

0.7 

1.3 

0.6 

1.1 


Std. Dev. 

3.9 

6.4 

8.7 

12.3 

8.3 

0*1 

0*2 

0.6 

0.2 

0*5 

5/31/73 

A 

17.1 

37.0 

35.6 

53.0 

29.3 

1.1 

0*7 

1.3 

0*7 

1.2 


Std. Dev. 

2.2 

2.3 

3.3 

5.4 

3.2 

0*0 

0*1 

0.1 

0.1 

0.2 


B 

22.9 

35.4 

33.2 

53.0 

30*1 

1.1 

0.7 

1.2 

0.6 

1.1 


Std. Dev. 

1.3 

3.5 

6.1 

2.7 

2*8 

0.1 

0.1 

0.2 

0*1 

0.3 


C 

27,6 

32.6 

29.2 

51.6 

30*6 

1*2 

0.6 

1.1 

0.6 

1.0 


Std* Dev. 

1.4 

3.9 

6.1 

3.1 

2.9 

0*1 

0.1 

0.2 

0.1 

0,2 


D 

36.9 

32.6 

28.9 

52.6 

31.3 

1*2 

0*6 

1.1 

0.6 

1.0 


Std. Dev. 

3.6 

3.9 

7.2 

2.6 

3*3 

0*1 

0*1 

0.3 

0.2 

0.4 

6/18/73 

A 

17.6 

42.5 

48.1 

53.1 

27.3 

0.9 

0.8 

1.6 

0.9 

1.8 


Std. Dev. 

1.9 

4.7 

6.3 

5.2 

2.7 

0*0 

0.1 

0.2 

0.1 

0.2 


B 

22.9 

41.4 

45.8 

53.3 

28.0 

0.9 

0.8 

1.5 

0.9 

1.7 


Std. Dev. 

1.3 

3.2 

2.8 

4*1 

2.4 

0.0 

0*1 

0*1 

0*1 

0*2 


C 

27.3 

38.5 

42.5 

52*3 

28.3 

0.9 

0.7 

1.4 

0.8 

l t 5 


Std, Dev. 

1.4 

3.3 

5.4 

3.0 

1.1 

0.0 

0.0 

0.1 

0,1 

0*2 


D 

36.7 

39.6 

43.8 

51*6 

27.5 

0.9 

0.8 

1.5 

0.9 

1.6 


Std. Dev. 

3.5 

3.0 

4.3 

2*5 

1.3 

0*0 

0.0 

0.1 

0.1 

0.2 
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Table 5,2* Means of MSS digital counts in relation to disease severity 


MSS Band 


Severity 
Date Group 

3/20/73 0 

Std. Dev, 


Std. Dev. 


Std. Dev. 


Std. Dev. 


5/13/73 0 

Std. Dev. 


Std. Dev. 


Std. Dev. 


Std. Dev. 


5/31/73 0 

Std. Dev. 


Std. Dev. 


Std. Dev. 


Std. Dev. 


6/18/73 0 

Std. Dev. 


Std. Dev. 


Std. Dev. 


Yield 

4 

5 

6 

26.2 

27.1 

24.4 

30,1 

6.8 

1.9 

3.4 

3. J 

32.4 

26,2 

22.6 

32, i 

4.9 

2.1 

3.6 

4 m i 

27.5 

27.0 

24.0 

30, ] 

7.1 

2.2 

3.7 

l.i 

20.8 

26.0 

22.6 

31. ( 

6.3 

0.8 

1.4 

1J 

28.2 

35.6 

31.6 

51.! 

6.2 

3.8 

6.3 

5.: 

31.8 

31.6 

25.3 

48.1 

5.2 

3.7 

4.8 

13.: 

28.0 

36.1 

30.9 

57.: 

7.9 

4.6 

6.5 

7.: 

22.7 

35.8 

30.8 

53. 

6.2 

4.2 

5.1 

7. 

26, 6 

32.5 

29.0 

52. 

7,1 

4.4 

6.8 

3. 

32.2 

33.8 

31.2 

51. 

5.1 

4.6 

8.1 

3. 

27.5 

35.1 

32,2 

53. 

7.8 

3.6 

6.4 

2. 

20.0 

35.8 

34,0 

51. 

3,8 

1.4 

2.3 

3. 

26.6 

38.8 

42.7 

51. 

7.1 

3.7 

5.6 

3. 

32.2 

39.7 

44.1 

52, 

5.1 

3.1 

4.3 

2. 

27.5 

40.8 

45.3 

52. 

7.8 

2.8 

2.9 

3. 

22.1 

42.5 

47.6 

53. 

5.7 

4.7 

6.0 

4. 



MSS Band Ratio 

4/5 4/6 4/7 5/6 5/7 

.8 1.5 

.1 0.3 

.8 1.5 0.7 1.3 

.1 0.2 0.1 0.3 


.9 1.7 

.1 0.2 

.8 1.5 

.0 0.1 


.8 1.5 


0.9 0.8 1.4 0.8 

0.0 0.1 0.2 0.1 

0.9 0.8 1.4 0.8 

0.0 0.0 0.1 0.1 

0.9 0.8 1.5 0.9 

0.0 0.1 0.1 0.1 

0.9 0.8 1.5 0.9 

0.0 0.0 0.1 0.1 
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Table 5.3. 

Correlation 

coefficients - MSS 

digital 

counts 

vs, wheat yields 





MSS 

Band 



MSS 

Band Ratio 


Date 

4 

5 

6 

7 

4/5 

. .4/6 

4/7 

5/6 

5/7 

3/20/73 

.18 

.14 

.32* 

.30 

-.02 

-.18 

-.18 

-.09 

-.10 

5/13/73 

-.02 

-.01 

-.07 

-.03 

.09 

.13 

.13 

.08 

.09 

5/31/73 

-.39* 

-.34* 

-.06 

.21 

.34* 

-.31* 

-.31* 

-.30 

-.29 

6/18/73 

-.23 

-.24 

-.16 

-.01 

.12 

-.17 

-.21 

-.20 

-.22 


*Significant at 5% level. 

Correlation coefficients with a minimum of 40 degrees of freedom; 1% = .393 

5% = .304 
10% = .257 


Table 5.4. Correlation coefficients - MSS digital counts vs. disease severity 


Date 


MSS Band 



MSS Band Ratio 


4 

5 

6 

7 

4/5 

4/6 

4/7 

5/6 

5/7 

3/20/73 

-.13 

-.10 

.01 

-.01 

.04 

-.13 

-.10 

-.12 

-.10 

5/13/73 

.10 

.02 

.20 

.16 

.04 

-.13 

-.11 

-.14 

-.12 

5/31/73 

.32* 

.27 

-.05 

-.28 

-.29 

.32* 

.32* 

.27 

.28 

6/18/73 

.35* 

.35* 

.22 

.03 

-.18 

.25 

.27 

.29 

.28 


*Slgnificant at 5% level. 
Correlation coefficients as above. 
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7.0 CONCLUSIONS 

A computer-based model of estimating evapotranspiration for irrigation 
scheduling requires a crop coefficient curve. The leaf area index (LAI) of 
wheat appears to provide a reasonable estimate of a wheat-crop coefficient 
and can be deduced from MSS digital data. The ratio of any two MSS bands 
(digital counts) are linearly correlated with LAI (r > .70); the ratios of 
MSS4:MSS5 and MSS5:MSS6 appear to best simulate LAI (r > .90), In addition, 
these linear regression equations are useful for estimating LAI of wheat on 
fields other than which the equation were derived. 

The soil moisture in the roo zone of wheat was estimated with reason- 
able success from MSS7 and MSS4:MSS5; however, the relationship is not unique 
and would depend upon soil type. This preliminary work needs further investi- 
gation. 

The early detection of wheat streak mosaic virus using the available 
ERTS-1 imagery for 1973 was not possible. At the time of greatest visual 
difference between healthy and disease wheat, no ERTS-1 imagery was 
available because of cloud cover. However, we attempted a later detection 
of disease severity using the May 13, 1973 Imagery but no significant cor- 
relation of disease severity to band digital count was found. The May 31, 

1973 imagery showed a higher degree of significance (5% level) of digital 
counts versus severity and digital counts versus yield. This relationship 
was attributed to the premature yellowing and thinning of diseased wheat. 

Our use of data collection platforms (DCP) to predict disease 
severity was successful and appears to have potential applications. The 
implementation of the DCP did not create a major maintenance problem. 
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8.0 FUTURE RESEARCH NEEDS 

Currently we are at the stage of development where it is feasible to 
monitor large agricultural areas and Identify the crop type and stage 
of development. However, further study is required for the detection and 
evaluation of disease and water stress. These physiological stresses are 
of major importance in their effect on crop yield. As with most agricultural 
programs, the timeliness of the observations is extremely critical. 
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ERTS-1 DATA COLLECTION SYSTEMS USED TO PREDICT 

WHEAT DISEASE SEVERITIES 1 

E. T. Kanemasu, H. Schimmelpfennig, E. Chin Choy, 

2 

M. G. Eversmeyer and D. Lenhert 


ABSTRACT 

The feasibility of using the data collection system on Earth Technology 
Satellite-1 to predict wheat leaf rust severity and resulting yield loss was 
tested. 

Ground-based data-collection platforms (DCPs) , placed in two commercial 
wheat fields in Riley County, Kansas, transmitted to the satellite such 
meteorological information as maximum and minimum temperature, relative 
humidity and hours of free moisture. Meteorological data received from 
the two DCPs from April 23 to 29 were used to estimate the disease progress 
curve. Values from the curve were used to predict the percentage decrease 
in wheat yields resulting from leaf rust. Actual decrease in yields was 
obtained by applying a zinc and maneb spray (5.6 kg /ha) to control leaf 
rust, then comparing yields of the controlled (healthy) and the noncontrolled 
(rusted) areas. In each field a 9% decrease in yield was predicted by the 
DCP-derived data; actual decreases were 12% and 9%. 


Contribution No. 1387, Agronomy Department, Eyapotranspiration Laboratory, 
and Contribution No. 595, Department of Plant Pathology, Kansas Agricul- 
tural Experiment Station, Kansas State University in cooperation with the 
United States Department of Agriculture. National Aeronautics and Space 
Administration provided partial support for this research. Date received 


Assistant Professor of Microclimatology, Electronic Technician, Research 
Associate, Research Plant Pathologist, USDA, ARS, NCR, and Associate 
Electrical Engineer, Department of Electrical Engineering, Engineering 
Experiment Station, Kansas State University, Manhattan, Kansas 66506. 



INTRODUCTION 


Epidemiological investigations have shown that the severity of wheat 
leaf rust (Puceinia recondita Rob. ex Desm f . sp. tritici) and subsequent 
loss in yield can be predicted (Eversmeyer and Burleigh, 1970; Burleigh 
et al., 1972a; Burleigh et al. , 1972b). Such predictions would be most 
important for determining curative measures to reduce economic loss. 

In the above mentioned investigations, stepwise multiple regression 

techniques were used to identify biological and meteorological variables 

useful in explaining variation in wheat leaf rust severities 7, 14, 21, and 

30 days after the date of prediction (DP) and the relationship between 

those predicted severities and yield loss. Equations in the form 

Y = K. + b.x, . + ... + b x . were formulated and tested. Variables which 
i i 1 li n ni 

they reported to be most significant in the successful prediction of wheat 
leaf rust development were: leaf rust severity on DP, growth stage of 

wheat on the date predicted, average hours of free moisture during seven 
days prior to DP, number of days or precipitation greater than or equal to 
0.25 mm during seven days prior to DP, a fungal growth function, and fungal 
infection function. The equations predicted leaf rust severity in test 
plots within ±1, 3, and 12%, 14, 21, and 30 days in advance, respectively. 
They studied the relationship between leaf rust severity at several wheat 
growth stages and yield loss and constructed general equations to predict 
percent loss. 

Successful prediction of disease losses for large remote areas would 
require continuous gathering of meteorological and biological data on 
widely separated fields, which using routine instrumentation would require 
an enormous maintenance capability. The recent launching (July 23, 1972) 



of Earth Resources Technology Satellite-1 CERTS-l) has permitted the use 

of spacecraft to collect data from ground-based transmitters placed in 

remote areas; in turn the satellite can retransmit data to one of three 

prime receiving stations: Goldstone, California; NASA Test and Training 

Facility; and Fairbanks, Alaska. (Detailed information can be found in 

3 

the Data Users Handbook) . 

The satellite’s data collecting capability offers a unique opportunity 
to test and evaluate the use of information gathered by data collection 
platforms to predict epidemics of wheat leaf rust. 

MATERIALS AND METHODS 

Two data-collection platforms (DCP) furnished by NASA as part of an 
ERTS-1 experiment (site 1, 39 Q 07 T N, 96 Q 35’W; site 2, 39°08’N, 96°35'W) — 
were located in two Riley County (Kansas) commercial fields (40 acres each) 
of wheat (Triticum aes tivum L. cv. Scout) . Site 1 and site 2 fields were 
a silty clay loam and silt loam, respectively. Roth fields were planted 
to wheat in late September in 20-cm rows. The DCPs were installed by 
December 1, 1972 and were operational until July 23, 1973, 

The data collection platform (DCP) is an automatic, data-relay ter- 
minal that accepts 8 channels of either analog or digital data from user- 
furnished electronics (sensor interfaces) . Every 3 minutes the DCP 
interrogates the 8 input channels and transmits the data, ,rggardless of 
the satellite’s position. The satellite passed close enough to the Riley 
County DCPs to receive the transmission with a ± 1 hour period at 1030 

3 

Available through General Electric, Space Division, Valley Forge Space 
Center, P. 0. Box 8555, Philadelphia, Penn. 19101. 
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and 2230 local time. Because many of our sensors required interrogation 
at other times, we designed a sensor interface that interrogates sensors 

at the proper time, then digitally stores the data for transmitting later. 

4 

(A detailed description of the interface appears in a NASA report) . 

The power supply (storage batteries) was enclosed in a box separate 
from the rest of the DCP to preyent corrosion by acid fumes. The trans- 
mitter and interface were enclosed in a double wooden box. 

Table 1 gives sensor data transmitted and received by the satellite 
during the two periods of each day that the satellite was within DCP range. 
Information on channels 1 through 4 was obtained primarily for input into 
the disease prediction equations; the visible and near infrared reflectance 
data (channels 5-8) were used to analyze crop growth. 

Relative humidity was measured by a sulfonated polystyrene sensor 
(PCRC-11, Phys-Chemical Research Corp.), and soil moisture was estimated 
by gypsum soil moisture blocks (CEL-WFD, Beckman Instruments). The signal 
conditioning consisted of an AC ohmmeter and a logarithmic amplifier for 
linearization. Free moisture was detected by measuring the AC resistance 
change of a bifilar array exposed to the atmosphere. When wetted, the 
electrical resistance of the array decreases and the resistance of the array 
was determined with a level-detecting AC ohmmeter. Air temperatures 
(maximum, minimum, and instantaneous) were determined by thermalinear 
thermistors (YSI 700, Yellow Spring Instruments). Maximum and minimum 
temperatures stored in the interface memory were compared with the current 
temperature every 3 minutes and updated. Visible (590 to 720 nm) and 

4 

Report No. 2263-3, Kansas Environmental and Resource Study; A Great 
Plains Model. January, 1973. NASA. 
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near infrared (730-1000 run) radiation streams were measured using silicon 

photocells filtered on their respective wavelengths ranges. (Details on 

4 

sensors and signal conditioning can be found in a NASA report) . 

The relative humidity and temperature sensors were located in small 
ventilated weather shelters. The weather shelters were positioned at least 
35 m from the edge of the field. Early in the season, the sensors were 
maintained at a height of about 30 cm above the soil surface but, as the 
plants developed, the shelters were raised to keep the sensors near the 
top of the canopy. The free moisture sensors were maintained at mid-canopy 
height. The photocells were located on stands approximately 2.5 m above 
the soil surface. 

Beginning April 25, weekly applications of a zinc and maneb spray 

2 

(5.6 kg/ha) were made on four randomly selected 2.2 m areas (healthy) 

near each DCP for control of wheat leaf rust. Grain yields were obtained 

2 

by harvesting four 1.5 m plots from the sprayed areas and four adjacent 
2 

1.5 m areas on which leaf rust was not controlled. Leaf rust severity 
estimates were made at the time of spray application using the modified 
Cobb scale (Peterson 1948). 


RESULTS 

The DCP can be used to collect and transmit data with minimal main- 
tenance. Two 12V storage batteries, used to power the system, were 
replaced every four weeks, when they needed recharging. At the same time, 
sensors were routinely checked. Occasional down-time was experienced due 
to an electrical storm or rodents chewing transmission cables. 



Data (IBM cards) were received from NASA 5 to 14 days after transmission 
to the satellite. Normally, 8 to 12 transmissions were received by the 
satellite each day from each DCP (Fig. 1). Figs. 2 and 3 show typical 
meteorological information acquired by the DCP. 

Maximum and minimum temperatures and hours of dew occurring each day 
as recorded by the two DCPs for the seven day period April 23 to 29 were 
used together with other biological data taken in the test fields to 
predict leaf rust severities that would be expected on May 6, 13, 20, 29, 
near each DCP. These predicted severities were used in the leaf rust loss 
equations to predict the percent reduction in wheat yields to be expected 
due to leaf rust development. Using meteorological data obtained from 
the DCPs a 9% decrease in wheat yields due to leaf rust was predicted for 
each site, and compared favorably with actual decreases in yield of 9% 
and 12% (Table 2) . Actual yield reductions were obtained by comparison 
of yields of the controlled (healthy) and the noncontrolled (rusted areas). 

The manpower required for obtaining data in remote areas could be 
minimized by use of DCPs* Our results indicate DCP-derived data can be 
effectively used in existing disease prediction equations. 
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Tabic 1. Data from DCPs in Riley County wheat fields transmitted at the 
two periods. 


9:00-11:00 TRANSMISSION 


21:00-23:00 TRANSMISSION 


1 

2 

3 

4 

5 

6 
7 


instantaneous relative humidity 
cumulative hours free moistures 
minimum temperature 00:00-11:30 
instantaneous temperature 
instantaneous incoming visible 
instantaneous reflected visible 
instantaneous incoming infrared 


instantaneous soil moisture 
cumulative hours free moisture 
maximum temperature 11:30-00:00 
instantaneous temperature 
14:00 incoming visible 
14:00 reflected visible 
14:00 incoming infrared 


8 


instantaneous reflected infrared 


14:00 reflected infrared 






table 2. Leaf rust severity estimates, wheat yields and percentages loss observed in commercial tields iron which ol i 
data were obtained. 


^Location 


^pril 26 


If 

%DS- 


May 18 May 25 June 4 

GS 7.DS GS 7-DS GS %DS 


June 11 
GS %DS 


kg/ha 


% Los s 

actual 


ERICHS EN 
| Healthy 


Rusted 


{farther 

Healthy 


Rusted 


M 

1 

SD 

3 

2913 

M 

1 

SD 

3 

2887 

M 

1 

SD 

3 

2911 

M 

1 

SD 

3 

2975 

2921 

M 

20 

SD 

40 

2627 

M 

20 

SD 

40 '• 

2690 

M 

20 

SD 

40 

2424 

M 

20 

SD 

40 

2511 

2562 

M 

.5 

SD 

.5 

3135 

M 

.5 

SD 

.5 

3050 

M 

.5 

SD 

.5 

3227 

M 

.5 

SD 

.5 

3143 

3139 

M 

15 

SD 

30 

2882 

M 

15 

SD 

30 

2762 

M 

15 

SD 

30 

2391 

M 

15 

SD 

30 

2854 


—^Growth stage; EJ (early joint); H (heading); A (anthesis); M (milk); SD (soft dough) 


12 / 

—disease severity 
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SEASONAL CANOPY REFLECTANCE PATTERNS OF WHEAT, 

SORGHUM AND SOYBEAN 1 
by 

E. T. Kanemasu^ 

ABSTRACT 

Reflectance characteristics of agronomic crops are of major importance 
in the energy exchanges of a surface. In addition, unique reflectance 
patterns may be an aid in crop identification by means of remote sensing. 
Our study suggests that the ratio of the reflectances of the 545-nm to 

the 655-nm wavebands provides information about the viewed surface, regard- 

1 

less of the crop. The reflectance ratio is less than unity early and late 
in the growing season. For all crops studied, the ratio closely followed 
crop growth and development and appeared to be more desirable than the 
near- infrared reflectance as an index of growth. 


Contribution No. 1385 Evapo transpiration Laboratory, Agronomy Department, 
Kansas Agricultural Experiment Station, Kansas State University, Manhattan, 
Kansas 66506. This study was partially supported by the National Aero- 
nautics and Space Administration. Received 

► 

'Assistant Professor of Microcliraatology , Evapotranspiration Laboratory, 
Department of Agronomy, Kansas State University, Manhattan, Kansas 66506. 
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INTRODUCTION 

Canopy reflectance patterns are important to the radiative balance of 
a crop and as possible discrimination features for remote sensing appli- 
cations. Individual leaf reflectance can provide valuable information 
after a canopy cover becomes complete. However, in most cases the condi- 
tion of the viewed surface is not known and discrimination analysis must 
be performed on canopy-reflectance data that may not be easily interpreted 
from leaf-reflectance data. 

This study focused on determining canopy- reflectance patterns that 
would allow the surface condition to be determined. 

METHODS AND MATERIALS 

TWo fields each of wheat (Triticum aestivum L. cv. Scout), sorghum 
( Sorghum bicolor L. Moench Pioneer var. 846) and soybeans ( Glycine max 
L. cv. Clark 63) were selected in a bottomland area where one field was a 
dark-colored, silty clay loam (lat. 39°08 , N, long. 96°37.5'W) and the other 
a light-colored, silt loam (lat. 30°08 , N, long. 96°37 , W). 

Growing conditions were considered normal (compared to previous 
seasons) and adequate soil moisture was maintained by precipitation and 
irrigation. Specific growth stages and leaf area indices (leaf area to 
ground area) were recorded periodically. All the plants within a half 
square meter area were taken to the laboratory and leaf area was deter- 
mined with an optical palnimeter. 

To determine the spectral hemispherical reflectance of the canopy, 
the sensor head of a portable spectral radiometer (LI-187, Lambda Instrument 
Co.) was pointed upward and downward. The spectroradiometer has 9 full-scale 
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ranges from 0,3 to 3000 watts m (ym ). The sensor head (8*2-cm diameter) 
consists of seven miniature sensors covering the visible and near infrared 
wavelengths (Table 1). The D, F, and G sensors (545, 655, and 750 nm) 
closely correspond to Earth Resources Technology Satellite (ERTS-1) bands 
4, 5, and 6* During the measurements the sensor head was positioned 
approximately 1*5 to 2.0 meters above the canopy. Azimuthal direction was 
kept constant relative to the sun; the observer always faced the sun. 
Measurements were taken only on clear days. 

RESULTS AND DISCUSSION 

During a significant part of a crop’s growing season, bare soil is 
exposed. As the plant grows, less soil is exposed and the soil’s reflec- 
tivity becomes less important in overall canopy reflectance. In our study 
soil reflectance was strongly influenced by the surface moisture (% by 
weight) of the silty clay loam (Fig. 1). The near-infrared wavelength 
band was the most sensitive to surface moisture. The longer the wave- 
length, the higher was the reflectance, a relationship which was consistent 
with the findings of other investigators (Bowers and Hanks, 1965). 

Fig. 2 shows the midday spectral reflectances for wheat, sorghum, 
and soybeans at growth stages early, middle, and late in the season. 

The highest reflectance was in the near infrared at midseason. In 
addition, at that time reflectance was greater at 545 nm than at 655 nm; 
the reverse was true early and late in the season. That suggests that 
the ratio of the reflectances at 545 and 655 may be an indicator of 
soil exposure early in the season and of crop maturity late in the 


season. 



The effect of solar elevation on reflectance usually hinders the 
interpreting of surface conditions from reflectance data (Suits, 1972). 
Therefore, the effect of sun angle on the reflectance ratio must be known 
before reflectance data can be interpreted correctly. Fig. 3 shows the 
ratios of 545 and of 655 nm reflectance with solar elevation for wheat, 
sorghum, soybeans, and bare soil. For each canopy, the ratio remained 
constant with increased solar elevation. For a mature crop, the ratio was 
about 1.3; for a bare soil, about 0.8. The near-infrared reflectance for 
wheat and sorghum decreased with increased solar elevation; but that for 
soybeans varied somewhat, perhaps because leaf angles changed with solar 
elevation (Fuchs et al. 1972). 

Because the reflectance ratio apparently is not influenced by solar 
elevation, reflectance ratios for wheat, soybeans, and sorghum can be com- 
pared over a large portion of the growing season without serious error due 
to changes in sun angle. The results of such measurements are shown in Fig. 4. 
Fig. 4a shows that the wheat on the light-colored soil had higher near-infrared 
reflectance than that of the dark— colored soil (early in the season) ; the 
near-infrared reflectance did not start to increase until the late-joint 
growth stage. The reflectance ratio apparently followed the leaf area 
index curve (Fig. 5a). The ratio increased above unity at a leaf area 
index of about 1.0 and remained above unity during maximum growth, then 
decreased below unity at maturity (leaf area index < 1.0). The reflectance 
ratio was greater for the field with the greater leaf area index. Similarly, 
Fig. 4b shows the same trends in the near-infrared reflectance and the 
reflectance ratio for soybeans. At 120 days after planting, the soybean 
leaves yellowed from an infection of bacterial pustule and the reflectance 
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ratio decreased to less than one on the infected soybeans. Thus, the 
reflectance ratio may serve as an indicator of physiological stress , such 
as brought about by disease, insects, drought or by normal maturation of 
the plant. 

Measurements on narrow— row (46 cm) and wide-row (92 cm) sorghum (grown 
on dark-colored soil) which were made over the entire growing season (Fig. 
4c), showed the effect of canopy cover more clearly. Plant density was 
maintained at 17 plants per square meter in each field (2 ha). The close- 
row spacing closed its canopy early in the season while the wide-row canopy 
never completely closed; both fields had similar leaf area indices. Near- 
infrared reflectance varied greatly early in the season, presumably because 
of changes in surface-moisture content. The reflectance ratio of the 
narrow- row sorghum increased to above unity at a leaf area index of about 
1.0 (Fig. 5b) which also corresponded to near 90% cover (visual estimate); 
the wide row sorghum did not reach a ratio of unity until a leaf area of 
2.5 (approximately 85% cover). The reflectance pattern for the sorghum 
(7 6- cm rows) on the light- colored field (not shown) closely followed that 
of the wide-row sorghum, illustrating that reflectance ratio may follow 
percentage cover more closely than leaf area index. The reflectance ratio 
decreased to below unity late in the season, even though the leaf area index 
was greater than one* The percentage cover at that time was about 60%. 

Table 2? shows the linear regression equations derived from Figs. 4 
and 5. They were obtained from single-field measurements for soybeans and 
sorghum (because leaf-area measurements were incomplete on the light-colored 
field); for wheat, data from three fields were pooled. The percentage 
cover was continuously estimated only for the wheat fields. Neither 
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reflectance ratio nor near- inf rared reflectance (NIR) offered a unique 
equation for relating reflectance to leaf area index (LAI) for all crops. 

The correlation coefficients were highest for soybeans and lowest for 
wheat. Where data from several wheat fields were examined the correlation 
coefficient was greater for the ratio than for the near-infrared reflectance. 

This study suggests that the reflectance ratio of the 545- to 655-nm 
wavelengths may serve as benchmarks for crop growth and possibly for 
indicating percentage cover. When the ratio is less than unity, soil 
reflectance dominates canopy reflectance. When the crop matures, the 
ratio decreases to less than one. The ratio apparently is a better indica- 
tor of crop growth than is the 750 nm reflectance. The reflectance ratio 
alone does not appear to discriminate between crop species, but should be 
a valuable parameter when used with other recognition processes. The 
two wavelengths involved in the ratio correspond to multispectral scanner 
bands 4 and 5 (MSS4 and MSS5) of ERTS-1. A study is underway to test the 
feasibility of using the ratio of MSS4 and MSS5 for wheat growth and disease 


detection. 
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Table 1. 

Optical characteristics 
radiometer 

of the spectral 

Sensor 

Center wavelength 
(nm) 

Band width 
(nm) 

A 

415 

40 

B 

450 

52.5 

C 

500 

45 

D 

545 

35 

E 

600 

40.2 

F 

655 

45 

G 

750 

80 
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Table 2, Linear regression equations and correlation coefficients for wide-row 
and narrow-row sorghum, wheat and soybeans, LAI is leaf area index; 
NIR is near-infrared reflectance. 


Linear regression equation Correlation coefficient 


Sorghum (wide-row) 
Sorghum (narrow row) 
Wheat 

Soybean 


LAI = 10.93 x (ratio)- 
LAI = 0.26 x (7.NIR) - 

LAI = 9.18 x (ratio)- 

LAX = 0.23 x (7.NIR) - 

LAI = 5.06 x (ratio)- 

LAI = 0.13 x (7.NIR) - 


8.37 

0.89 

2.70 

0.84 

6.90 

0.80 

3.24 

0.87 


4.07 0.75 

1.67 0.64 

7. cover = 109.88x(ratio)-63.71 0.87 

% cover = 2.85x(7.NIR) -19.24 0.72 

LAI = 13.67 x (ratio)- 11.28 0.96 

LAI - 0.296 x (7.NIR)- 5.10 0.98 





flEFLEC TANCE 



2a. Spectral reflectance of wheat on February 23, 1973 
(dormancy). May 16 (heading), and June 11 (hard 
do 
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g. 3. Effect of solar elevation (O-horizon) on the 

reflectance ratio of the 545 to 655 nm waveband 
and on the near infrared reflectance (750 nm) 
of wheat (boot stage), sorghum (early heading) , 
soybean (late podding), and dry silty clay loam 
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Fig, 4a * Seasonal trends in the reflectance ratio and near 
infrared reflectance of a silty clay loam wheat 
field (dark-colored) and a* silt loam wheat field 
(light-colored ) , 
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Fig* 4b. Seasonal trends in the reflectance ratio and 
near infrared reflectance; of a silty clay 
loam soybean field (dark) and a silt loam field 
(light) . 
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FLEXIBLE DCP INTERFACE 1 
by 

2 

H. Schimmelpf ennig and E. T. Kanemasu 
ABSTRACT 

A user of an ERTS data collection system (DCS) must supply the sensors 
and signal-conditioning interface. The electronic interface must be 
compatible with the NASA-furnished data collection platform (DCP). We 
describe here a "universal" signal-conditioning system for use with a 
wide range of environmental sensors. 

The interface is environmentally and electronically compatible with 
the DCP and has operated satisfactorily for a complete winter wheat growing 
season in Kansas. 


Contribution No. 1397, Agronomy Department, Evapo transpiration Laboratory, 
Kansas Agricultural Experiment Station, Kansas State University, Manhattan, 
Kansas 66506. This study was partially supported by the National Aero- 
nautics and Space Administration. Received . 

2 

Electronic Technician and Assistant Professor of Microclimatology, Evapo- 
transpiration Laboratory, Department of Agronomy, Kansas State University, 
Manhattan, Kansas 66506. 
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INTRODUCTION 

The Data Collection Platform (DCP) is a field-deployable, automatic, 
data-relay terminal that can be located in remote areas to gather infor- 
mation for specific applications or to complement imagery information 
received from the ERTS system. The DCP (which consists of an electronic 
unit, an antenna assembly, and an interconnecting cable) accepts sensor 
input data. The accepted data can be in the format of eight analog inputs; 
eight 8-bit, serial-digital inputs; eight, 8-bit, parallel-digital inputs; 
or combinations of these formats. To gather these data, the user must 

supply a power source, sensors and signal-condition system for his specific 
_ 3 

requirements . 

The data collection platform interface (DCPI), described here, is a 
"universal" signal-conditioning system that accepts inputs from 8 sensors 
of almost any type and interfaces them to the DCPs The DCPI contains 
power supplies, control logic, memory, and signal-conditioning modules 
for each sensor. An alteration in a signal-conditioning function can be 
easily affected by inserting the proper module in the DCPI. Modules can 
perform a variety of signal conditioning functions such as amplifying, 
linearizing, integrating, totaling counts, sample-hold, and comparing 
previous readings. Appendix A lists the sensor and signal-conditioning 
system used in our wheat study. 

TECHNICAL DESCRIPTION 

Figure 1 is a functional block diagram of the DCPI. Two 12-volt car 
batteries power both the DCP and DCPI. Time between recharging the battery 

3 

Earth Resources Technology Satellite Data Collection Platform Field 
Installation, Operation and Maintenance Manual. NASA. Goddard Space 
Flight Center. 



is dependent on the particular function modules; for our study, it was 
about 2 months under average conditions. 

The 24-hour clock sends time data (in 10 minute increments) to the 
modules. The main controller sequentially interrogates the signal- 
conditioning modules and controls the data storage cycle. The analog 
to digital converter converts analog to digital data for storage in the 
memory. 

The 64-bit memory stores the data from the 8 sensors (8 bits for each 
sensor). Data are transferred to the DCP for transmission, or back to 
the modules for comparison to current data. 

Light emitting diodes on the instrument panel monitor the data and 
provide information during routine field inspections. The DCP transmits 
at 3-minute intervals 24 hours a day. The TRANSMIT CLOCK from the DCP 
shifts 64 data bits (8 per sensor) from the memory to the transmitter. 

After the transmission, the DCPI scans the eight channels and stores any 
new data in the memory. The logic flow in a scan cycle is best illustrated 
by specific example, like maximum-minimum temperature. 

Example Logic Plow . At the start of a scan, the main controller turns on 
the power supply and, after a 5-second delay, interrogates channel 1. 
Following interrogation, logic control shifts from the main controller to 
the interrogated module, which then connects itself (via a set of CMOS 
switches) to the 10 control-buss lines. Assume that the sensor is a 
temperature sensor and the module is designed to find the maximum air 
temperature during the day. The maximum temperature module converts the 
resistance of the temperature sensor to a properly scaled, analog voltage, 
and then generates a START OF CONVERSION (SOC) command, causing the ANALOG- 
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TO-DIGITAL CONVERTER (ADC) to convert the analog voltage from the module 
(ADO) line to digital data* After the digital conversion, the ADC generates 
an END-OF-CONVERSION (EOC) signal informing the module that the conversion 
process is complete. 

Previous maximum temperature is stored both in an 8-bit, parallel 
output shift register on the module and also in the main memory. These 
digital data are converted (at the module) to analog and compared with 
current temperature. Two cases present themselves: (1) the present 

temperature is higher than the maximum or (2) the present temperature is 
lower than the maximum. 

Case 1 . The module sets the DATA, REPLACE- SAVE (D, R-S) line to 
REPLACE and generates a START-OF- STORE (SOST) command. The main controller 
accepts the command and generates an 8-pulse, shift-clock train (SHFT CLK) . 
The SHFT CLK shifts the digital data from the ADC to the input of the memory 
and shifts the memory 8 places, thus transferring the new data into memory. 
The DIGITAL DATA in (DDI) line transfers the new temperature reading from 
the ADC to the module using the SHFT CLK signal for synchronization. 

Case 2 . The module sets the D, R-S line to SAVE and generates an SOST 
command. The main controller again generates SHFT CLK* The replace-save 
switch at the memory input is now set at SAVE and the memory cycles the 
stored data, leaving them unchanged in memory. 

At the end of the storage cycle the main controller generates an 
END-OF-STORE (EOST) command, which resets a latch on the maximum temper- 
ature, and then interrogates Channel 2. The same type of control chain 
now occurs with module 2. This sequence continues until all 8 channels 


have been scanned. 



Digital data also can be gated from the module to the main memory. 
Module 2, hours of free moisture, would set the DATA, ANALOG-DIGITAL (D,A-D) 
line to digital. Digital data gated by SHFT CLK proceeds from the module 
on the DIGITAL-DATA-OUT (DDO) line, through the analog-digital switch, 
directly to the memory. The replace— save function is operable in the 
digital as well as the analog data mode. 

The sequence of INT, SOC, EOC, SOST, and EOST is repeated for each 
module. After the last module is interrogated, the interface turns off 
everything except the continuous power. 

Control Cards Al and A2 (Figs. 2 and 3) . 

The DCPI control section includes two circuit cards, Al and A2. Al 
contains the 64— bit memory and most of the control logic. A2 contains 
the analog-to-digital converter and display drivers. 

At the start of a DCP transmission the DATA GATE drops low and remains 
for its 80 ms warmup-transmit cycle. Ql, Q2, and Q 5 interface between the 
5-volt DCP TTL logic and the 12-volt DCPI CMOS logic. U10D turns on V x> 
a 5-volt supply for Q5 (needed during transmit only). 

If a manual scan is not in progress, the LOW at the output of U9D is 
gated through U13B to the mode control on the 64-bit memory (U12). U12 is 
now in the recirculate mode and, when clocked, its data bits will leave 
Q, through U17 to the DCP. 

At the end of the transmission U18 latches turning on a 5-volt supply 

(V c ) and the main ± 15-volt supply (V ) . After a 5-second warmup, counter 
5 b 

U26 is reset and advanced to channel 1, sending a HIGH interrogate signal 
to module 1. Module 1 returns its analog data (If any) to the analog-to- 
digital converter on control card A2. Module 1 next returns a high SOC 
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command to the ADC. The network of U3, U4, and U5 provides a delay, after 
the analog signal reaches the LH0042 buffer amplifier, before the A to D 
conversion can start. This network also assures that SOC commands from two 
consecutive modules will be u see" as two HIGHS, not one long continues 
HIGH. 

Upon completing the A to D conversion, the ADC returns an EOC pulse 

to the module. When ready, the module sets the D, A-D line to 1 if it is 

to send digital data for storage. It sets the D, S-R (DATA, SAVE-REPLACE) 
to 0 or 1 depending if the new data (either from the ADC or digital data 

from the module) are to be retained. Next, the module returns a START OF 

STORE (SOST) command. The 8-pulse, shift gene ator (U19, U15, U20) sends 
8-clock pulses to the memory (U12) ; the modules; and on card A2 to Ul, the 
ADC to serial shift register, and to U5, the memory-to display shift 
register. U24 gates the digital data from the ADC shift register or the 
module to the memory. When the storage cycle is complete, an END OF STORE 
(EOST) signal is sent back to the module and U26 is advanced to interrogate 
the next module. After the last module has been interrogated, U14B resets 
latch U18A, turning off the power supplies and ending the scanning cycle. 

The scanning cycle can be run under manual control to observe data 
and control states at each important step, thus facilitating trouble 
shooting. 

Pressing the MANUAL SCAN START button sets latch U18B, which blocks 
out interferring signals from the DCP transmitter and SOST signals from 
the modules. U18A latches on and U26 advances to interrogate module 1. 
Panel LED's now display the output of the analog to digital converter, 
whether the data in memory will be saved or replaced, and whether analog 
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or digital data from the module will be stored. Pressing the STEP button 
starts the store cycle. The LED’s now display the memory contents for channel 
1: either the new channel 1 data, or the previous channel 1 data which have 

been retained. Pressing the STEP button again advances U26 to channel 2, 
Switching the DCAN, MEMORY CYCLE switch to the MEMORY CYCLE position locks 
the memory in the recirculate mode, thus allowing a review of the memory 
contents unaltered by a scan cycle. 

TIME CLOCK (Fig. 4) . 

A crystal-controlled clock sends 24-hour time information to each 
module. This enables the modules to operate on time dependent data. 

Crystal-controlled oscillator U1 runs at 27.96 KHz. U4 and U8 divide 
this frequency down to 1 pulse per 10 minutes. U3, U6, and U10 give time 
outputs in 10-minute increments through 24 hours. The time outputs are 
bussed in parallel to all the modules. 

U7 gives a 14:00 signal to the four radiation modules to avoid adding 
identical decoding circuitry to the four modules. 

MODULES 

RADIATION MODULES (Fig. 5) . 

The radiation module is designed for silicon photocells. TJ1 is a 
current-to-voltage converter. Rl sets the gain. U2 is a buffer for the 
100-second, RC network R7, C2. 

A time signal to Pin 9 causes the 14:00 data to be retained in memory 
for the evening transmission. From 00:00 to 13:59 the DCP transmits 


instantaneous data. 
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MAXIMUM-MINIMUM- INSTANTANEOUS TEMPERATURE (Fig. 6) , 

This module is divided into two circuit cards and supplies data to 
two channels. The module outputs instantaneous temperature on one channel. 
The other channel is the maximum temperature between 11:30 and 23:59, or 
the minimum temperature between 00:00 and 11:30, Unless a front passes 
through, the true maximum and minimum temperatures are transmitted. 

U8 and U9 form a linear thermistor thermometer using a YSI (Yellow 
Spring Instruments) thermalinear network. When Pin 19 is interrogated, 
instantaneous temperature is sent to the controller, which converts the 
instantaneous analog signal to digital data. The digital data are passed 
back on 10 parallel lines to the input of 10-bit latch. 

The returning EOC signal strobes U26, latching SAVE-REPLACE LATCH 
U6A and U6B if a new maximum or minimum is present. This enters the new 
maximum or minimum in U3. The DAC converts this to analog data. When 
the next channel is scanned, U4B passes the maximum or minimum to the 
controller. 

RELATIVE HUMIDITY AND SOIL MOISTURE (Pig. 7) . 

This module transmits relative humidity from 00:00 to 13:59 and soil 
moisture from 14:00 to 23:59. 

Ul is a 1200-Hz oscillator. U1C gives a high pulse during the last 
1/4 of the 300-Hz square wave at the output of U2B. U3's output is a ± 1 
volt square wave which drives the EH and SM sensors through DC blocking 
tantalum capacitor pairs. U4 and U5 convert current through the sensor 
to voltage. 

The output of U4 and U5 is a square wave with a high spike on the 
front. The spike width is proportional to the lead wire capacitance and 


is an error term. 
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U6 and U7 are precision full-wave rectifiers. U8 C and D and U9 form 
a sample-hold circuit that eliminates all but the last 1/4 of the wave form. 
This technique gets rid of the lead-wire capacitance error. U8A and -B 
gate the RH or SM data to the next stages, depending on time of day. 

LOG AMP is a logarithmic amplifier, which straightens the RH and SM 
curves. U10 and Ull add offset and gain to put the signal in the final 
form. 

The relative humidity sensor works over a 2,000 ohm to 2 megohm range. 
Soil moisture is read over a 200-ohm to 20,000-ohm range. 

HOURS FREE MOISTURE (Fig. 8) . 

The presence of dew or rain is sensed by the lowered resistance of 
a bifilar grid. U1 drives, a 27 Hz square-wave through the sensor and Rl, 

R2. Cl, C2, and R4 rectify and filter the output. Operational amplifier 
U2 compares this signal with a fraction of the logic supply. 

If the sensor is wet, a 36.621 ms/cycle signal from the time clock 
is gated into 20-stage ripple counter U3 and 114. The last 8 bits of the 
counter total 256 bits at 5 minutes per bit. Hours of free moisture are 
obtained by subtracting successive transmissions. 

POWER SUPPLIES (Fig. 9) . 

Two 12-volt storage batteries are the main power source. The DCP 
requires 24 volts and the DCPI 12 volts (Fig. 9). 

V, 

Continuous 12-volt power runs all the CMOS logic in the DCPI. A 12 
volt to ± 15 volt converter supplies continuous operational amplifier 
power if needed (Fig. 9). During a transmission , a 5-volt supply turns 
on to interface the CMOS to TTL. During a scan, a high power 12-volt to 
± 15-volt converter supplies power to the operational amplifiers, while 
a high power 5- volt supply powers the ADC and some TTL logic in the DCPI. 
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CH 1 


CH 2 


CH 3 


APPENDIX Cl 

DCP SENSOR AND SIGNAL CONDITIONING CHARACTERISTICS 

Relative humidity and soil moisture 
Position: in the canopy 

Sensor: Type, Relative Humidity FCRC-11 sulfonated polystyrene 

(Phys- Chemical Research Corp.) 

Span: 0-100% RH ± RH 

Accuracy: ± 1% RH 

Position: 30-cm depth 

Sensor: Type, soil moisture block CEL-WFD 

(Beckman Instruments) 

Accuracy: 2% of reading 

Signal Conditioning: Type, Lead-wire, capacitance-eliminating, 

AC ohmmeter with log amplifier for lineari- 
zation. 

Hours of free moisture (dew and rain) 

Position: above canopy 

Sensor: Type, Bifilar array on printed-circuit board (G-1G epoxy 

base) 

Signal Conditioning: Type, level-detecting, AC ohmmeter 

Accuracy, ± 20 MS per change of sensor state 
with no additional cumulative error. 

Maximum and minimum temperature 
Position: at top of canopy 

Sensor: Type, YSI series 700 thermalinear thermistor probe. 

(Yellow Springs Instruments) 
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CH 4. 


CH 5. 


CH 6 


CH 7 


Signal Conditioning: Type, digital storage of max (min) and analog 

comparison with present temperature. 

Accuracy: ± 0.25°C 

Instantaneous temperature 
Position: at top of canopy 

Sensor type and accuracy: same as channel 3. 

Signal conditioning: Type, thermalinear thermistor bridge 

Accuracy: ± 0.15°C. 

Incoming visible radiation 

Position: approximate 2m above soil surface 

Sensor: Type, silicon photocell SBC 255 

Instrument: (1) cosine corrected head 

(2) 6 mm heat adsorbing glass (KG-3) 

(3) diffusing plastic 

(4) wratten 26 filter 

Response: 590 to 720 nm 

Construction: Built by E. T. Laboratory 

Signal Conditioning: Type, Signal averaging filter 

Accuracy: ± 0.3% 

Response time: 10 to 90% - 220 seconds 

Reflected visible 

Same as Ch 1 except sensor is faced downward 

Incoming near-infrared radiation 
Position: 1.5 m above soil surface 

Sensor: Type, Silicon photocell 
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Instrument: (1) cosine corrected head 

(2) diffusing plastic 

(3) wratten 88A 

Construction: Built by E. T. Laboratory 

Response: 730 to 1000 nm 

CH 8. Reflected near infrared 

Same as Cb 3 except sensor faced downward. 
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APPENDIX C2 

PARTS LIST POR CONTROL BOARDS A1 AND A2 

U2 - Noninverting buffer, RCA CD4050AE 

(schematic shows this as an inverting buffer) 

U12 - 64-bit shift register, RCA CD4031AE 

U17, U18, U20 - dual flip-flop, RCA CD4013AE 

U19 - counter, RCA CD4022AE 

U26 - counter, RCA CD4017AE 

Rl, R4 - 180K and R2, R40 - 120K 

R3, Rll, R12, R14, $17, $19, $27, $28, $35, $37, $38 - 100K 

R5, R6 - 390K 

R7, R8, R9 - 82K 

RIO - 39K 

R13 - 12K 

R16, R18, R26, R29, (134, R39 - 1.2K 
R21 - 470K 
R23 - IK 
R24 - 68K 

R25, R30, R31 - 270K 
R32 - 560K 
R33 - 220K 
R41, R42 - 1.5M 

Cl, C6, C7, C8, CIO, Cll, C12 - .OlyF 
C2, C3 - lOyF 
C4, C5 - 470 pF 
C9 - .002 yF 



Notes : 


015 


(1) All digital logic parts are RCA COSMOS except where noted 
differently on schematics 

(2) COMOS CD4016AE transmission gates are shown as 



(3) NPN transistors are 2N222A 
PNP transistors are 2N29Q7A 
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Fig. 5. Radiation Modules 
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INTRODUCTION 


An expanding population has brought about an awareness that there are 
only limited resources on the Earth. This realization comes at a time 
when resource use is greater than ever before. Adequate informational 
techniques are necessary for improved resource development. These tech- 
niques can aid in wise resource management. 

The magnitude of the data required for improved resource management 
has led to the development of automatic recognition techniques for agricul- 
ture. These systems utilize remote sensing from aircraft and spacecraft. 

Earth Resources Technology Satellite program is a major step in combining 
space and remote sensing technologies into a system for developing and 
demonstrating the techniques for efficient management of the Earth’s re- 
sources (NASA Earth Resources Technology Satellite Data Users Handbook, 1972). 

Over 400 million acres of land are irrigated in the world (Israelsen 
and Hansen, 1967). Some of the water applied is needlessly lost by excess 
applications. Irrigation scheduling can help to better conserve this valu- 
able resource. One method of scheduling irrigation requires the determination 
of crop water use (evapotranspiration) . Actual evapotranspiration is depen- 
dent upon potential evapotranspiration and a crop coefficient. One possible 
approach to predicting the crop coefficient is the use of a plant’s actual 
growth which may be determined by its reflection of solar radiation from 
the plant canopy (Myers et al . , 1966). If this method is to be used, the 
relationship between reflectance, soil moisture and vegetative growth must 
be established. 

The purpose of this research is to evaluate reflectance for prediction 
of soil moisture and vegetative growth, and to determine the. feasibility 
of using vegetative growth to evaluate the winter wheat crop coefficient. 
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REVIEW OF LITERATURE 
Remote Sensing 

Remote sensing refers to the acquiring of data at a distance by detecting 
the radiant energy which the object either reflects or emits . Detection 
devices can be field spectrometers and cameras or instruments designed for 
installation in aircraft and space vehicles. 

Albedo is the ratio of the entire solar radiation spectrum reflected 
from a body to the total incident radiation (Ashburn and Weldon, 1956), 
while reflectance is the ratio of reflected radiation to the total incident 
radiation at a specific wavelength. At any specified wavelength. Reflectance 
+ Absorptance + Transmittance = 1. Transmittance of any opaque material is 
zero; thus a decrease in reflectance will cause an equal increase in absorp- 
tion. 

Physical Properties that Affect Reflectance 


Soil Factors 

The albedo of various soil surfaces was compiled by Kondrat’yev (1965). 
The soils had extremely variable albedos. The variability was attributed 
to the different soil color, soil moisture content, organic matter and 
particle size. The soil moisture content was considered the most important 
factor. He pointed out that a decrease in albedo with an increase in mois- 
ture was due to water's low albedo. Bowers (1971) indicated that the 
relationship between soil moisture and reflectance is precise enough to 
utilize reflectance techniques to measure surface moisture (Fig. 1). However, 
due to the soil color, a calibration is necessary for each soil type. 



reflectance 



Fig. 1. Reflectance from Newtonia Silty Clay Loam at 
Different Soil Moisture Percentages (Figure 
reproduced from Bowers, 1971)* 
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Allen and Sewell (1973) concluded that the use of infrared films and 

electronic scanner detectors could detect fallow soil moisture over a range 

of 1 to 24 percent dry weight* Their prediction equations for both the 

surface soil moisture and soil moisture at the 4 inch depth had regression 
2 

coefficients (R ) of at least 0.94. 

Organic matter also influences reflectance. A study by Bowers (1971) 
shows that an oxidized soil sample compared to the check or control sample 
has a greater reflectance. He also states that some of the change could 
have been due to oxidation of the carbonates , although in one soil no 
carbonate was detected. 

Bowers (1971) and Myers and Allen (1968) also reported that particle 
size has an effect on reflectance. In most cases an increase in particle 
size decreased the reflectance. This was due to the fine particles filling 
the volume more completely, thus a more even surface. Coarse aggregates, 
having an irregular shape, formed a large number of pores and cracks in the 
surface* When the soil surface was wet and pulverized there was very little 
difference in reflectance from soils, instead the real contrast was at a 
low moisture content. 

Vegetative Factors 

The main factor that causes variation in reflectance from crop canopies 
is leaf density or leaf area index. Leaf area index is defined as the ratio 
of the leaf area to soil area. Stanhill ej: al. (1968) reported that leaf 
area index is linearly correlated to albedo or shortwave reflection. The 
plant albedo increases with increasing plant development to a maximum at 
full plant canopy. The suggested model indicates internal trapping of 
radiation, which decreases albedo. Internal trapping is almost complete 
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after the second reflection with hardly any effect by height after a minimum 
value. In the near infrared region, reflectance increased 17 percent with 
two leaf layers and only slightly more for each additional leaf layer. When 
the crop cover is incomplete all of the soil factors mentioned previously, 
including soil color, soil moisture, particle size and organic matter, caused 
variation in reflectance. In addition, leaf reflectance also is affected 
by stand geometry and leaf morphology, most significantly in the near infrared 
region (Gates , 1965) , as well as the variety and relative maturity of the 
crop (Remote Multispectral Sensing in Agriculture, 1970). 

A comparison of different varieties of a crop by Interpretation of 
Remote Multispectral Imagery of Agricultural Crops (1967) and Remote Multi- 
spectral Sensing in Agriculture (1967) indicated that the spectral responses 
were statistically different. These differences could also have been attributed 
to variations in crop canopy or leaf area index and crop maturity. In 
mid-season it could have been due to weed infestations, diseases or farming 
practices . 

Variations of reflectance were found with spectral bands. In the visible 
region, the striking feature of the leaf spectrum was the high absorptance 
from 0.4 to 0.5 y, the reduced absorptance from 0.5 to 0.6 y, the high 
absorptance from 0.6 to 0.7 y and the low transmittance in the entire region 
(Fig. 2). This was mainly due to the chlorophyll and carotene absorption 
that predominates in this region (Remote Sensing, 1970). Sinclair, et al . 

(1973) reported that cell walls scatter the light diffusively, but the 
chlorophyll or other pigments are present to absorb the light. The 
absorbing process is a dominate factor in influencing the spectral res- 
ponse in the'visible region. If water deficits occur, the metabolic 




Fig. 2. 


Characteristic Spectral Reflectance Curve of a Green Leaf 
from Remote Multispectral Sensing in Agriculture, 1970). 


(Figure reproduced 
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processes slow down resulting in the breakdown of carbohydrates and protein 
within the plant cell. As the stress becomes more severe, accelerated 
migration of soluble leaf phosphorous and nitrogen compounds to the stem 
occurs. The loss of chlorophyll accompanying the breakdown and migration 
results in higher reflectance (David, 1969). Therefore, reflectance is 
related to the amount of plant pigments. Other factors may result in the 
loss of chlorophyll such as leaf -maturity, salinity, disease or mineral 
deficiencies. Severe nitrogen deficiences increase reflection (Remote 
Sensing, 1970), but differences in available nitrogen produce differences 
in vegetative growth (Bhangoo, 1956, Bolaria, 1956, and Monteith, 1959). 

In the near infrared region (0.7 to 1.3 p) reflectance is caused by the 
lack of pigment absorption and by the lack of absorption by liquid water 
(Remote Sensing, 1970). Sinclair, et al. (1973) suggested that reflectance 
had to occur at interfaces within the leaf where total or critical reflec- 
tance was possible. The requirements for total or critical reflectance are 
that the radiation pass from a material with a high index of refraction to 
a material with a low index of refraction and that the angle of incidence 
must be sufficiently large. The increase in reflectance as the leaves become 
more nitrogen deficient suggests that the leaves are thicker since reflec- 
tance increases exponentially as leaf thickness increases. Moisture stress 
causes physiological changes in the leaf that cause the infrared reflectance 
to decrease with an increase in moisture stress. The low absorption or 
high reflectance in this region is a distinctive feature of vegetative. 

Remote Sensing (1970) reports that of the total incident radiation which 
strikes a leaf, about 50 percent is reflected, 45 percent is transmitted 
and the remaining is absorbed. Sinclair et al. (1973) provide a more 
detailed explanation of the reflectance of an individual leaf in both the 
visible and near infrared regions. 
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Sun angle and attenuation are two factors that affect reflection, from 
an object# At low sun angles the reflectance of an object increases com- 
pared to a large sun angle# Attenuation is defined by Remote Sensing (1970) 
as including losses from a beam of radiation by either atmospheric absorption 
or scattering# In the visible region absorption plays only a minor role 
compared to scattering. Scattering is caused by interaction between radi- 
ation and small particles (dust or water droplets usually in the form of 
a cloud or haze)# 

Estimating Soil Moisture 

A large amount of time and effort has been expended in the research of 
transpiration and evaporation with only recent applications in the modeling 
of evapotranspiration for management of irrigated land. This comes at a 
time when studies indicate that the timing of irrigations and the amount of 
water applied have changed very little (Jensen et al. y 1971). If a model 
is to be used on a practical basis for irrigation scheduling, necessary 
information must be relatively simple to obtain. 

Jensen et al, (1971) have developed a computerized model to estimate 
soil moisture depletion# One of the model f s primary objectives is the 
orientation for the user instead of the researcher. To calculate the poten- 
tial evaporative flux, the Penman combination equation is used (Penman, 1963)* 
The meteorological data necessary to evaluate the equation include minimum 
and maximum daily air temperatures, daily solar radiation, dew point temper- 
ature at 8 AM and daily wind run. 

The crop coefficient used in the computer model represents the effects 
of the resistance of water movement from the soil to the evaporating surfaces, 
the resistance to the diffusion of water vapor from the surfaces to the 
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atmosphere and the amount of available energy compared to the reference crop 
(Jensen, 1968). Thus the crop coefficient is limited by the available soil 
moisture as well as the daily meteorological conditions and stage of plant 
growth. For each separate crop a coefficient must be developed for the 
model. A more detailed explanation can be obtained from Jensen at al . 

(1971) . 

Ritchie and Burnett ( 1971 ) an d Ritchie (1972) determined a nonlinear 
relationship between the leaf area index of a crop and the ratio of the 
plant’s evapotranspiration to the potential evapotranspiration. They 
reported that while an adequate supply of water is available in the soil, 
plant factors influence evapotranspiration rates. 
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INVESTIGATION 

Objectives 

This work was concerned with problems dealing with utilizing remote 
sensing data* The objectives of the study were: (1) to evaluate reflec- 

tance for prediction of soil moisture and vegetative growth, (2) to deter- 
mine the feasibility of using vegetative growth to evaluate the winter 
Wheat crop coefficient, and (3) to evaluate the winter wheat crop coeffi- 
cient in the mathematical model by Jensen <rt al. (1971) for irrigation 
scheduling* 

Equipment 

ERTS-1 satellite revolves in a circular orbit around the Earth every 
103 minutes at 914 km above sea level* The satellite travels over the 
research area in luidmoming in a north to south direction. It passes over 
any location on the Earth’s surface once every 18 days at the same time of 
day. 

The Multispectral Scanner (MSS) is a line-scanning device that operates 
in two bands of the visible spectrum and two in the near infrared. Band 4 
included the spectrum between 0.5 and 0.6 y, band 5 between 0.6 and 0.7 y, 
band 6 between 0.7 and 0.8 y and band 7 between 0.8 and 1.1 y. Fig. 3 
shows the 4 bands with the energy emitted in the solar and thermal spectrum. 
An oscillating mirror in the MSS causes light energy from a 185 km swath to 
be swept across the focus of a small telescope. At the focus is a four-by- 
six array of 24 optical fibers (6 for each band). The fibers carry the 
energy from the light through spectral filters to detectors that convert it 
to an electrical signal. An area of 79 meters square is contained in each 
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fiber. The MSS image covers 185 km square with 4 images per area. The 
imagery is relayed to ground stations and then is processed into photographs 
at Goddard Space Flight Center in Greenbelt, Maryland. The resolution 
capability reveals surface features at a scale of 1:250,000 and information 
at a scale of 1:30000. Further details of the equipment aboard the ERTS- 
1 satellite are given by NASA Earth Resources Technology Satellite Data 
Users Handbook (1972). 


Methods of Procedure 

The research was conducted on winter wheat fields approximately 30 
kilometers northwest of Garden City, Kansas. Two soil moisture treatments, 
one dryland wheat field (A) located 38“ 9.6* North latitude and 101 c 5.9' 

West longitude and one irrigated field (B) 38° 8.5' North latitude and 101° 
4.9’ West longitude, were used with approximately 60 hectares in each. 

Field B was irrigated by a center pivot sprinkler system. The two fields 
were located within 3 km of each other. The area's normal annual precipi- 
tation is 43.6 cm with about 70 percent of the precipitation during September 
through June. 

The two fields were located on Ulyssess-Richfield silt loam with an 
average organic matter of 1.5 percent and soil pH of 6.9. The exchangeable 
potassium was in excess of 560 kg per hectare. Available phosphorus in 
field A was 117 kg per hectare and in field B was 64 kg per hectare. Particle 
size analyses revealed that both field's soils contained an average of 50 
percent silt and 20 percent clay. 

Field A had been in fallow the previous year. Scout wheat was planted 
at a seeding- rate of 29 kg per hectare on September 15, 1972. The grain drill 
used had a 25.4 cm spacing between rows. By May 24, 1973, the wheat was 
completely headed and was harvested on July 5. 
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Since field B had been in wheat the previous season, the field was 
preirrigated. Anhydrous ammonia at a rate of 90 kg of nitrogen per hectare 
was applied to the field. On September 22, 1972, Eagle wheat was seeded 
at a rate of 50 kg per hectare with a row spacing of 30.48 cm. According 
to Variety Tests with Fall-Planted Small Grains (1971), Eagle wheat was a 
selection of Scout with nearly identical vegetative characteristics. Water 
was applied with the center pivot irrigation system on May 23 (3.05 cm) and 
June 2 (3.05 cm). Harvest of the wheat was completed on July 5. 

Data Collection 

Both fields A and B were divided into four square equally sized plots 
with a sampling area in the center of each plot. An additional sampling 
area was also set up in two of the plots in field A where the corners had 
been double drilled. This gave a total of six sampling areas in field A 
and four in field B. By the use of random sampling techniques, the areas 
were broken down into one meter squares, where the leaf area index and soil 
moisture were measured. 

The soil samples were gathered at the surface and at intervals of 0 
to 15, 15 to 30, 30 to 60, 60 to 91, 91 to 121, 121 to 152 and 152 to 182 cm 
with a soil sampling tube. The samples were later dried in an oven at 105 C 
until they reached a constant weight. Then the soil moistures were calculated. 

The leaf area was determined by measuring the length and breadth of 
each leaf from randomly selected plants in the one square meter and using 
the following equation (Teare and Peterson, 1971): 

LA = -0.64 + 0.813 X (1) 

where : 

2 

LA = leaf area (cm ) 

2 

X n product of length times breadth of leaf (cm ). 
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The leaf area Index is the total leaf area divided by the land surface area. 
Both soil moisture and leaf area index data were obtained within one day of 
the flights over. 

The meteorological data were from the Garden City Experiment Station. 

These data included maximum and minimum temperatures, dew point temperatures 
and wind run. Also the field capacity, permanent wilting point and bulk 
density for Ulyssess-Richfleld silt loam were obtained from the experiment 
station. This information was determined by laboratory measurements and may 
not describe the test fields accurately. Solar radiation was obtained from 
the Dodge City Weather Service while rainfall readings were taken near the 

research area. 

Data Analysis 

Using a negative transparency from ERTS-1, the general area of fields 
(A and B) was located. Then the specific fields were found by the use of 
computer printed gray scales. From the gray scales the coordinates were 
located and the numerical values were stripped off the magnetic tapes. To 
prevent any overlapping outside of the research area, one row of data points 
around the edge of the fields was eliminated. The mean and standard deviation 
of the remaining data of the four bands were calculated (Tables 1 and 2). 

Also the mean and standard deviation of point by point ratios were determined 
(Tables 1 and 2). Stepwise Deletion Multiple Regression (1973) was used 
to evaluate the relationship between reflectance, soil moisture and leaf 
area index. 

The meteorological data, as well as the soil moistures on March 22, 
were used in the computer model of evapotranspiration (Appendix, Table 11) 
developed by Jensen et al. 0-971). The original wheat crop coefficient 



Table 1. ERTS-1 Data for Field A 


Date 


MSS4 

MSS5 

MSS 6 

MSS 7 

MSS4/5 

MSS4/7 

MSS5/7 

9/22/72 

Mean 

34.75 

37.89 

38.64 

19.55 

0.918 

1.779 

1.939 


S.D.* 

1.41 

1.90 

2.18 

0.86 

0.040 

0.068 

0.080 

3/22/73 

Mean 

33.26 

32.29 

45.87 

25.25 

1.031 

1.318 

1.280 


S.D.* 

1.28 

1.58 

1.74 

0.69 

0.040 

0.055 

0.069 

5/14/73 

Mean 

29.74 

24.50 

48.11 

28.08 

1.218 

1.064 

0.877 


S.D.* 

1.69 

2.12 

1.79 

1.66 

0.066 

0.101 

0.104 

6/1/73 

Mean 

33.43 

29.48 

52.32 

29.87 

1.138 

1.121 

0.990 


S.D.* 

1.72 

2.42 

1.84 

1.04 

0.062 

0.083 

0.104 

6/19/73 

Mean 

41.14 

49.33 

55.26 

28.70 

0.835 

1.436 

1.722 


S.D.* 

1.62 

2.07 

1.49 

0.92 

0.033 

0.074 

0.090 

7/7/73 

Mean 

59.46 

78.53 

77.68 

36.36 

0.758 

1.636 

2.161 


S.D.* 

2.14 

4.25 

2.72 

1.49 

0.030 

0.061 

0.115 


*Standard deviation. 
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Table 2. ERTS-1 Data for Field B 


Date 


MSS 4 

MSS5 

MSS 6 

MSS7 

MSS4/5 

MSS4/7 

MSS5/7 

9/22/72 

Mean 

37.05 

40.41 

40.96 

20.78 

0.919 

1.786 

1.947 


S.D.* 

1.62 

2.54 

2.37 

1.02 

0.038 

0.094 

0.128 

3/22/73 

Mean 

33.54 

32.99 

41.47 

22.57 

1.019 

1.488 

1.463 


S.D.* 

1.09 

1.91 

2,15 

0.96 

0.049 

0,073 

0.088 

5/14/73 

Mean 

27.63 

19.22 

56.66 

36.78 

1.454 

0.760 

0.532 


S.D.* 

1.60 

2.68 

3.56 

3.18 

0.132 

0.109 

0.129 

6/1/73 

Mean 

26.93 

20.03 

48.66 

31.61 

1.355 

0.858 

0.638 


S.D.* 

1.32 

2.23 

3.43 

2.54 

0.111 

0.083 

0.094 

6/19/73 

Mean 

36.68 

37.94 

52.00 

29.97 

0.971 

1.227 

1.270 


S.D.* 

1.21 

3.11 

2.05 

1.56 

0.060 

0.079 

0.131 

7/7/73 

Mean 

54.46 

73.87 

77.48 

38.24 

0.739 

1.425 

1.932 


S.D.* 

2.30 

4.37 

3.39 

1.31 

0.033 

0.060 

0.100 


^Standard deviation 
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curves were evaluated first. Then curves developed by regression analysis 
from the leaf area index data were used as the crop coefficient curves. From 
the computer model, soil moisture depletions were predicted. 
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RESULTS 

Prediction of Vegetative Growth 

ERTS-1 passes over any location on the Earth's surface once every 18 
days at the same .time of day, but some dates had high percentages of cloud 
cover. Neither aerial nor ground data were collected on those days (Table 3). 
These data (Table 4) were used as a means for determining vegetative growth 
with Stepwise Deletion Multiple Regression (1973)* The July 7 data were not 
used because of the alteration of the natural vegetative growth by harvest- 
ing the wheat. The wheat threshed straw provided a stubble mulch com- 
pared to the uncut wheat. The equations that best describe vegetative 
growth were: 


LAI » 

2 . 92MSS4/5 - 2-63 , 

R 2 

= 0.95 

(2) 

LAI = 

-0.065MSS5 + 2.66 , 

R 2 

1! 

O 

• 

00 

o> 

(3) 

LAI = 

-1. 22MSS5/7 + 2-08 , 

R 2 

= 0.85 

(4) 


where 

LAI = Leaf area index 
MSS4/5 - Ratio of band 4 to band 5 
MSS 5 = Band 5 

MSS5/7 = Ratio of band 5 to band 7 
2 

R = Regression coefficient. 

For the predicted values of leaf area index to have meaning, it is nec- 
essary that a minimum or maximum* value of MSS4/5, MSS5 and MSS5/7 be set 
so that the predicted leaf area index is never negative. 

The general trend from equation 2 indicates that as the ratio of 
band 4 to band 5 increases the leaf area index increases linearly. This 



Table 3. Weather Conditions at Flight Time Over Test 
Fields . 


Date 


Weather 

Condition 

Data 

Acquired* 

September 

4, 1972 

Cloudy 

' 

September 22 , 1972 

Clear 

X 

October 

10, 1972 

Partly Cloudy 


October 

28, 1972 

Cloudy 


November 

15, 1972 

Cloudy 


December 

3, 1972 

Partly Cloudy 


December 

21, 1972 

Partly Cloudy 


January 

8, 1973 

Cloudy 


J anuary 

26, 1973 

Cloudy 


February 

13, 1973 

Rain 


March 

3, 1973 

Foggy 


March 

21, 1973 

Clear 

X 

April 

8, 1973 

Heavy Snow 


April 

26, 1973 

Rain 


May 

14, 1973 

Clear 

X 

June 

1, 1973 

Clear 

X 

June 

19, 1973 

Clear 

X 

July 

7, 1973 

Clear 

X 


*Indicates both ERTS-1 and field data taken 



Table 4. Leaf Area Index Data for Fields A and B 


Date 

Field A 

Field B 

Mean 

Standard 

Deviation 

Mean 

Standard 

Deviation 

9/22/72 

0.00 

0.00 

0.00 

0.00 

12/21/72 

0.33 

0.00 

0.12 

0.07 

3/22/73 

0.37 

0.10 

0.44 

0.07 

5/14/73 

0.97 

0.26 

1.53 

0.39 

6/1/73 

0.89 

0.25 

1.23 

0.36 

6/18/73 

0.00 

0.00 

0.00 

0.00 

7/7/73 

0.00 

0.00 

0.00 

0.00 
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means that reflectance due to plant growth in band 4 Increases faster than 
band 5 since the vegetation reflects less radiation in band 5. Equation 2 
(Fig. 4) best describes leaf area index because of its high regression 
coefficient. The ratio appears to have cancelled any soil moisture variations 
Equation 3 shows a linear relationship between leaf area index and 
'band 5. From the equation it appears soil moisture is not significant in 
band 5. Of the three equations presented, an error in band data would have 
the least effect on leaf area index as represented by the low coefficient 
of the band in equation 3. Equation 4 uses the ratio of band 5 and 
band 7 to evaluate leaf area index with no significant variation from soil 
moisture. The reflectance due to vegetation of band 7 increases at a much 
faster rate than band 5 as plant growth continues, causing a decrease in the 
ratio. 


Prediction of Soil Moisture 

The Stepwise Deletion Multiple Regression (1973) was used to help 
interpret the aerial and ground truth data available (Tables 5 and 6). 
The information for field B on March 22 was eliminated since rain fell 
before the soil moisture could be measured. Again the July 7 data were 
not used due to the stubble mulch caused by harvesting the wheat crop. 
The equations determined were: 


SM2 = 164.44 - 

■ 4.00MSS4 - 

■ 24.08LAI 

, R 2 « 0.93 

(.5) 

SM2 = 80.70 - 

1.41MSS6 + 

10.00LAI 

, R 2 = 0.80 

(6) 

SM2 = 77.92 - 

2.56MSS7 + 

20.36LAI 

, R 2 = 0.79 

(7) 




Table 5, Soil Moisture Percentages^ for Field A. 


Date 




Soil Moisture at 

Increments (cm) 



Surface 

0-15 

15-30 

30-61 

61-91 

91-122 

122-152 

152-183 

9/22/72 

Kean 

10.43 

22.97 

23.70 

21.35 

17.65 

14.35 

12.85 

13.55 


S.D. * 

2.18 

0.99 

2.58 

0.93 

3.10 

1.87 

0.75 

0.97 

12/21/72 

Mean 

34.30 

30.40 

27.70 

26.50 

24.30 

21.10 

15.70 

13.80 


S.D.* 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3/22/73 

Mean 

8.13 

22.98 

25.42 

23.82 

20.70 

16.63 

14.62 

14.98 


S.D.* 

2.59 

1.27 

1.81 

1.26 

2.02 

2.70 

2.61 

2.61 

5/14/73 

Mean 

3.82 

16.92 

17.20 

19.12 

20.35 

20.03 

18.52 

16.45 


S.D.* 

0.75 

1.87 

1.61 

1.15 

1.61 

1.52 

1.99 

2.40 

6/1/73 

Mean 

2.87 

11.12 

13.45 

15.22 

15.72 

15.95 

16.38 

15.97 


S.D.* 

0.84 

1.50 

0.63 

1.22 

2.04 

2.19 

2.00 

2,28 

6/19/73 

Mean 

0.85 

6.35 

9.80 

11.12 

10.13 

11.08 

12.47 

13.32 


S.D.* 

0.44 

0.73 

1.55 

1.31 

1.68 

1.41 

2.06 

2.01 

7/7/73 

Mean 

1.77 

15.65 

11.45 

12.92 

14.15 

15.20 

16.33 

17.33 


S.D.* 

0.21 

1.73 

0.89 

0.39 

1.00 

2.30 

2.65 

2.51 


Soil moisture percentages on dry weight basis 


*Standard deviation 



f 

Table 6. Soli Moisture Percentages for Field B. 


Soil Moisture at Increments (cm) 


Date 


Surface 

0-15 

15-30 

30-61 

61-91 

91-122 

122-152 

152-183 

9/22/72 

Mean. 

8.35 

22.40 

20.75 

16.93 

13.43 

13.00 

14.70 

16.05 


S.D.* 

1.81 

1.30 

0.34 

2.99 

2.81 

4.13 

3.62 

2.40 

12/21/72 

Mean 

16.28 

30.10 

26.27 

24.90 

19.25 

14.83 

15.70 

16.70 


S.D.* 

3.71 

5.60 

2.11 

3.19 

2.83 

3.41 

2.75 

2.88 

3/22/73 

Mean 

19.00 

27.05 

23.97 

24,15 

20.02 

14.78 

15.00 

15.97 


S.D.* 

6.44 

4.11 

1.68 

4.16 

3.49 

4.56 

3.43 

3.18 

5/14/73 

Mean 

5.58 

21.20 

16.33 

18.15 

18.38 

17.73 

16.70 

16.90 


S.D.* 

0.93 

3.27 

2.06 

3.53 

3.85 

3.75 

3.85 

2.58 

6/1/73 

Mean 

25.10 

25.47 

18.93 

15.48 

13.88 

13.68 

15.55 

16.62 


S.D.* 

12.43 

4.10 

5.16 

3.70 

2.85 

2.40 

4.12 

2.81 

6/19/73 

Mean 

2.28 

9.63 

9.23 

11.98 

10.78 

11.18 

11.48 

13.50 


S.D.* 

1.13 

2.19 

1.73 

3.88 

3.02 

1.72 

1.68 

1.91 

7/7/73 

Mean 

2.60 

17.20 

9.40 

11.43 

11.10 

10.50 

10.38 

13.15 


S.D.* 

1.25 

1.81 

2.23 

1.53 

0.67 

1.39 

1.27 

1.64 


4 * 

Soil moisture percentages on dry weight basis. 

*Standard deviation. 
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SM2 = Soil moisture dry weight at 0 to 15 cm (%) 

LAI = Leaf area index 
MSS 4 = Band 4 
MSS 6 - Band 6 
MSS 7 - Band 7 

MSS4/5 = Ratio of band 4 to band 5 
2 

R ** Regression coefficient. 

The soil moisture equation 5 indicates that an increase in leaf area 
index, with soil moisture remaining constant, decreases the reflectance in 
band 4. This could be caused by the reflectance of the soil being greater 
than the plant reflectance. Thus as the leaf area increased, more surface 
was covered by the plant canopy causing a decrease in reflectance monitored. 
The fact that soil moisture increases absorption is reaffirmed by equations 
5, 6 and 7. Equation 5 is the best equation due to its high regression 
coefficient. 

Equations 6 and 7 indicate that the reflectance of the plant is greater 
than the reflectance of the soil. An error in band reading or leaf area 
index would cause the least change in soil moisture in equation 6 due to 
the small coefficients. 

Upon substituting equation 2 into equation 5, soil moisture at 0 to 15 
cm depth became: 

SM2 - 101.11 - 4.00MSS4 - 70.31MSS4/5 (8) 

Table 7 and Fig. 5 show a comparison of soil moisture predicted by equation 
8 with the measured soil moisture. Equation 8 was developed for soil factors 
pretaining to the fields. Different soil factors would require a new equation 
to be developed for soil moisture. These factors include soil type, organic 
matter, particle size and cultural practices. 



Table 7. Predicted Soil Moisture Percentages at 0 to 15 cm from ERTS-1 Data 





Field A 



Field B 


Date 

MSS4 

Predicted 
MSS 4/5 SM2 

Actual 

SM2 

MSS 4 

Predicted 3 
MSS4/5 SM2 

Actual 

SM2 


9/22/72 

34.75 

0.918 

24.24 

22.97 

37.05 

0.919 

14.80 

22.40 

3/22/73 

33.26 

1.031 

22.25 

22.98 

33.54 

1.019 

21.85 

27.05 b 

5/14/73 

29.74 

1.218 

23.09 

16.92 

27.63 

1.454 

14.91 

21.20 

6/1/73 

33.43 

1.138 

13.86 

11.12 

26.93 

1.355 

24.69 

25.47 

6/19/73 

41.14 

0.835 

c 

6.35 

36.68 

0.971 

12.66 

9.63 


Calculated by SM2 = 101.11 - 4.00MSS4 - 70.31MSS4/5. 

^Precipitation fell after ERTS-1 flight but before measurement. 

c 1 

A negative value is predicted which has no meaning. 
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ACTUAL SOIL KOISTUK 



PREDICTED SOIL MOISTURE (>£) 


Actual and Predicted Soil Moisture Percentages at 0 
to 15 cm. 


Fig. 5. 
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Soil Moisture Model 

The original wheat crop coefficient curve developed by Jensen al. 
(1971) was: 

Y = 0.233 - 0.0114X + 0.000484X 2 - 0.00000289X 3 (9) 

Y = 1.022 + 0.00853D - 0.000726D 2 + 0.00000444D 3 (10) 

where : 

Y *= Wheat crop coefficient 
X = Percent of crop cover 

D ■ Days after 100 percent crop cover. 

Equations 9 and 10 in conjunction with climatic data (Table 8) and soil 
moisture information (Tables 5, 6 and 9), were used in the computer model 
developed by Jensen et al. (1971). The soil moisture depletion for both 
fields in most cases was overestimated (Table 10) . 

Regression analysis of leaf area index data for field A (Fig. 6) was 
used as the new winter wheat crop coefficient curve (Fig. 7). The equations 
of the curve were: 

Y = 0.005 + 0.0165.X - 0.000467X 2 + 0.00000402X 3 (11) 

Y = 0.998 - 0.00297D - 0.000747D 2 (12) 

where: 

Y “ Wheat crop coefficient 
X 55 Percent of crop cover 

D D Days after 100 percent crop cover. 

Fig, 8 and Table 10 represent the results from the computer model with 
equations 11 and 12 on dryland (Field A) compared to the actual measured 
values. The' actual soil moisture values compared very closely with predicted 
values of the model until near maturity of the wheat crop on June 19. At 
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Table 8. Climatic Data. 


Month 

Day 

Minimum 

Temp. 

Cf) 

Maximum 

Temp. 

<'F) 

Solar 
Radiation 
(cal/cm 2 day) 

Dew Point 
Temp . 
CF) 

Wind 

Run 

(miles /day) 

Rainfall 

(Inches) 

March 

20 

26 

52 

561.5 

28 

144 



21 

31 

53 

492.6 

36 

113 



22 

39 

56 

505.9 

36 

327 



23 

39 

63 

33.8 

39 

157 

1.10 


24 

36 

49 

108.2 

36 

16? 



25 

35 

40 

90.2 

35 

415 



26 

35 

42 

468.6 

35 

284 



27 

35 

55 

205.5 

35 

127 

.70 


28 

42 

51 

163.4 

42 

200 



29 

32 

50 

91.1 

32 

160 



30 

32 

39 

47.7 

32 

166 

1.00 


31 

32 

37 

214.9 

32 

325 


April 

1 

30 

45 

588.5 

30 

239 



2 

30 

58 

428.7 

34 

79 



3 

31 

54 

429.0 

32 

187 



4 

31 

48 

642.3 

31 

274 



S 

22 

54 

634.0 

30 

163 



6 

31 

63 

627.3 

28 

164 



7 

37 

70 

44.7 

35 

124 

0.25 


8 

24 

37 

381.6 

24 

378 



9 

17 

33 

596.2 

17 

262 



10 

19 

35 

664.7 

19 

192 



11 

26 

53 

625.2 

37 

123 



12 

31 

64 

400.6 

38 

68 



13 

37 

62 

595.8 

40 

102 



14 

46 

66 

470.8 

57 

273 



15 

58 

78 

216.5. 

58 

387 



16 

25 

61 

642.3 

31 

220 



17 

35 

60 

652.1 

39 

158 



18 

46 

76 

643.9 

48 

209 



19 

45 

77 

596.8 

42 

336 



20 

36 

60 

693.4 

23 

219 



21 

38 

72 

672.7 

37 

259 



22 

36 

65 

612.9 

38 

98 



23 

33 

67 

655.0 

42 

78 



24 

46 

73 

162.5 

50 

117 

0,80 


25 

44 

57 

107.9 

45 

115 



26 

34 

48 

221.9 

36 

181 



27 

31 

48 

200.5 

36 

122 



28 

38 

66 

666.9 

43 

153 



29 

45 

82 

635.9 

45 

167 



30 

49 

78 

368.8 

49 

134 
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Month 


Hay 


June 


Minimum Maximum Solar Dev Point Wind Rainfall 

Day Temp* Temp* Radiation Temp- Run (Inches) 

(«F) (°F) (cal/cm 2 day) <°F) (miles/day) 


1 

44 

73 

156,6 

45 

176 

2 

35 

46 

633.8 

37 

178 

3 

32 

58 

704.2 

36 

90 

4 

40 

71 

688.0 

41 

165 

5 

50 

79 

503.6 

45 

319 

6 

47 

79 

702.9 

47 

207 

7 

48 

77 

520.8 

50 

185 

8 

42 

68 

682.7 

44 

160 

9 

48 

79 

706.0 

45 

109 

lb 

44 

77 

698.4 

45 

106 

11 

50 

77 

681.7 

48 

140 

12 

46 

70 

674.5 

40 

144 

13 

42 

68 

672.1 

37 

67 

14 

38 

65 

728.4 

42 

59 

15 

38 

66 

727.4 

38 

84 

16 

45 

78 

718.4 

39 

123 

17 

42 

71 

568.8 

42 

127 

18 

48 

88 

708.2 

46 

77 

19 

54 

87 

705,3 

49 

102 

20 

54 

84 

633.4 

52 

109 

21 

57 

86 

689.5 

61 

201 

22 

51 

85 

611.9 

50 

148 

23 

51 

68 

672.9 

53 

79 

24 

54 

80 

738.7 

50 

79 

25 

47 

72 

641.0 

47 

143 

26 

55 

80 

488.7 

56 

249 

27 

46 

68 

107.2 

37 

266 

28 

50 

53 

624.3 

48 

490 

29 

40 

73 

674.8 

42 

208 

30 

46 

70 

406.9 

44 

125 

31 

42 

62 

751.1 

44 

43 

1 

48 

77 

623.8 

60 

133 

2 

57 

82 

659.7 

56 

247 

3 

53 

86 

645.8 

53 

192 

4 

53 

80 

599,6 

54 

94 

5 

47 

68 

667.5 

48 

133 

6 

50 

79 

736.5 

46 

70 

7 

51 

88 

729.0 

49 

78 

8 

56 

94 

719.4 

53 

87 

9 

57 

97 

739.9 

56 

98 

10 

60 

92 

734.1 

58 

182 

11 

62 

90 

707.2 

59 

277 
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Tub le 8. Continued. 


Month 

Day 

Minimum 

Tetnp. 

CF) 

Maximum 

Temp. 

C?) 

So lor 
Radiation 
(cal/crc" day) 

Dew Point 
Temp. 
(°F) 

Wind 

Run 

(milea/day) 

Rainfall 

(inches) 

June 

12 

64 

91 

498.0 

60 

210 



13 

64 

82 

627.2 

66 

102 



15 

59 

89 

737.1 

52 

216 



16 

57 

94 

743.0 

46 

215 



17 

53 

84 

740.3 

49 

121 



18 

48 

95 

757.4 

32 

239 



19 

54 

78 

695.6 

35 

170 



20 

45 

79 

738.7 

41 

113 



21 

52 

85 

663.0 

54 

85 



22 

55 

86 

725.5 

51 

62 



23 

57 

91 

723.3 

51 

84 



24 

64 

98 

729.4 

46 

126 



25 

61 

98 

663.5 

49 

186 



26 

63 

101 

701.2 

51 

130 



27 

62 

102 

690.7 

50 

156 



28 

63 

93 

594.5 

61 

106 

0.90 


29 

64 

87 

646.0 

66 

97 



30 

65 

88 

647.8 

68 

82 


July 

1 

66 

94 

668.4 

68 

146 



2 

70 

102 

613.5 

63 

197 



3 

67 

92 

639.6 

63 

75 



4 

66 

102 

661,3 

65 

160 



5 

62 

95 

702.4 

62 

94 



6 

65 

97 

715.1 

62 

102 



7 

68 

101 

714.1 

64 

170 




Table 9. Soil Moisture Information.* 


Depth 

(cm) 

Field 

Capacity 

.(%) 

Permanent 
Wilting Point 
(%) 

Bulk 

Density 

(gm/cnr 

0-30 

28.5 

14.5 

1.29 

30-61 

28.0 

14.0 

1.37 

61-91 

27.5 

13.5 

1.39 

91-122 

27.0 

13.0 

1.16 

122-152 

26.5 

12.5 

1.16 

152-183 

26.0 

12.0 

1.16 


*Obtained from the Garden City Experiment Station 



Table 10. Soil Moisture Depletion Using the Model Developed by Jensen et al. 


Date 


Field A 

(cm) 


Field B (cm) 


Actual 3 

Jensen^ 

Q 

Revised 1 

Actual 3 

Jensen^ 

c 

Revised 1 

Revised 2^ 

3/21/73 

17.65 



19.28 




5/14/73 

19.84 

25.07 

19.35 

23.44 

26.14 

20.80 

26.52 

6/1/73 

27.74 

31.24 

27.86 

27.15 

29,24 

26.56 

32.66 

6/19/73 

37.52 

32.16 

34.65 

37.77 

27.61 

31.52 

37.90 

7/7/73 

28.68 

31.29 

34.51 

34.21 

26.75 

31.70 

38.07 


a Actual field measurements of soil moisture depletion. 
^Original wheat crop coefficient suggested by Jensen et al. 
c V7heat crop coefficient using leaf area index of Field A. 
^Wheat crop coefficient using leaf area index of Field B. 
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this date soil moisture depletion was underestimated, but still the difference 
in values were insignificant compared to the available moisture. After the 
June 19 date, comparison became difficult due to the discrepancy of actual 
soil moisture increasing 8.84 cm while rainfall only totaled 2.29 cm. 

Fig. 9 and Table 10 show the results of the irrigated Field (B) using 
equations 11 and 12. The computer model consistently underestimates the 
evapotranspiration. For the time period up to June 1, the differences were 
not significant in relation to the available soil moisture, which included 
an irrigation on May 23 of 3*05 cm. By June 19 the two had considerably 
different values with another unexplained increase of 3.56 cm in soil 
moisture and only 2.29 cm of rainfall. 

Regression analysis was used to develop a third wheat crop coefficient 
curve from the leaf area index of Field B (Fig. 10). The equations for the 
curve were: 

Y = 0.0109X - 0.000288X 2 + 0.00000 33 3X 3 (13) 

Y = 1.52 - 0.000834D 2 (14) 

where: 

Y « Wheat crop coefficient 

X - Percent of crop cover 

D = Days after 100 percent crop cover. 

The computer model's results using equations 13 and 14 indicate that the 
soil moisture depletion was overestimated meaning the crop coefficient used 


was too large. 
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DISCUSSION 

The computer model of irrigation scheduling developed by Jensen ,et al . , 
(1971) uses a crop coefficient which represents the effects of the resistance 
of the water movement from the soil to the evaporating surfaces, the resis- 
tance of the diffusion of water vapor from the surfaces to the atmosphere 
and the amount of available energy compared to the reference crop. The 
model predicts percent effective cover by assuming that it is equal to days 
after planting divided by the days from planting to heading for small grains. 
This proves to be a poor assumption for winter wheat. 

An alternative to this method of crop coefficient determination would 
be the direct use of wheat vegetative growth or more specifically leaf area 
index. If a leaf area index versus the crop coefficient curve was developed, 
vegetative growth would then indicate a specific value for the crop coefficient, 
at a certain point in time. This would eliminate problems due to seasonal 
variation of weather conditions such as an early fall or late spring. 

From this study it appears that a further step can be taken to utilize 
remote sensing. The winter wheat leaf area index has been described, with 
high correlation, by reflectance readings. These readings could be used as 
a direct input into a computer model instead of the original percent of 
effective cover. 

If remote sensing data were available within hours after flight over 
an area, the following procedure might occur. Data direct from the remote 
sensing device would be fed into the computer containing an irrigation 
scheduling model. Meteorological data and a weather forecast for the pre- 
diction period would be the other inputs. From a leaf area index curve 
averaged over many years and the value from the remote sensor, the growth 
of the crop could be estimated for the prediction period. Knowing the 



growth or water use, the computer model would then be able to predict the 
irrigation requirement necessary* This process could be handled by one 
manager for large areas of irrigated wheat land. 
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CONCLUSIONS 

Results from this study indicate: 

1. Vegetative growth was best predicted by a linear relationship 
between leaf area index and the ratio of band 4 to band 5. All 
significant soil moisture effects were cancelled by the ratio. 

2. Soil moisture at a depth of 0 to 15 cm, with specific soil factors, 
was predicted by band 4 and leaf area index with a high regression 
coefficient. 

3. Vegetative growth, measured by leaf area index, was one of the 
necessary inputs in evaluating the winter wheat crop coefficient 


from March to maturity. 
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SUMMARY 

A realization that wise resource management is necessary comes at a 
time when resource use is greater than ever before and the population is 
still increasing. With the use of remote sensing large quantities of data 
are available for ^resource management. These large quantities of data 
have led to the development of automatic recognition techniques in agricul- 
ture. Earth Resources Technology Satellite program provides a system for 
developing and demonstrating the techniques for efficient resource manage- 
ment. 

With the large amount of irrigated land in the world, excess irrigation 
applications means large quantities of water needlessly lost. This valuable 
resource could be better utilized through the use of irrigation scheduling. 
Irrigation scheduling predicts the consumptive use (evapotranspiration) . 

The actual evapotranspiration is dependent upon potential evapotranspiration 
and a crop coefficient which may be predicted by the plant’s actual growth. 
The plant’s growth can be determined by reflection of solar radiation from 
the plant canopy. 

The objectives of this study were to evaluate reflectance for predic- 
tion of soil moisture and vegetative growth; and to determine the feasibility 
of using the plant's actual growth for use in determining the winter wheat 
crop coefficient curve and using it in a computer model developed by Jensen 
et al. (1971). 

The study was conducted on winter wheat fields located northwest of 
Garden City, Kansas. Two soil moisture treatments were used, one dryland 
wheat field and one irrigated wheat field. Both fields were on Ulyssess- 
Richfield siit loam. 
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ERTS-1 satellite passes over any location on the Earth’s surface once 
every 18 days at the same time of day. The satellite contains a line scanning 
device (Multispectral Scanner) that operates in two bands of the visible 
region and two in the near infrared region. Band 4 includes the spectrum 
between 0.5 and 0.6 y, band 5 between 0.6 and 0.7 u, band 6 between 0.7 and 

0.8 y and band 7 between 0.8 and 1.1 U« 

The ground truth data were gathered within one day of the aerial flights 
by ERTS-1. The ground truth data included soil moisture at various depths, 
leaf area index measurements and rainfall readings. The meteorological data 
were from the Garden City Experiment Station with the exception of solar 
radiation which was obtained from the Dodge City Weather Service. 

Stepwise Deletion Multiple Regression (1973) was used to formulate 
equations with the use of reflectance data for vegetative growth and soil 
moisture. The equation that best described the relationship between reflec- 
tance and vegetative growth was: 

LAI » 2.92MSS4/5 - 2.63 , R 2 *» 0.95 (2) 

where: 

LAI = Leaf area index 

MSS4/5 = Ratio of band 4 to band 5 
2 

R = Regression coefficient 

Soil moisture at a depth of 0 to 15 cm was best predicted by; 

SM2 « 101.11 - 4.00MSS4 - 70.33MSS4/5 (8) 

where: 

SM2 = Soil moisture dry weight at 0 to 15 cm (%) 

MSS4 « Band 4 

MSS4/5 ='Ratio of band 4 to band 5. 
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The best winter wheat crop coefficient curve was developed by regression 
analysis on the leaf area index data of the dryland field (A). The crop 
coefficient curve was: 

Y * 0,005 + 0.0165X - 0.000467X 2 - 0.00000402X 3 (11) 

Y « 0,998 - 0.00297D - 0.000747D 2 (12) 

where: 

Y = Wheat crop coefficient 

X = Percent of crop cover 

D = Days after 100 percent crop cover. 

Meteorological data, starting soil moistures and crop coefficient curve were 
used in the computer model by Jensen et al . (1971). From results obtained, 
vegetative growth provides a feasible method for evaluating the winter 
wheat crop coefficient from at least March through maturity. Within the 
limits specified by Jensen et al , (1971), the model and modified coefficient 
proved to be a good estimator of soil moisture. 
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SUGGESTIONS FOR FUTURE RESEARCH 

The research on evapotranspiration modeling and determining the crop 
coefficient by leaf area index should be expanded to include other crops 
and the whole growing season as well as increasing the number of test 
fields. More frequent sampling of soil moisture and leaf area index may 
be helpful. The. neutron probe method for determining soil moisture measure- 
ment would provide a more representative indication due to the increased 
area of sampling. Continued research in using remote sensing for predicting 
vegetative growth with an emphasis on its use as an input in evaluating the 
crop coefficient in an evapotranspiration model may prove beneficial. 

Additional research in the area of detecting soil moistures at depths 
greater than 15 cm with thermal energy could prove productive. 
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Table 11. Computer Model of Evapotransplration by Jensen et al. 
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DATA SUMR t fc'T,DPL/ , 30* , 0.0 ,30*0, 0, 30*0.0/ 

PF AD DATF 

RPAO( 5,14) { DA T £(K ) ,K-1 » 4 ) 

N 4- 1 

F=o.g 

on 100 I =1 , NR EG 

URITFC6 , t 3 ) ( A( T ,.)) , J»i,S) 

13 FORM AT 1 1H1 , • PcGlON: *,5A4,//) 

14 FORMAT < 5X, 1 SA4 ) 

RcAOl 5 , 1 0 ) L L 

10 FORMAT ( 2 5X, 15) 

M=NU) +3 


NN=N( I ) 


DO 100 L=l » LL 

»t40( 5 , 10 ) ijF N 

Rtf AD (5 ,14) (0: SC IK), K* 1,5) 

W 3 I T r (6 ,1 5 ) { f )f. 5 C ( K ) ,K = 1, 5) , ( DA T E { K ) , K = i , 4 ) 

15 FOF MAT ( • 1FAP" : * , 5A4 , 3X, * OATF OF COMtHH AT K'ijs »,4A4,/> 

4 3 I T if (6,16) 

16 POVUT ('0*,T11,'| • |***~* -OIL ^ISTU’fc Of. PUT I C>: 

l IP*. I'. AT IONS I inCML-S 1 * , / , * *‘i T 1 1 , * I 1 , T 1 8 , * ! ' 

n T 3 7 , ’ I ' , T 4 1? , ' | * , T 5 5 , ' 1 * , T 6 4 , *1 IF 

3 1 TO (',/, 

4 i cRfn-rui l coijF | to ^ate I typf-o I oft i r ' u -'- 1 l patf i 

5 I RApJ^O | CAT') I APPLY | ^r, FM no* ) 

24 on no ‘,f=i,.mfm 

if ( p>gpo.FO. f'rTHU no in 1 



I WITH 


LAST* , * 


00 TO .2 

1 no 26 J=l* 1 

26 PI J)=PR«0< I ,.)) 

2 F.KA H5,17nr.l',r.RO|>(l ), CP >*M 2 ) , CROP ( 3) , Mu I* , NDr , l.illl , £ , A VM 
IT FORMAT (5X, 12,2A4,A2,3I5,2F5.2) 

IF ( PHID. >-0 • M T H 1 ) 00 Tn 23‘ 
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Table 

90 

91 

92 
9 } 

94 

95 

96 

97 
9 ? 

99 

100 

101 

102 

103 


104 

105 

106 
107 
109 

109 

110 
111 
112 

113 

114 

115 

116 
117 
1 19 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 
130 

137 

138 

139 

140 

141 

142 

143 

144 

145 


11 4 Continued. 

rtFAO( 5,1 U ( A I ? J1, J=) , 2) ,0P»,M5, (R< J] ,J-4,M ) 

1 13 FORMAT ( 5X ,243 ,^4 .1 , 14 , 1 K ,.2/20'4.2 ) 

RTft'M 5, 19) 0»L( 9- ) , <U 'i Cir ) , ( p - < J ) ,.1=1,3) 

10 FQ3M\T (S e 19.?l 
9'jT'i 22 

23 9::4‘) (5,20) (MRR(J), J= 1 , 2 > ,0»A, N5 , I PR 
20 F'lRNiT (8X,2A3, F4.1 , I4,F4.1 ) 

4?A3 (5,19) nr>L ( NF ) , SU - '-- (NF I 

IF (N5.GC.il « (N5F3 J =f (N5+3M-I 99 
22 C P -4 T I MJ- 
440=0.0 
AKC 1=9. 0 
°C T = 9 .0 
OT =0. 0 

C** J = 4 REPRESENTS F 1 9 ST OA Y OF THE PEFIOD FOR WHICH ANALYSIS IS BEING 

C** RUN 

on 98 J=4,M 
FT! J) =0.0 
FTR=0.0 
RX= MJ) 

SU**(NH- SUMP (NFJ+R (J) 

IF( J-N5-3) 76, 75, 76 

C** DPL AND SUM* ARt SET TO ZERO ON THE DAY HF IRRIGATION 

75 OPL(‘JF)=0.0 
SU^R ( NF ) =0 • 0 

• GO TO 99 

76 IF('ineU )-NDF)109,176, 176 
176 !F(NPO( 1 J-N9H)29, 29, 109 

29 IF|NCP(I1*J-4-iNDE) 30,30,31 

30 PCT = 1 00 . 0 7 * ( NNB< l > +J-4-N& 0 ) /C NTf-MOP ) 

AKCl = C(NCR,mC(HC=,2)-PCT + C(NCr,3)^CT^2^(NCP.,4)vPCT^.3 
I p ( AKCl-i>l( NCR )) 231,232,232 

232 AKC 1=01 ( NCR ) 

231 4V = ( 1 ,0-OPL ( NF ) /AVM) *100.0 
I F ( AV ) 1 30, 131 , 131 

130 «V=0. C 

131 AV3=1 .O+AV 

A4C=AKC1*AL0G( AV3)/AL0G( 101.01 
GO TO 32 

31 DT=NOB( I J+J-4-NDE 
PC T=1 00. 

AV= ( 1 • O-HPL ( NF > / A V M mOO . 0 

AKCi-C.lHCR,5)»r.|HCR,6J*im-C(WCR,TI*l>T**2*CJNM,8»*0T**3 
I F ( AKC l- r >( NCR 1 )9rt,235,2‘i5 
235 I F ( AKC1-01 ( NO R )) 242 , 241 , 241 

241 AKC1 = D1(M*.S) 

GO TO 242 

9.3 AKCl=f;( NCR ) 

242 I F(AV) 233, 234,234 

233 A V = 0 . 0 

234 A V 3= 1 . O-fAV 

A < C * A KC 1 * 4 (. or, ( AV3 ) / ALC'G (101.0) 

32 C T( J)=AKC+X( 16,1 , J) \ 

1F(AKC-FI 38,121,121 

31 JF(F(J-1 1)42,42,43 
43 FTW=0.9>MF-/.k.C)4X( 14, J ,J) 

RC J-l)=r (j-i)-FTr 

jr (4 JJ-D) 49, 121,121 
49 P< 4-I)*RC J-n ♦«.( J-l J 
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146 

147 
14R 
14') 

150 

151 

152 

153 

154 

155 
154 

157 

158 
154 
160 
161 
162 

163 

164 

165 

166 

167 

168 
160 

170 

171 

172 

173 

174 

175 
1 76 
177 
173 

179 

180 

131 

132 
183 
134 

185 

186 
137 
183 

189 

190 

191 
19? 
193 
104 

195 

196 

197 

198 

199 
7 JO 
201 
202 


203 


9 ( J — 1 |*0.0 

45 1 <= P I J-2 > 146.121 , 121 

4 'i I < J-3)aKM-3J+4( J-2) 

L ’ ( J-2 ) = 9 .0 

43 KPIJ- 3)) 33, 121*121 

53 FT7* F TK +R ( J-3 ) 

f- ( J-3 ) =0.0 
GO TO 121 

62 I“(H(J-2))44f 44. 47 

47 £ T 3 =0. 5* (F-.'Kf )*XI16,1 »J) 

8 ( J-2)=4 (J-2J-ETR 
GO TO 45 

44 1 17 ( 7 ( J-3 1 ) 121 , 12 1 ,40 

43 FT?a0.3*(F-AKCl*X(16»!»Jl 
ru J-3PR ( J-3 ) - FTP 
GO Tp 40 

121 I F( STR ) 50, 51*51 

50 ETR^O.O 

51 ETC J) = FTt JI+ETR 

91 DPH NF I -GPL (*\F )*ET(J )— RX 
I F ( OPL ( NF I 1115,99,99 
115 DPL(NF) =0.0 

99 CONTINUE 

FTF.SETINF, JJ=ETR 
FTSET(NF,J1=ETIJ) 

• AKC11 (NF, J) =AKC1 
&KCSETCMF, J)=AKC 
DO 390 NM=1,4 

890 RSETI N'SNF * J) =8 ( J-NM+1 1 
DPLSf T (NFt JI=yPLI NF) 

98 CONTINUE 
S!JM=T = 0.0 
04 57 J=4,K 

57 ,SU02T- SUME T + ET(J) 

R 0 I FsK — 3 
A ET=SU«cT/ROIF 
AcTFin(NF)=AET 
On 880 J =1 * 3 
830 C90PSTINF,J|=CRJP(J) 

M,3l)=NCI? I I)+N I M 

IF (fJD!UII*NI JJ+3-NDF) 250,250,255 
250 PCT=100.CM None I) +?:< I) *2-NDP) /( tvGF-NOP) 

AKC5 = C(:jCR,l)+C(N'CP,2>+PCT + C(f'CF,3>*PCT**2+CINCR»4>* , >CT**3 
GO TO 260 

255 0T=M08( I)+NU 1+3-NOE 
PC T = 1 00 . 0 

A KC5 = C ( NC P » 5 ) + C ( NC R , 6 1 * DT + C ( f R , 7 1 *-n t 4 2 + C I M C F. » 8 ) * 0 T* * 3 
2o0 IF (AXC5 .IT. 0(MC C 1) *» K C 5 = 0 1 G 9 ) 

I* (AKC5 .GT. ni(NCf II AKC5*01( \CF) 

A J JS-NOOI I<J+M{ I 1+3 

IF (AJJ5 .GT. TP(II) Gu T-’l 7034 

OLT^DTIM) 

G-J TO 7 341 
7054 niT=0T2( n 

73 41 FTP5= ( FT API I) / C iiXlM I ( A J J5-T P ( I > ) /OIT ) * * 2 1) J * TCT (I) 

ETAS = AKC 5*cTP5 
OPLA = HPI. I F’P ) 

C** SUBSCRIPT J=1 1 5 20 79 — J= 2 IS 30 *0 — J-3 IS 4070 

C*» ■* J = 4 IS 5T'I) -- J =5 IS 60t0 

NPCTs-IOO.OM I l >+'•:( 1 ) + 2~N r 'Pl/ ( J r E + 33.-';i>P 1 
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2 04 

205 

206 
207 

2 0* 
210 
211 
212 

213 

214 

215 

216 
217 
2U 
210 
220 

221 

222 
223 
2 24 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 
2 37 
230 

239 

240 

241 

242 

243 


244 

245 

246 

247 
243 
249 
? 50 
251 


I e (SPOT- 100) 243,243,249 
?49 mpCT* 100.0 
243 C«WT1NIJP 

DO 108 J = 1 , 5 


RJJ=J*1 

OPAK$U< J J»i\PCT*AV w *K 

l D r= (j*i i *io 

AVWsC'PAKS'H j)-c>pl(nf ) 

CALL SCMtD' ( *‘!?9, 4 VW , M'.H, N XT , NXDD» I t fPLA , 4V M » ) 
CALL CAT f £ < NXP, l X * I Y» NDH ) 

CALL OAT p F INXOP, JX,JY,NUHI 

59 IF (OPAKSIUJ) - DPL ( NT j ) 60,61*61 

60 A T 9 = l)PU NFJ/E 


GO TO 63 


61 AIR = OPAK$U(J)/E 

63 WRITE (6,64) * CROP, AKC5, OPL (MF) ,0P AKPUC J I . ETA5 , AIPR f PCM I X),1 V.MHM 

1 ( JX ) , JY , A I R , !»l,NF 

64 FORMAT < • 0 • ,2A4,A2»F5.2,F9.2*5X, *20* P',2F9.2,' I • ,2A3,2 (2X , A4, I 
23),* I • ,F4.1, 17,214) 

00 TO 108 

65 WRITE (6*6 3) IPC ,0P AKSO ( J ) * tTA 5 » A IF ft , MOM ! X ) , I Y* ■ v SO:>J < JXJ ♦ JY , A Ik 

63 FORMAT (' '»T31»12f*7, O’ » 2 F 9 . 2 , ' I ' » 2 A 3, 2 ( 2X , A4 , 1 3 I , 


1 » | • ,F4. 1,17, 214) 

103 CONTINUE 
109 CONTINUE 

W1 ( 1 » N4 I =D PL ( NF > 

W 1 ( 2 ♦ N 4 ) = S U MR ( NF ) 
WH3.N4) =K( M-2 1 


Wi(4,N4) =R(M-1) 


N 1 { 5 , N4) =R ( M) 

N4=M4«-1 
110 CONTINUE 

WK = ( N£8( I) + N t I » - 531/7 

op = 14.*(5(l ,11+fM !,2)*WK + ft( I,2)+WK**2+ P(I,4)*WX**3 + 

1 5(1*5) *w'K°*4 + D( I , 6) *UK**5> 

I p (PP * LT . 0.0) PP = 0.0 

WRITE (6*163) P^tlfL _ „ 

163 FORMAT ( * OPRORA OLfc RAIN NEXT TWO WEEKS = ' » F 5. 2 » 2X * ' I NC HE S 1 * 30X » fc 

1 ) 

WRITE < 6,8.01) 

801 FORMAT ( »-**-*T Atilt OF OAKY VALUES***') 

00 830 NF=1,NFN 

WRITE (6,303! (CRnPST(MF,K),K=l, 3) 

803 FORMAT ('O' ,2A4,A2,/, 

X « 0 DAY FTP ET 20 AKC1 AKC', 

2 PX R(J-l) R ( J-2 ) '» (J-31 DPI'*/) 

00 R20 J=4 , v 

WRIT'- (6,302) f!D ( I * J ) * ( TRSvT I NF , J) , i: TS t T ( N F , J ) ,X ( 16 , I * J ) , 

1AKC11 ( MF , J) * AKC S r T ( NF , J ) , ( f . S E T ( NM , JF , J I * f' M - 1 » 4 ) , DPI S ET ( f F * J ) 

802 FORMAT (' ' , 15, 2X, F 3, 4,f.a.3,dF 0. 2) 

820 CONTINUF 

WRITE (6,321) AETFUMMF) 

021 FORMAT ( 13X, * Af T = • , F 7 . 3 I 
P30 CONTINUE 
100 CONTINUE 


2 


t 


C 

C** 100(1)= NO. i*F FlEmr. FCT WHICH ANALYSIS WAS »UN 

•lATJ 1 ) *N4” 1 
81 RETURN 


252 

253 
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Table 


2 54 

255 

256 

257 


258 

259 

260 
26 1 
262 

263 

264 

265 

266 


267 

268 

269 

270 

271 

272 

273 

2 74 

275 

276 

277 
273 
2 79 
280 
281 

232 

233 

284 

285 

286 

207 

283 

239 

290 

291 
2 92 
2 97 

294 

295 

296 


11* Continued . 


end 


c 


SHRRiliJTTNF EVAP(1*K) 

SUBROUTINE Ti.l CALCULATE EVAP^TKANSPP AT ITM POTENTIAL 
- A L A* 5 T H 1 

C'lMM.IM A ( 4 1 5 ) * CTRI4 ). "XF (4) ,* '•( 4»30> , 

IX (16, 4,30) ,ntf,C<5> ,f'AT*' (4 ) ,Crr.P t 3) , AI cl -‘ ( ?) , F« PC (15 ) , 

2N ( 4) ,Ni)R{4) , ” S n ( 4 ) » "f'/.Y ( 4 ) , '•! 1 ( 3 * 1O0 ) * I.’ ( J » 6 ) , 9 !■ ( S ( 4 i ; 01 , r ( i0 ) 


09 10 J=4,K 
X ( o , ! , J 1 = (X(1,I,.1) 
15 X ( 7 , 1 , J ) - CTR 1 1 1 *< 
10 C 3NTTNUE 
RETURN 
END 


+ XI?,T, j>)/2.0 

Xto, I .J1-TXM I ) ) l 3 T I * J ) * 0.000673 


C 


SUBROUTINE VAP.3RU ?K,CW) 

SUB S ‘ , U T I \’ E TO CALCULATE HEAT FLMX.EQ POTENTIAL? NET RADIATION 
REAL METH1 

C OM AON A ( 4 ? 5 ) » CTR<4),TXM4},N0(4,30) , 

1X1 16 ? 4 t 70) , QESC I 5) » n ATh (4 1 , CRO^ H ) , A1 RR( 21 ? Ff'RC (15)?^ 

2NI 41 .NO 1 ? 14) »RSO( 4> .MflDAYl 4) ,Kt( 5 , 100) ,C< 8, 3) , P.REG.(4»30) , R f 30) 

COMMON /NEW/ W ( 4 ) ? MON (13)? I D,NCR » NDE ,N ; JP , R ( 4 » 6 1 , c TAP ( 4 ) r TP{4) , 
10TK4) ,DT2(4) ,FCT(4) ,ETP5 
DIMENSION C W I 4 1 
on 30 J=4,K 

• IF(X(4,I,J) .EQ.OIGOTO 35 

X(3,T,J)= X( 5. I ? JJ/24.0 __ 

VPS1= -0.6959+0. 2946 *Xl 2, ! ?J)-Q. 005 195* X (2. I ♦ J 1 **2+0. 000009* 


IX( 2, I , J)**3 

VPS2= -0.6959+0. 2946*X I 1, I , J I -n.005195* XU t I t J J **2 + 0. 0000*9* 
1X11,1, J 1**3 

X ( 9? I ? J ) * IVPS1+VPS2) /2.C 

XI 10,1 ,.l)=-0. 6959+ 0.2946*X(4,1 , J 1 -0. 005 195 *X ( 4, 1 , J)**2 +' 
10.000039* X(4? I « J 1**3 

X( 11 , 1 , J 1 ■= ( X ( 6 ? I ,J)-CX(1, I , J — 1 >+X(2,l , J-l )*X( l,I,J-2) + X(2,I,J 
1 — 2) + X{l,IfJ— 3)+X(2,I ,J— 3))/ 6.0)*5 
Tl= 0.041 + 0. 0125+Xl 6, I » J)-4.534*X<6, i , J) **2/10**5 
T2= 0.959 -0.0125=*X( 6, I , J 1+4. 534*X(6, 1 , J1**2/10**5 
X (12 i 1 , J ) = ((Xlltl.J 1-321/1. 3 + 2731/100.0 
X ( 13 ? T ,J1= 1 1 X ( 2 » I , J 1 -32 1/1.8 + 2731/100.0 
Y= XI 10, I ,J) 

JJ=ND&( II + J - 4 

EMT=0.3 25+ 0. 045* S INI 3 0^- ( .1 J/ 30.-1. 51*3. 1416/1 ?0. 1 
X ( 14.) , J 1 — I E NT -0. C44*S OPT ( Y 1 1 * 1 1. 719 ( x I l 3 , I , J 1**4+ X { 12 , I , J I 
1**4 1*0. 5 

X ( 1 5? I ? J 1 = 0.77*XI3,I.JI-I1.22* X (3, 1, J > /r SC (1 1-0.13 )*X( 14, 1 , I ) 
30 X( 16? I , J 1 = ( T 1 * ( XI 15,1 r l)-x(ll,I ,J> I +T2*15.36*(.75+CW( 1 1 * 
lX<5»I»J>) f: <X(9,l . J I —X t 10 , I , J 1 ) 1*0.0036 7 3 
A J J5*ND? U ) +M( I 1 +3 
IF (AJ.I5 .DT. TP (1)1 CO TD 34 

niT*onm 

03 T3 341 
34 0 L T » 3T 2 1 1 1 

341 f TP5= ( El APU 1/ ( EXI’ III A JJ5-T P| I) 1/TLT ) **2M)*FCT< I) 

3 5 PE TURN 
END 


SUPRDtJT 1 WF PP I NT R ( l ? K J 
C. SI*8«3UTINE TO PRINT REMPNAL DAT A 
REAL NET HI 
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2j 7 


298 

299 
390 
3 31 

302 

303 

304 

305 

306 
307 

30a 

309 

310 

311 

312 

313 

314 

315 

316 


317 

313 

319 

320 

321 

322 

323 

324 

325 

326 
3?7 
3 28 

329 

330 

331 

332 

333 

334 

335 

336 

337 
3 38 

339 

340 

341 

342 

343 

344 
34 5 
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Ci'SMnu M4,Mt CT9<4) ,TV- (4) ♦* . < 4,30) ♦ . , , . k 

1 X < 1 6 1 4 •*' o 1 , 'H SC i « ) ♦ i) AT • ( * I ,C -V ( » > , A I ' ' < 2 J , rr NC iio ) t 
2» ( M I , NN< : ( 1 1 , « S'U 4 ) , " >!’ '.Y ( 4 I , v. 1 ( f , 100) *L < :< I o ) i r R >. >• ( A , „ 0 ) » •• ( .01 
C V«yo\* /I'C .</ VM A) ,M<.> 1 1">) » ! C,UCdM0e. fi'UP.H (4,6) * t TAP ( A 1 , T » ( A ) , 
lDTlC A) t OT2( A) ,PCT< Al ,tlf5 
j.mvopi i ) 

CAM. r.ATCfc l JJ f M?J,-Nll)» 33o) 

W K I T f ( A » 1 0 1 ( A( I , J1 , .1=1 , *3 I , V GN< Y‘!1 t f;l r: 

10 FORMAT ( 1 H- i 5 X * * R f- G I l) : * » 5 A A * 5 X » 1 a f GINNING O AT E = * t A A , 1 3 J 
WR I TF ( 6 » 1 5 J 

15 Fjt R?1 AT C 1 H- f • DAY TAVG PS UA VPS VPU Pn G 

1 FTP FQ* ) 

WRITE 16, 27) 

27 Fi'JKMATUH J 
rvt 20 J =4 , K 

1 , X ( 1 5 , I»JJ» X ( 1 1 , l * J I , X ( 7 , I , J 1 < X ( 1 6 , I » J ) 

20 CONTINUE 
35 WRtTE16,40) FTP5 

AO FORMAT 1 1H * • TOP ? C AST * POTENTIAL ET NEXT 5 PA YS= * * F5 . ?. ) 

25 FORMAT C 1H ,I5»F7.1,F6.0»F6.1,F7.1,Fn.l,F7.0,FR.l t 2F7.2) 

RFTIJPN 

END 


S(IH°PIJTIN£ ETAVGC I I , FT A , MOD , 1 , 0 , 01 . tVP, DPI. ) 

•COMMON M 4 * 5 1 * CTR( 4 I , TXR 14 ) »NN ( 4 , 30 ) , 

1 X ( 16, A, 30), 0= SC I 5 >» 0 ATFI A) ,CR0r>(3) ,AIPP.(L) ,FNRC(.15), 

2N( A) , NO-31 A) , R SO ( 4 ) , W 0C AY ( A ) t Wi 1 5, 100) , C t 8* d ) t Rr -EG < A » 30 J , P. I 30 I 
COMMON /NEW/ W 1 4 ) , ‘dN 1 1 3 ) * 1 0 » NCR , NOr. ,NDR » P (4 ,6 J » t TAP ( 4 ) , T P ( 4 ) , 
1DT11 A) , IiT2( A) ,FCT(A) »CTP5 
DIMENSION DC 01 ,01 181 
A I = 1 1 

AV-t 1 .0-DPL/AVM)n00. 

IF (AV .GT. 0.0) -GO TO 200 
AV=O.C 

300 AV3® 1 +AV 

5 I? (II iGT. NOE) GO TO 2 
AP=NOP 
Afc=NDt 

PCT=100.«( AI-APJ / (A'-AP) 

AKC1 = C(NCP.,1)+CINCR, 2)*PCT+C<NCR.3)*PCT**2«-nNCK»4|*PCT**3 
GO TO 1 
2 OT-II-NOE 

A K C 1 = C ( OCR, 5)+C(NCF ,c)*DT+C( NCR»7)*OT**2 fC<NCR* 3)*OT< *3 
1 IF ( A K C 1 .LT. C'l’lCR)) AK C 1 = 0 ( NC R > 

IF (AKC1 . G T , OinCR.n AK.Cl = Di(NOU 
1FUI ,r.T. TP(II) GC TO 7 
DLT=OTim 
GO TO 3 
7 OLT=OT2( I) 

M A K C- t K.Cl*'LOG( AV3 1 /AU’0< 101.01 

fc T Aai,K C ‘ ( F7AP( ))/( RX<> ( ( I AI-Td 1 ) ) AM T 1**2) U 

IF III-M.O .IT. 5) E T A = L T A* FC T t T J 

RETURN 

FND 

SOARiiU 7 ( N£ OAT F £ ( II , '*0, II i) . N')M ) 

C CJlLC'll ATT', vr),NTH AN n DAY F’J'Vl JULIAN GAY 

DMFciSK'O • (121 

nata N’jn/o, 31,40*01 , a 2 1 , 15'2, 1R2»213» ZAA»2 74,205 ,335/ 
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34 6 
347 
34? 
349 
35') 

351 

352 

353 
3 54 

355 

356 

357 
35? 

3 59 

360 

361 

362 


363 
3 64 

365 

366 

367 
363 

369 

370 

371 

372 

373 

374 


375 

376 

377 
373 
379 
3dO 
3 31 
3 32 
333 
3B4 

385 

386 

387 
3 8fi 
389 

39n 
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m io J=7,12 

TP (IT .U . N*i>C J H 91 TfJ 12 
10 C )•'! I NIJC 

j*n 

12 MN=J-1 

110=1 I-NNOU-U 

If III .17 . NOH ) 00 TO 14 

M^n 
1 1 n = o 

14 ffTURW 
EiND 

CliSo') C /”.l/ WI4I .K.’.'CS.’nt.NOC.iK.ol.tUI'l'.l .TPC.I . 

1 'OCT ' t 'NOV' , 'DEC ' * 'NoNE '/ 
end 

SUSP f> LIT 1 NF SCHEO ( MBIT , AVU. NDH *f!X D . NX DP* 1 * D p L t AVMt D t 01 1 
COMMON 4(4,5), CTR(4>,TX<M4>,ND(4,30|. 

IX ( lu.4, aO» , HE $Ct 51,0 ATT (4) .CROP ( 3) , AIRRI2I , F n «C (1 5). 

COMMON / Hfc -if W ( 4 ) » i'TN ( 1 3 ) » I 0 ,NCR , NOti , NOP .B < *+ , 6 > . E TAP ( 4 ) * Tpm , 
IDT 1(4) »0T2l 4 1 , FCTt 41 .ETP5 

- ’ CH(-;rK S TO ^SEE^l^THE 1 FlrlO NEEDS IRRIGATING A T BEGINNING OP DAY 

IF (AVW.LE.0.01 GO TO 10 

C CALCUl AT liJG ESTIMATED DATE OF IRRIGATION WITHOUT PFOB RPf.ClP 

DO 1 I 1 a!/RD, NOH 

CALL ETAVG U I .FT A , MGO. 1 . 0 . n l , AVM , DPL ) 

AVW= A VW“ ETA 

.IP I A V <4 • L E • 0 . 0 ) nr TO 2 

1 C TNT I NHF 

c IF IRRIGATION IS NOT REQUIRED BEFORE HARVEST 

GO TO 12 

2 NXD=T I 
NXDP=NXD 

C CHECK IF RAINFALL PROBABILITY IS TO Bf USED 

C B ( I,1) = 0 IF RAINFALL PROBABILITY IS NOT 0E$1»E0 

IF | ABS ( 0 ( l , 1 J ) .LT. 0.DD001) G«J TO 11 
£ OETEPMPIL' NUMBER OF DAYS F0° EXPECTED PRt;C l p l TAT I JN 

WK={MBD-53)/7 
IS A I ^ I I 
T=A1-BD 
BtisMO+T 

IF IT ,LP. 14. ) GO TO 15 
rt DsBD-T +14. 

T = 14. 

13 6\M=*AVW*P AMT<TtW^tI) 

L = I I ^ 1 

O'l 3 ?I*L**Jr)H 

CALL ETAVG I I I ,f TA.UBC, 1 ,0.U! . AVM,D»L) 

AVV»»4VH-f TA 

IF CAVW . L E . 0.0 1 GC TO 4 

3 CJNTFf.'UF 

C IP?. I GAT I ON UOT REOUl'RFU BEFORE HA^Vi-St 

C ^CHECKING IF EACH F X T R • r D J ~ ? I G A T I ON -ATE USING PPORA", J t. I T 1 1 S 
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Table 11. Continued. 


( 


(. 

f 


f 


346 

347 
34* 
349 
300 
351 
353 
353 
3 54 
3 50 
356 


no io j = ? 1 1 z 

! s (M .ilr# J )) V T ri 12 

10 C i m T1MJL- 

j=n 

12 M N - J — 1 

iio= i i-uno< j-i ) 

IF (IT .IT. MDH) 00 TO 14 

l in = o 

14 P-r= TU*5W 


357 

35? 

359 

360 


nine*, data 

COMMON / \ u >J / W ( 4 ) i M0N{ l 3 1 i I D t M r 'E: ♦ MOP 

lUTl(4) f DT2(4), FCTl4)fFTr>5 

DAT M ON /•JAM 1 * 1 F r 15 1 • 1 K 1 i ‘APRS 1 w AY f t 

1 *L«CT • t 1 MOV 1 T 1 DcC 1 f •NitWS 1 / 

EMM 


f JK,6)»tTAP('tJ,TPC 

• JU:^ ■ , 1 JUL 1 t i 1 


4) t 


, c P f 


t 


361 

362 


363 
3 64 

C 

365 

366 

C 

367 

363 

369 

370 

371 

C 

372 
3 73 

374 

C 

c 

375 

C 

376 

377 

373 
379 
3dO 
3 ?3 1 

332 

333 

334 
385 
336 
387 
33R 
3B n 


SUBROUTINE; SCHEf) ( MhO f AV-^t NDM t NXD f NX(?P V I f f JPU. t 01 ) 

COMMON M 4 , 5 ) t CT^(4)tTX r M4> f Nn(4,30N . 

IX I 16t 4 i 301 t 03SC ( 0 ) ♦ DAT r ( 4 ) , C K-P O) t M RM 2 > t F^RC 1 15 > » 

2MU > , NUB (4 ) » ".SH( -V ) , vrjf.Av ( 4 ) , Wl < 5 , 100) ,C ( 8 , 3 ) * P.Wl 4 , *0) , ? *30 
CQMMOV /'-it a/ W<4) ,0‘O"<13) ,lf>,NCft,ri0tt'-inn,K(H,6J,cT4PtH) ,Tf (h) , 
10T1(4) ,nT2(4),?CT('+1,ETP5 

CH!*:f.K S Tn 'sEE 1 f°tMP FIfLD NEEDS IRRIGATING A T BEGINNING OF DAY 
IF ( AW.LE.O.O) GG 7 0 10' 

H0= M 3D 

CALCULATING FSTI MATEO DATE OF IRRIGATION WITHOUT PROH P r &Cl® 
on 1 H=w«o,NOH 

CALL ETAVG U I ,FT A, MBD, I ,0,01 ♦ AVM.DPL) 

AVW=AVW-£TA 

.1 F I a V :1. IE. 0 . 0 ) nr TO 2 

1 CONTINUE 

IF AN IRRIGATION IS NOT SPOUIREC BLFORS: HARVEST 
GO TO 12 


2 NXO= I I 
NXDP=NXO 

CHECK IF RAINFALL PROBABILITY IS 
B( 1,1 1 = 0 IF RAINFALL 
IF ( A R S ( :‘< ( I , 1 1 ) 

DETERMINE 
WK=( HBrj-53 > / 7 
IS A I = ] I 
T=AT-FO 
P.f) = *D + T 

IF <T . L fl. 14. 1 
HOsGil-T H4 . 

T = 1 4 . 

15 AVH = AVl->nAMT( T,WK, I ) 

L = I 1 + 1 

O'l 3 T 1 =L* “1DH 
CALL ETAVGUI »f TA,’1BC. ! 

A V Vi = A V w “ E T A* 

IF (AW) .U. 0.0 I GO 7 1 1 4 
0 CONTINUE 

irrigation mot fttcuisro nr-FOR-- 

•Gil TO 13 


in ef us c o 

IS MuT OESl»EO 
TO 11 

PRECIPITATION 



, AV M , L) P L ) 


HARVEST 


CHECKING IF EACH FXTC"''FM IRRIGATION O-iTf USING PPf'B AG II I T I : S 


PROBABILITY 
.LT. 0.000011 00 
NUMBER OF DAYS FP° EXPECT? 0 


GO TO 15 


, 0,01 


non 


D-67 


Table 11* Continued* 


301 

392 

393 

394 

393 

306 

307 

3on 
309 
409 
4 31 

402 

403 

404 

405 

406 

407 

408 

409 


410 

411 

412 

413 

414 
*15 

416 

417 

418 

419 

420 

421 


423 
4 24 

425 

426 

427 

42a 

429 

430 

431 

432 


C Mf I.- A 1 J ; * r> SUITS T '-' f=u> r M r > r: XT C K$ 1 1 Ni I; s M0ATlt T \ Pz*\]0" 

4 i* I1H •i’9. uxmn f >:. TO n 
WK - UK 4*1/7 
no s ii 

C SITUATK'M UHiiRt FI FLO NFF9S I f:P I r » ATI 0*1 AT THr h LulUUlt.fi DA f 

10 UXC-M80 
NXO^NXO 
^ n TO 11 

c SITUATION WHfcRF AN' IRK 1 GATI 0 N IS MOT KCWnCD 'ir.HW M.'^VJST 

12 

13 *.!<nP=Nf)H 

11 RETURN 
END 

■ funct ion °‘ :j it(t, wK| i > 

C F-J'CT I ON FOR PFDBABLF PRECIPITATION 

CJKMnN /TEW/ Ml O.^NIllI ,IP,NCR.Nn=,N'JP.a(4»6l .tTAP(A) ,TP(4) , 

IDT l ( 4 ) , 0T2 ( 4) t c CT (4) , F T^S 
PA^T =T*(SU .1)+KM,2)*WK+P( !i3)*WK.*WX + ft(I #41- 

1 V.'K**3 + R( I » 5 ) *WK.**4 + P.( 1 , 6 ) *W<**5 ) 

RETURN 

ENO 

SUP. ROUT INF PRINTS ( NR e G t METH1 , RNRP ) 

C SUBROUTINE TO RETAIN INFOf M AT I n N IN "SAVE" FCR NEXT PUN 
REAL »‘LT HI 

CQvsvjOK A { 4 ♦ 5 ) * CTR(4),TXP(41,N0(4,30), 

1 X( 16 » 4»30) ,C r :SC( 5 ) » 0 ATF (4 ) iCROP ( S I , Alf R ( 2 ) * FC’r.C ( 15) * 

2N(4) . t-;DKm , - S 0 ( 4 ) ,M’.ifAY( 4 ) ,'•'!( 5,100) ,C( 0, 3) ,RKcG,4, r>0) ,F ( rOI 
C0'1'«f N /HEW/ W ( 4 ) f WON (13) , I D , NCR , NOE, NOP , r i ( 4, 6 ) , (-TAP (4 ) * i «’( 4 ) , 
10T1(4),OT2(4) ,PCT<4) *E TP5 
WRITE (ft, 11) NOB ( 1 1 

11 FORMAT ( lHl * ' NO. OF FIELDS ='*I5) 

00 40 1 = 1 i NR EG 
Kl*l«( I ) + 3 
K=K 1-2 

IF ( R NftO . C 0 . M t T H i ) GO TO 15 

WRITE (7,10) (N0( 1 , J ) , X( 1 , I , J ) ,X (2, 1 , J) *X( , I , .) ) , X(4, I , J) , 

IX ( 5 , 1 ,J) , J = K , X 1 ) / 

VI R \ T r (6*10 \ (UP( IfJJ iXlli ItJ) »X(2|IiJ)iX(3t ItJ) fX(4f.IfJl» 

IX ( 5,1 ,J> ,J=K,R1) 

10 FiJR*'AT (5 X, I 5,5F5.0) 

G 1 TO 40 

15 WRITE (7,20) (ND l I , J ) ,X( 1 , I » J) , X(£, I * J I ,X( 3i I , J) , X(4, I , J ) *X( 5, l , J « 
i ,RPf.('.< 1 , J) , J=K,K 1 ) 

;jo j Tp ( o , 20 ) ( MO ( I , J ) , X( 1 , I , J 5 , X ( l , I , J ) , X ( 3 , I , J ! , X< 4 ■ i , J 1 , X ( > , I , J 
1 , RP I G ( 1 , J ) , J=K ,K 1 ) 

^0 CONTINUE 

20 FORMAT (5X, I5,5F5.0,r5.2) 

K.*M0P. ( 1 ) 

PI *»0 J*1,K 

i.i p j t r- ( 7 , 5 '3 | w U 1 , J I , N i ( 2 , J ) » ■« 1 1 ? , »U , '•> 1 ( 4 , .) ) , w 1 ( 5 , .1 ) 

WR IT t. (6,55) Will , .1 1 ,Ai ( ^ , J I ,H 1 ( 3 , J) , -i 1( 4 , J i , •1 ( 5, J ) 

50 CONT I f IJI 
55 FOpf'Mt 5F10.2 ) 

RETURN 
END 



♦ENTRY 
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ABSTRACT 

Wise resource management techniques are necessary if the population 
of the Earth is to continue to expand. The Earth Resources Technology 
Satellite program combines remote sensing in space with efficient resource 
management. Water is a valuable resource needlessly lost by excessive irri- 
gation applications. If needless loss of water is to be lessened, deter- 
mination of evapotransplration will be necessary. Actual evapotranspiration 
is dependent upon potential evapotranspiration and a crop coefficient. One 
method of predicting the crop coefficient is to use the plant’s vegetative 
growth which may be determined by reflection from the plant canopy. 

The relationship between soil moisture, vegetative growth and solar reflec- 
tance was studied. Vegetative growth was evaluated by leaf area index with 
the equations 

LAI = 2.92MSS4/5 - 2.63 , R 2 = 0.95 

where: 

LAI = Leaf area index 

MSS4/5 = Ratio of band 4 (0.5-0. 6 p) to band 5 (0.6-0, 7 p) 

2 

R ** Regression coefficient. 

It appears that the ratio eliminated soil moisture effects. At a depth 
of 0 to 15 cm soil moisture was predicted by: 

SM2 « 101.11 - 4.00MSS4 - 70.31MSS4/5 

where: 

SM2 = Soil moisture dry weight at 0 to 15 cm (%) 

MSS4 = Band 4 (0.5-0. 6 p) 

MSS4/5 = Ratio of band 4 (0.5-0. 6 p) to band 5 (0.6-0. 7 p). 
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The equations of the wheat crop coefficient for the evapotranspiration 
model of Jensen and associates s developed by using leaf area index of 
dryland wheat, were; 

Y » 0,005 + 0.0165X - 0.000467X 2 - 0.00000402X 3 

Y » 0,998 - 0.00297D - 0.000747D 2 

where s 

Y ■» Wheat crop coefficient 

X ~ Percent of crop cover 

D ® Days after 100 percent crop cover. 

This method of evaluating the crop coefficient provided reasonable estimates 
of soil moisture depletion. 



Appendix E 


Computer Program to Generate the Mean and Standard 
Deviation for the Interior of a Field 



" DIMENSION LI 

2 D1 MENS TOM 3( 

M F 

01 

(OIO)»NCFtL(010),SLOPE(010),THETA(010),ALINtt010) 

0) . ACH.L ( O10 ) 

3 TNTFGF°-4 DA 

r r - 

? tapf 


ms ton cell ( 7 i 

MS I DM TFRFOT 123,161 

VST ON A M T N ( 1 6 ) ,AMAX(16)*S 
MSI 0*4 RAVF.(16),RST0(16),R 

FIELD/’ ALL 1 / 

0=9 
M = 13 

0 

0 

0=0 

= n 

= 64 
= 64 

0 444,0, NSM, MSWI 

AT< p 5o 2,71,11) 
(5,500,FWn=999) NOATE , MTAPE 

_l!5J 




TOO CDNTIN'IP 

TFINSVII.FQ.l) GO TD 101 

Eg APT 5, 1000, END = 999) SFO, VP T S , (L INFM ) «MCELL( T ),I*1,4),NF!LE,NC, 
1M T APEA,FI C L0 

1000 FDR MAT(I1*T1, 4(14,13) , IX, 1 1 , 1 1, 2X, 1 4, 3X, A3) 

PS IMT 1102 , SEG_, NPTS . ( L T NF (1) . NCEt L (II ,I=l,4),NFILF,NC,NT4PgA,FIFLD 


2 FORM4T(lH , «?FG= 1 , 12, *PTS = ' , I 2, 4( 75, IX, 14, IX) , »FILS= 1 , I 
1I2»*TAP«; = " , 15, 'FTELD=', A4) 

GO TO 202 .. . 


101 CONTI WMF 

READ(5, 1003, END= 999) SEG, NPTS, (LINE (I),NCELL( I) , I = 1 , 8 ) , NF IL E , NC , 
1NTAPFA,F1EL0 ' 


1003 FORMAT m, 1 1,8< I4,T3)f2X, 11,11, IX, 14, IX, 6 3) 

POINT 1101 »SFO» NPTS , ( l INF ( I ) f N0ELL(I ) ,1 = 1, 8) »NFILE,NC,NTAPEA, FIELD 
1 101 FORMAT (1H o >SFG=> , T2. < PTS = 1 ,12,fl(I5,lX,I4,lX), 1 FILF = , ,I2, 1 NC = 1 , 


IT2, 5 TAPe=' ,15,* FIELD* 1 »A4) 

202 CONTINH c 

TF(NTAOFi.NEoNnATE) GO TO 991 


I F ( NO N. c Q. 1 ) GO TO 10 
N0=0 

nn 220 T=1,N T ERM 


Do 19P J=1 ,128 
1 90 IFRF0(J,n=0 
A M I N ( ! 1=128 


AM4X( T ) = 0 = 

S'JM f T ) =0 o 
00 21P .1=1 , NTFP M 
PPODn,-! 1 = 0. 

CONTINUE 

CONTINUE 


IFIS^GoEO.NSEG) GO TO 870 
30 R* AD ( NR F AD ) DATF ,TAP6, IDl , T02 ? L INFT,NCEI.LT 
t F( 0" T r .A!= .VOAT" ) GO T n 990 
IF(TAProNFoNTAPF) GO to 990 
1 F ( I D 1 o NE o N 1 9 1 ) GO TO 990 

I’=H r *2«NF»NID?i GO TO 9«0 





56 

57 
5a 

59 

60 
61 
6 ? 
6?- 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
8? 

83 

84 

85 

86 

87 

88 

89 

90 

91 

“92" 

93 

94 
-7?6- 

96 

97 
99 
99 

100 

Terr 

19? 

103 

TW 

105 

106 
TOT 
105 
1 09 

TT- 
111 
15 7 

TTT 

114 

115 


READ(NOEAO) xdata 
370 CO NT T NUE 

NCN=NC 

NSFG=S c G 

IFdTNFC D-LT. L1NET) GO to 994 

I 8 ( L I NF C 1 > • GT . { JT NET *64 ) 1 GO TO IQ 

NCR*’nCELL8>647< S C G- 1 )«-64*NC.ELL< 1 > til 

PR 1 NT , MC ELL B 9 SE G i NC c Ll (1 1 »NCB 

T F ( NC FLL T . GT .TGI GO T 0 9 9 5 

I F (NC ELI T . L T . ( NCB- 63)1 GO TO 10 
HO 300 T=1,NPTS 

L I N c ( n = LI NE ( 1 1 - L T NCT-O-I 

C’FLL ( I )*64.*NCFLL ( I 1 /77. 

800 CONTINUE 

on 900 t = i qNpt f 

II = 7-0-1 

I F ( I- o FOo NPT 5 ) I T =1 

0X= L I NET T 1 l-L INF ( ! 1 : 

PVl^C FLL (I!) -CELL! I) 
r>Xl = PL0AT(0X) 

ip I OX „ EQ .0 ) GO TO 880 

SLOPE ( I ) =071/0X1 
THETA n ) »ATAM?( OY1 ,0X1) 

GO TO 890 

880 SLOP" ( 1 1 =999. 

THFTM 1 )=ATA\'2( DYl, .0011 

«90 CONTINUE . 

3fn=CFLL(!J-SL0PFUl-LTNE< 1 1 -0 /COS ( THET A ( I 11 
900 CONTINUE 

AAHIN = 100000. 

AAMftX=0. 

00 910 I=1,NPTS 

n*?»i 

T 9 ( I .EO.NPtS) 11=1 
ns = Sl OPE ( T 1-$L0PP(1 I 1 

ALIN? (II = < 8 (IT )- 8U) I/OS 

AC ELL l I 1 = f'8"f T f T* SITOPTI lTD ( I 1 7SLGPF ( II) 1 /OS 

IF{ AAMIN.LT. ALTNEtT )) GO TO 909 

M 1 NT = T 

AAN1N=ALINF( I 1 

909 I F ( A A M A X o L T « A L I NE l T) ) A A M AX= AL I Nr ( I 1 

910 CONTINUE 

LINRTN=IFIX(AAMin«-1. I 
LTN9AX = I FI X ( A A MA X 1 
DM1NI-1 

1 1 -MIN 1*1 

TF(IoEO.O) I=NPTS 
IF( I I .GT.NPTS ) 11*1 

rrrrrTTTTTTnT Tn ctwwt 

TF ( L 1 NNAX.GT .M A XL 1 L INNA X=NAXL 
IF IL IN«T N.GT oL 1NM4X 1 GO TO 2100 

no 1 1 co" j ^mrwr N , l i nmav 

IF(JaGT 0 I c TX(AL 1N£( T I 1 1) I I=*T 1*1 

T P t . GT . T r T X ( Al INF(Ill) I * t— 1 

— pTTTo.m ~r=M pt ? , 

1PUI.GT.NPT8) ! T-l 

ITDT-o-1 

rPTTTTT^Tt^PT'sr'''! TT^TfrANPTo; 

LOCFI I *TPlX(SLOP':(l I Il'jo-Ff IT 1 1 o-l 1 
HKFLL*! P! X (SLGPv U I U I«-8( U 1 1 





116 

117 


CC Ft. L = S1.0 t>F < IT + I I) 

RC ELL-SLOP?. 1 1 H ) *J + B< TTT J +1 
^n.nr r lL«,LTolJ LOCFLl-1 


1F(HJCELL.»0T.M/'XC ) H I CEL L = MAXC 
IF ( LOCEL L . GT .HICrLL >. CD TO 1105 

pn 1 1 04 L = l.nC£LL»Mir.FL L 

XOATAfl,!. , J»s1 
N0=NO+1 
pn 140 K = 1 ,4 


PK=K«-1 

NA=XDATA(KK»L,JI 
AN ( K ) = FLOAT IN A) 
fOMTINLF 
R ( 1 1 = AM ( 1 ) 



R ( 3) = ANi 3) 
R ( 4 1 = AN ( 4 ) 


R < 6) = 2=*AN(4)-ANtl} 
R ( 7 1 = AN ( 3) -AN ( 2 I 
R(8)=2.*ANt41-ANI2) 


R{9)=AN< 1 ) / AN I 2 1 

Piio)=awm/ANi3) 

P( 11I=AK|(1I/A.M(4) 


R(12)=AN(2l/AN(3) 
R< 13)- AN <21/ AN (4) 
on 2 50 K = 1 , NTF RM 


t>A = R (K) 

NA=IMT(R A) 

IF t R A „|_T . AM IM( K I > AM IN ( K I = P A 


I F ( S A o GT « A M AX ( K 1 ) AMAX(K1=RA 
IFf K„l T.9) GO TO 242 

N A = IMT C 1 0 . * R A I 


1F(NA»L?oOJ NA = 1 
IF(MA.GT.l28) MA = 128 

T FOFQIMft »!Q = IFRFQ(NA t KH-l 


SU M ( K ) = SUM ( K I 4-R A 
DO 250 KA=!,K 

ppnn(KA,X1 = PRQOKA»KUR (K I*R (K AI ' 

CONTINUE 
CONTI Ml! F 

CONTINUE 

CONTINUE 

CONTINUE 

IF(NO.LE.l) GO Tp 992 


TF (NCNoEO.l ) GO TO 100 
PO=FLOAT(NO) 

00 320 T=1,NTFRM * 

RA»'E< n = I I/Pp 

R STDI I ) = SORT ( ( PR 00 (I *1 ) -RO*RAVF ( I )*R AVF ( T))/(R0-1.}) 

nn 321 .1 = 1,1 

F C OR ( J , I >=0. 

I F ( R STD (1 ) ,F0=0,) on TO 321 

tF(RRTQ(J) qFQ.O.) GO tq 321 

PCO° { J„ II =! °ROp( J, 1 ) /PO-PAVE < J I *P AVEII II / l RSTD (II *RSTO< J »» 

CnNTTNUF 

CONTINUE 

FORMATCIMI , “IMMRFR of o^TNT<; = ',15«MN HOP '«A4 t 'WITH JN^ET O c 
1 « » c 3.1, 'UNITS, TAP L' *' , T4, , I 5,2A4, • INITIAL LINE = ',I4,» CEL 
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WRIT r (6, 10?) N0,FTFL0,0,PATF ,TAPF»ini» I02,L INET ,NCELLT 
PR T NT 103 

FORMAT! 4 0 4 ,10X, ‘ M $ R 4 1 ,4X , »»SSA ■ ,4X, 1 MSS7 4 , 3X, 4 MSS6- 4- 

IS $7-4 MSS6-5~ 2MSS7-5 MSS4/5 M SS4/6 M5S4/7 MSS5/6 MSS5/7 4 ) 
PRINT 104* IRAVn T > , I «1, 1 3 > 

FORMAT ( ' C 4 i:; 4 MEANS' , 3 X * 1 6 p _8_tU 

PRINT 105f(AMTN'( 1)* 1-1*13) 

PHrmATI •0« , 4 NT mi mu^S • .16F8.3) 

PRINT 106, ( AMA XI I) , 1 =1 , 13) 

FORMAT ( *0** * MAX T M'JNS ' ,16F8.3) 

PRINT 107* ( RSTPf I ) ,1 =1*13) 

FORMAT ( 4 O 4 * 4 5. TD DEV *,16F8.3) 

PRINT 116 

FORMAT! 4 - ' , 4 C n R R FLAT I ON CO"FF TCI FNTS 4 ) 

RRINT 111, (RC np <01 * I ),T = 01*MT r .PM) 


FORMAT! 4 0* , ' MSS4* ,4X,16F8.5) 

PRINT U2, (PCnc ( 02, I) , I = 02*NTRRM) 

FORMAT! 9 0 * . 8 MSS5 ' , 12X.15F8.5) 

print 113, (RC OR (03,1 ) , f =03 * NTF RM) 

FORMAT! '0 s , 4 M$$6 4 »20X,14Ffl.5) 

PRINT 114, (RCQR (04,1), I = 04,»TFRM) _ 


FORMAT! "O', •MSS7S28X,13F8»5> 

PRINT 115, (RCD»l<05, I ),T = 05*MTERM) 

FORMAT ( 4 0* , 4 MSS 6- 4 4 *34X,12rR.51 

PRINT 121, (RCOR (06, I ), 1-06, NTF RM) 

FORMAT! 4 0 4 , 4 2MSS 7-4* ,41X* 11F8.5 ) 

PRINT 12?, (RCOR! 07 ,T ) ,1-07, NTF RM) 


FORMAT! 4 0 4 , ’ MSS 6- 5 4 ,50X,10F8.5) 

PRINT 123, (RC0IU08* I) * 1=08, NTFRM ) 

FORMAT! 4 O 4 , » 2MSS7-5 4 , 57X,QqF8.5) _ 


PRINT 124, (RCOP 109* I )♦ I=09»NTERM) 
FORMAT! 4 0* , 4 MSS4/5 4 , 66X»8F8«5) 
PRINT J2M, ( RCQR 1 10, I ), I =1 0, NTF R M ) 


AT 

T 

AT 

( p 0 1 , **SS4/6’ ,74X f 
126, ( PC PR C 11 f n ,1 
{ ’OS c MS $ 4 / 7 9 ? B2 X t 

7PB* 

^11 

6F-B, 

T 

127, (PC0IU12, Dtl 

= 12- 

4T 

( «0 4 , 1 MSS5/6 4 ,90X, 

5F 8 , 

T 

128, !RC0fcU3,T),I 

- 1 3 


"WRITE (7, 131) £ I C LD,DAT-, ( R S T 0 1 1 1 ,1=1,13) 

WRITE (7, 132 ) FI ELD, DATE, ( AMIN (11*1 = 1,131 
WRITF (7, 134) FIR LD.DATP, ( AMAX! I ), 1=1,13) 

FORMAT! • 1 ' , A 3 , 1 4 , RF 5 • 2 , 5F4. 2 ) 
rnQMflT ( 4 2 4 , A?, 1 4,8F5.2,5P4.2) 

FORMAT! 4 3" , A3,I4,4F5<.2,4F6«2*5 c 4i>2) 


FORMAT! 4 4 1 , A3, 14, 8F5 »2, 5f4.2) 

CONTTNUF 

PRINT 1 08,FI r I.O,nATE ,TAP<^ , I 01 » I P? » L I NFT , Nf. c I.L T 
~ FORMAT! 1 HI, 4 HISTOGR A M 4 ,3X, »F IFI.O = 4 , A4, 4 TAPS, = 4 , 14,'-' ,15,2/ 
l 4 INITIAL LINE = 4 , 1 4, 4 CELL = 4 ,T4) 

L MI N= MI N 1 l am IV ( 1 ) * A M INI 2 ) , AM IN ( 3 ) , / M I N( 4T» A H1AU5) , AMIN! 6) , 
1 A M fN !' 7 ) , flM I N { 9 ) , /I Ml" N"<9) =10. , AMI N ( 1O)*10.*AM?NU1)*IO. , 

PAM IN! 12)*10.,AM1M(13)*1C.1 

l MAX = MAX1! AMAX! 1 ) ,A,M AX! 2) ,AMAX( 3) , AH ftX ( 4 ) ( 5) ,AMftX(6) * 

| m A X ! 7 J ft X ( 8 ) 'A M A X <9)^10, , A.v»x( 10 ) * lO. , A/*AX ( 11 )*10. , 

2AAIAXJ 12) *10., A M AX! 13)* 10. ) 

t f ( LM I N . I F«0) Lf^lN=1 


227 
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228 TMLMAX.GT.128) LMAX=128 

229 HO 400 I*LMIN,l.MAX 

230 400 PPINT 1 09.1 , < I fr ^-0>i , J) , J = 1 .NTF.BMJ 

231 100 FOPMATt !H t I5,16I 3X,T 5 ) ) 

232 WRTTE(6»1001 I H NM IN , FI EL D, OAT E, TAPE* 101, 102 

233 lorq F^RMATM 1« . « MINT MtJM LINE NO = » ♦ 1 8. 3X , j F 1ELP t A 4 t * T APP = S 

1I4, , - , »T5*2A4) 

234 HO 1110 J=1 » MA XL 

234 WP T TF < 6« ICO?) ( XPATft ( 1 . 1 , J ) , I»l, MAXC ) 

236 1002 FORMAT (• *,6411) 

237 1110 CONTINUE 

233 NNTMO^O 

239 GO TO 100 

240 991 WRI TF ( 6, 120) NTft PEA , NOATt 

241 I?n FfjPMaTI *0* . » F I FI 0 CARD PFQL’FSTS TAPE , .«L4,» RUT TAPE USED IS *_IAL 

242 GO TO 100 

243 990 WR I T F ( 6 » 503 ) NOA TE , N T A^F *NI 01 ,N 102, DATE , T APF , ID1 * I D2 

244 903 FORMA T( » 1 « . j, T APF RFO'JFS T E Q1 * . 1 4 i * - 1 . 1 5 . 2 A 4 ± ' ) — OTONOT MATCH TAPE. M Q 

lUNTFDl • »I4,*-' , I5,2A4,*) ’ ) 

245 GO TO 999 

246 99? I M NC N. FQ.l I GO TO 100 

247 WRTTF(6»119) F1FLD 

248 119 FOPMATJ 1H1, • FIFU 0 ' , A4, • TOO SMALL • ) 

249 GO TO 100 

250 994 TFINCELLT.LT. NCELL1M 00 TO 10 

251 995 WRTTF(6,il01 L INF( 1 V,L INET,NCB, NCELLT 

252 nr FORMAT! 1H t *FRRQR TAPF TOP FARtLTNE - * , T 5 , * T APEL * » 1 5 t * CELL g *t 

ITSc’TAPF CELL = SIS) 

253 1 F(NWIND.E0.1 I GO TO 999 

254 PFWTND NREAO 

255 MWIN0=1 

256 GO TO 10 

25 7 999 RETURN 

ENO 


250 
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t JDR 

i 01 ^ C MS I ON LINFIOIO) ,NC Fit. 1010), SLOPE I 010},THFTA( 010) , All NEI010) 

_Z DIMENSION 9 t 01 O> .Ar-LL ('-’101 ' 

3 T M T FGF<? rt 4 DATE, TAPE 

4 I NT FGF9* 4 SCG,PX 

F TNT g G c P*2 XDATAI 5,64,64) _ 

ft INTEf,F9*4 HICFL1 

7 0! MENS TON CELL (7 ) 

_? DIMENS I ON IFREOf 128, 16) . 

0 DIMENSION AMIN ( 1 6 ) » A^AXI 16 ) , $UM( 16) ,PR3Dt 16,161 ,9 I 16) 

|C DIMENSION PAVE ll6)»PSTDfl61,RCQP(16,16),ANI4) 

JJ PEAL FIELD/’ ALL’/ 

12 NPEAn=9 

13 NTCRM=13 

14 0 = 1.0 

15 NCN=0 

!ft NWINO-O 

JJ NSFG^O : 

18 maxC=64 

1ft MAXL=64 

20 READ 444,D,NSN, NSW! ; 

21 444 FORM A Tl c 5« 2 , T 1»T1) 

22 5 READ! 5, 500, FND-999) NDATE , NTAPF ,NIO 1 , N I 02, NCELl 5 

23 800 FPRMATII4,IX,I5,2A4,2X,I5) 

24 1.00 CnNTIN'JF 

25 TF(NSWI.PO.l) GO TO 101 

2ft DFiDf 5,1C-P0,FH0=999) SFG,NPTS, (LINFI1), NCELL IT), I=1,4),NFTLE,NC, 

1 N T AP FA , F I ELD 

2 T 1000 FORMAT! II, T1*4CI4*I3),1X, 1 1 , 1 1* 2X, I 4, 3X , A3) 

_29 PP IMT 11Q2,SEG,NPTS, 1 1 INF t Li » NC C LL ( 11 , 1=1 .4) * NF ILF, NC ♦ NT APE A , F I FLO 

29 1102 FftRMATClH , ’ 5 FG = • , 1 2 , * PTS = • , I 2, 4( I 5 * IX, T 4, 1 X) * * F I L F= ' « I 2 ♦ 'NC = ' , 

1 ! 2, *T APE*', 1 5, 'FIELD*', *41 

JO GO TO 202 

31 101 CONTI NOE 

32 PH AD I 5* 1 003 * END= 999) SEG, NPTS , ( L INF ( T ), NCELL ( ! ),I = 1,8),NFILE,NC, 

1NTAPFA, FIELD 

33 1003 FORMAT (J 1,1 1,S{ I4»T3),2X» II, 1 1, IX, I 4, IX, A3) 

34 PPTNT 1 1 oi , SFG, N ft TF » II. INF (I) ,NC FLU I) , I = 1 , 8 ) , NF IL E , NC,NT APFA, FIFLD 

35 IIP 1 FORMAT 1 1 H , »$FG*», T2,«PTS*',I2,8(15,1X,I4, IX) ,'FILE*» ,I2,»A'C=», 

1T2,«TAPE=',I5, , FIEL0=',A4) 

36 202 CONTINUE 

37 TFINTAP FA »NE. NDATE) GO TO 991 

38 I F ( NCN. c Q. 1 ) GO TO 10 

39 N 0=0 

40 no 220 - ! 

41 DO 190 J=1 , 128 

42 190 I FPFOI J , I ) = 0 

4 ? AMIN! I 1 -128 * 

44 ANAXIT)=0. 

45 8"N(T)=0o 

46 00 21° J^loNTOP'A 

47 ppnnu,i)*o. 

4R 210 CONTINMF 

4" ? 2 n CONTINUE 

50 IFIS^G.FO.NFEG) GO fn R70 

51 ]0 RPAniNPFADI OAT." ,T/n>t , jpi, TD2,L !NFT t NCFI.LT 

52 TF|OAT'*.N,* „Mil AT”) r-O TO OQO 

"53 'fp | T A PC .N c .NT A t>c f rio' Tn qqo 

54 IF( IDKNF.MIOl) GO TO 990 

55 T c I I Q2.NE.N 10? ) GO TO 9QQ 




Appendix V 


Computer Program and Flow Chart to Create 
Contour Plots on Calcomp Plotter 



CONTUR 


CREATE CONTOUR PLOT ON CALCOMP PLOTTER 

Programmer: Jay Alloway, Kansas State University Computing 

Center, May 1073. 

Language: FORTRAN IV for 3G0/370 with calls to CALCOMP plot 
subprogram PLOT. 

Purpose: This subroutine prepares a contour plot from a 

rectangular grid of uniformally spaced values. 

Calling sequence: 

CALL CONTUR (GRID, ROW, COL, CU, NV, XLEN, YLEN) 
where 

GRID contains the values (REAL*4) to be plotted. 

These could be reflection intensities, levels 
of radiation, etc at equally spaced points. 

ROW is the number of rows to be plotted down the 
paper and the size of the first dimension of 
GRID. ROW is integer. 

COL is the number of columns to be plotted across 
the paper and is the second dimension of GRID. 

COL is integer. 

CU Is a vector (REAL*4) containing the 

desired contours. 

NV is the number of contours in the CU vector. 

XLEN is the length ( REAL*4 ) of the x-a;tis (down 

the paper) in inches. 

YLEN is the length (REAL*4 ) of the y-axls (across 
the paper) in inches. 

NOTE: CONTUR assumes the pen is located at location 

(0.0, 0.0) when called and does not reorigin before 
returning. Location (0.0, 0.0) is GRID(1,1). 

Location (XLEN, YLEN ) is GRI D (ROW, COL ) . 

Method: CONTUR looks at the grid point by point forming a 

square at each point consisting of the points to the 
right, below, and right and below the point is question. 
Except for two coses where two lines nay be drawn, only 
one contour line may appear in a square. CONTUR assumes 
n uniform media within a square. The line drawing 
technique outlines flat (equal contour) areas. There 
Is a total of ( ROW- 1) * (COL-1) squares In GRID. 
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SUBROUTINE CDNTUR < X , N , M f CU , NC , XS T ZE , YS I ZE ) 


C 

C 

C 

C 

c 

c 

c 


c 

c 

c 

c" 

c 

c 

“c 

c 

c 


SUBROUTINE CONTIJR CREATES 
FROM AN ARRAY OF POINTS. 


A CONTOUR PLOT ON THE CALCUMP PLOTTER 


WRITTEN BY JAY 
IN MAY 1^73. 


ALLOW A Y f KANSAS STATE UNIVERSITY COMPUTING CENTER 


DIMENSION X ( N t M ) , CU(NC) 

" LOG I CAL *1 ISW 

INTEGER UP, DOWN 

DATA T 0 L r UP, DOWN /IE-7,3,2/ 


CALL GUNTUR (GRID, 
WHERE GRID 

' "ROW 

COL 


c 

c 

c 


cu 

NV 

-*LW' 

ylen 


ROW, COLt CU, NV, XUENt YLEN) 

CONTAINS THE VALUES TO BE PLOTTED. _ 

I S A N I NT EGER G I V T NG” ThE* NU M R Eft OF ROWS TO 
PLOTTED AND THE 1ST DIMENSION OF GRID. 

IS AN INTEGER GjVlNGjHB NUMBE R OF COLUMNS 
PLOTTED AND THE 2ND DIMENSION GO. GRID. 
CONTAINS VALUES OF THE DESIRED CONTOURS • 

IS THE NUMBER OF CONTOURS IN CU. 

T S T H FT^TTGTTTTIF - 
is the length of 


BF 


TO RE 


TFTF 

THE 


Y-AXIS 


(COLUMNS) , 


CNTROOIO 
CNTR0020 
C NT R003 0 
C NT R 0040 
CN TR0050 
CNTROOIO 
CN T R 007 0 
(TNTROORO 
CNTK0090 
CNT RO 100' 
CNTROllO 
__CNTR01?0^ 
CNTR0T30 
CN T RO 1 40 
CNT RO ISO 
~CFTR0160 
CNTR0170 
CNTR01R0 


CNT RO 1 9 0 
CNTR0200 
CNTR02 l 0 
cntroftf 

CNTR0230 

CNTR0240 


~F~ 



MHQli&SslSM 

C 

DEFINE BASIC FUNCTIONS NEEDED 


CNTR0260 

c 



CNTR0270 


“TCLOC I A n = ( A I /'AMr*TST7E 


CNTR02B0 


YLOC(AJ) = ( AJ/AM )*YSIZE 


CNTR029 0 


A INC ( AS ?AE ) = ( C-AS > ! ( AE-AS ) 


CNTR0300 


e’OU'JTUTPT) = ABS(C-RT) " 


CNTR0310 


BETWEN ( AS * AE I = ( C-A S ) * ( C-AE ) 


CNTR0320 

c 



CNT RO 33 0 

C' 

‘ INITTAUTTE’ CONSTANT'S 


CNTR0340 

c 



C NT RO 35 0 


NX =* N-l 


CNTR0360 


MX = M-l 


CNTR0370 


AN - NX 


CNTR03R0 


AM st MX 


CNT RO 390 

u~ 



CNT ft 0400 

c 

LOUP THROUGH ROWS 


CNTft04l0 

c 



CNTR0420 


' DO 2 T5“TX"= 1TNX7Z 


CN 1 R043 0 


I = IX 


CNTR0440 


ISW = .FALSE. 


CNTR0450 


-T,cr-T0- T5 


LNIKU4 hU 


10 I = IX+l 


CNTR0470 


IF < I . GT .NX) GO TO 215 


CNTR04R0 

' 

ISW = .TRUE. 


CWTRU?? 1 ?!} 


tfl 

►— « 

T) 

11 

1— 1 
+ 


CNTR0500 


AIM! = f - L 


CNTR0510 


Ai = I 


CNTRD‘570 

c 



C NT RO 53 0 

c 

LOOP THROUGH COLUMNS 


CNTR0540 

■ cr 



CTO R~0"5"5 0 


DO 210 J X= 1 , M X 


CNTR0560 


J = JX 


CMTR0570 

— ^ 

“FOR~‘FASTER PLOTTING. P ITT! EVER Y - OTHER ROW BACKWARDS 


CNTR05"R"0 


IF (ISW) J = M-JX 


CNTR0590 


JP1 = J+l 


CN TR Of> 00 






i 


nno 
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AJM1 = J-l CNTROM 0 

AJ = J CNTR0020 

C GET VALUES OF CORNERS OF SQUARE CNTR0630 

A 1 = XI 1 ' J I “ ~CNTR0i540 

A2 = X ( I ,JP) ) CNTR0A50 

A3 = X< I P 1 ,J) C^TROAAO 

A4 = 'X(lPU.iPl) ’ CNTR0A70 

cntroaro 

LOOP THROUGH CONTOURS _ _ _ CNTROA9Q 

"" CNTR07 00 

DO 205 R=l,NC CNTR071 0 

C = CU(K) CNTR07P0 

C ■ PLOT BEST CURVE THROUGH THIS "SQUARE * CNTR0730 

IF (EOUALIAD.GT.TOU GO TO 70 CNTR07 AO 

IF (EQUAL! A2I.GT.T0l.) GO TO 35 CNTR0750 

IF ( EQUAL ( A3 ) ,GT .TOL j GO TO 20 ' ' CNTR0760 

IF (EQUAL! A4) .LE.TOU GO TO 205 CNTR0770 

, CALL PLOT (XLUC(AI)fYLOC(AJMl). UP) CNTR0780 

CALL PLOT ( XLQC ( A I M 1 ) » YLOC ( A J ) . OOWN ) CNTR0790 

GO TO 205 CNTROROO 

20 XX = XLOC(AIMl) CNTR08I0 

" "“Call plot "Txx , yloc I ajm'i ) , up ) cntrob 20 

YY = YLOC(AJ) CNTR0830 

IF (EQUAL! A4) .LE. TOL) GO TO 30 CNTR0840 

I F- ( rt£TWEN'( A3 t “A4 jTLE“.OTO ) GO TO 25 CNTR0850 

CALL PLOT (XX, YY, OOWN) CNTR0860 

GO TO 205 CNTR0870 

2 5 ~C A L L PLOT ( X L OC < A I ) , Y LOG ( A J Ml + A I NC ( A 3“,A4) ),DOWN> CNTRORRO" 

CALL PLOT ( XX , YY , DOWN ) CNTR0890 

GO TO 205 CNTR0900 

30 CALL PLO T”i XLOC ( A f),~Y Y , OOWN j CnTROOUT 

GO TO 205 CNTR0920 

35 XX = XI.OC ( A I M D CNTR0930 

Y V = ~Y L OC !~A j M 1 T CNTR0940" 

IF ( EQUAL ( A4 ) . GT .TO!. ) GO TO 45 CNTR0950 

IF (EQUAL (A3) .GT .TUL) GO TO 40 CNTR0960 

CALL' PLO T ' ( XX , YY » UP’)' CATfR09Tcr 

CALL PLOT ( XLOC ( A I ) . YLOC < A J ) , DOWN ) CNTR0980 

GO TO 205 CNTR0990 

40" I F ( BETWEN ( A? , A3 IVST r<570T C'fT TrTT05 CNTR1000 

CALL PLOT ( X X , Y Y , U P ) CNTR1010 

D]FF = A I NC ( A2 » A 3 ) CNTR1020 

~C A L C P Co T (X c OC ( A I Ml + D ITFT, Y L 0“c ( S J -0 I F F )" ,U OWN) CN T RIOTO - 

CALL PLOT ( XLOC t AD ,YLOC(AJ) , DOWN) CNTR1040 

GO TO 205 CNTR1050 

4 5 " I F ( E QU A I. ( A 3 ) . L F . T 0 L FT, 0" T 0 "7)71 — CWTRKTW 

!F (RETWENI A2,A4) .GT.O.O) GO TO 50 CNTRI070 

IF ( BET WEN ( A3 , A4 ) .GT . 0 . 0 ) GO TO 55 CNTR1080 

“CALL PLOT CXLCIC'C A I I ,YI_OC ( AJM1VA I NCCA37F41T7UP'I WTITTOW 

CALL PLOT ( XLOC( AIM1 ) ,YLF1C( AJM1 ) ,OOWN) CNTRllOO 

CALL PLOT ( XLOC ( A I Ml + A INC ( A2 t A4 ) ) , YLOC ( A J ) , DOWN ) CNTRU10 

go to 205 c wrmr ?ir 

50 IF ( BETW6N( A?, A3 ) .GT.0.0 ) GO TO 205 CNTR1I30 

IF (BETWEN(A3,A4). GT.0.0) GO TO 205 CNTRU40 

CALL " PLOT ( XX » YY »UP I CNTfrYlSO" 

DIFF = A I NC ( A2 , A 3 ) CNTR1I60 

CALL PLOT ( XLOC ( A IM1+DIFF ) , YLOC ( AJ-U IFF) .DOWN ) CNTR1170 

CALL" PLOT ( XLOC ( 41 ) , YLOC"! A JM 1+AT NC (A3, A4TY, "DOWNY CTTTRTY ATT 

GO TO 205 CNTRL190 

55 IF (_HETWEN(A2,A3) .GT.O.O)_ GO TO 208 CNTR1200 






V 
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CALL PLOT I XX » YY, UP ) 

□IFF = AINC<A2,A3) 

CALL PLOT ( XLOC( AIM1+01FF) , YLOC < AJ-DlFF ) .DOWN > 


CALL PLOT ( XLOC ( A I M l + A INC H A2 , A4 > ) ,YLI)C I A.i f.OOWNT 
GO TO 705 

60 CALL PLOT |XX,YY,UP) 

IF (BETWENI A2,A4) .LE. 0.0) GO"TO 65 

CALL PLOT IXL0CIA1) .YY.OOWN) 

GO TO 705 

65 CALL PLOT ( XLOC I A IM1 +A 1 NC ( A2 , A4 ) ) , Y LOC t A J ) . DOWN ) 
CALL PLOT ( XLOCI A f ) ,YY,OOWN) 

GO TO 205 

70 IF I EQUA L ( A2 ) . GT . TOI. ) GO" TO '1 05 
IF < EQUAL < A4 ) .GT .TOL ) GO TO «5 


CNTR1210 
CNTR1270 
CNTRl230_ 
CNTR 1 7 AO 
CNTR] 260 
CNTR 1760 
CNTR 1270" 
CNTR 17R0 
CNTR1290 
"CNTR 1300 
CNTR1 310 
CNTR 1320 
‘CNTR 1330* 
CNTR 1340 
CNTR1360 


' " CALL PLOT IXL0UA1 ) fVLOCI AJMl) ,UP) 

CALL PLOT 1 XLOCI AIM1 ) ,YU1C ( AJ ), DOWN ) 
GO TO 205 


CTmT3"60 _ “ 

CNTR1370 

CNTR1380 

'75' CALL PLOT I Xl.OC < A I ) , YLOC I A J ) f UP ) 

IF (BETwENfAl ,A3 ) .GT ,0,0) GO TO RO 

CALL PLOT { XLOCI AIMl + AINCI A1 , A3 )), YLOCI AJM1 ) ,DOWN ) 


CNTR 1390 
CNTR 1400 
CNTR 1410 

BO C ALL °LOT ( XLOCI AIM 1 ) , YLOCf AJ) ,DOWNj 

GO TO 205 

B5 IF ( EQUAL ( A3> .LF.TOD GO TO 100 


CTTTR 14?0 
CNTR 1430 
CNTR 1440 

!F""( BfTWENt A1 , A3 ) . GT . 0. 0 ) GO TO 90 
IF (BETWENI A3, A4) .GT.O.O) GO TO 95 
CALL PLOT ( XLCIC ( A I Ml +A I NC ( A 1 , A3 ) ) , YLOC ( A JM l ) , UP ) 


CNTR 1450 
CNTR 1 4 60 
CNTR 1470 

CALL PLOT { xLOCTAl M 1 ) , YLOC ( A jT ,"DOWN T 

CALL PLOT ( XLOC< AT ) , YLOC 1 A JM 1 + A INC ( A3 , A4 ) ) , DOWN ) 
GO TO 205 


CNTriTRO 
CNTR 1490 
CNTR 1500 

90 IF I BETWENI AT, A4) .GT.O.O ) GO TO 205 

IF (BETWENI A3, A4) .GT.O.O) GO TO 205 
CALL PLOT t Xl.OC ( A I Ml ) , YLOC ( A J ) ,UP ) 


CWTR 15T0 
CNTR 1520 
CNTR1530 

• "0 I FF" = ’ A I NC ( A 1 , A4’ ) 

CALL PLOT ( XLOCI A IMI+DI FF ) , YLOCI AJMI+OIFF ) ,DOWN) 
CALL PLOT (XLOCIAI ), YLOCI AJMl+AINCI A3, A4) ) .DOWN) 


CNTR l 5^0 
CNTR1550 
CNTR1560 

- 6IJ T0 205 

95 IF I RETWENI A1 ,A4) .GT.O.O) GO TO 205 

CALL PLOT I XLOCI AIMl+AINCI Al, A3) ), YLOCI AJMl >, UP) 


CfTTRT5TC5 
CNTR 1 6 RO 
CNTR 1590 

- "DIFF -="AINC(Al,A4) 

CALL PLOT ( XLOC ( A I Ml + I) ] FF ) , YI.OC ( A.IM 1 +0 I FF ) , DOWN ) 
CALL PLOT (XLOC(AIMl), YLOCI AJ), DOWN) 


Cn ) k 1600 
CNTR1610 
CNTR 1 6 ?0 

■ ~~ "GO TO 205 ' 

100 IF (BETWEN(A1,A4) .GT.O.O) GO TO 205 
CALL PLOT I XLOCI A I) , YLOCI AJM1 ),UP) 


CNTR 1640 
CNTR 1650 

— ~OIFF~=' AINC t A1 , A4) 

CALL PLOT ( XIJ'JC (A IMI+DI FF) , YLOCI AJMI + OIFF ), DOWN ) 


lN i k 1 b 6 u 
CNTR 1 670 - 

CALL PLOT (XL0CIAIM1) , YLOC ( A J ) , DOWN ) 


CNTR 16 80 

GO"TO" 205 


CATTRTTT'TO 

105 IF ( EQUAL l A3) .GT .TOL ) GO TO 130 
IF < EQUAL (A4).LE.T0L> GO TO 120 


CNTR1700 

CNTR1710 

“ “I F TBETWEN|A1,A2) .GT .0,0 ) GO TO 110 


CWTRT7T0 

• IF I BETWENI A2 , A4 ) .GT . 0. 0 ) GO TO 115 


CNTR 1730 

CALL PLOT ( XLOCI A1M1) , YLOCI AJM1+A1NCIA1 ,A2 )) , UP) 


CNTR 17 40 

- CALL PI-flT I XLOCI AI) , YI.OC TAJMl ), DOWN ) 


CNTR 1T5TJ 

CALL PLOT <XL0CIAIM1+AINC(A2*A4) ), YLOCI AJ) , DOWN) 


CNTR1760 

GO TO ?05 


CNTR 1 770 

— — 1 IT) "I F (BETWENI Al , A4 ) ,GT .0.0) "GO TO 205 

IF I BETWENI A2 ,A4 ) . GT.O.O) GO TO 705 


CNTR 1TTT0 
CNTR1790 

CALL PLOT I XLOC I A I ) » YLOC I A JM 1 ) , OP ) 


CNTR1RO0 
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OlFF = AINC(Al,A4) 

CALL PLOT <XLOC(AIMl+D!FF ) ,YLOC( AJMl+DIFF) .DOWN) 

CALL PLOT ( XI.OC ( A I M1 + A INC { A2 , A4 ) ) , YLOC ( AJ ) .DOWN I 

DO TO 205 

IF ( RETWENt Al ,A4) .GT.O.O) GO TO 205 

CALL PLOT (XLOCI AI > ,VLOC( AJMU.OP) 

D1FF = A I NC (Al. A4 ) 

CALL PLOT (XLOCI AIM1+OIFF) ,YLOC( AJM1+OIFF) .DOWN) 
CALL PLOT ( XLOC< A I Ml ) , YI.OC < A JM 1 + A I NC (. A 1 , A2 ) ) , DOWN )_ 
GO T 0 205 
XX = XLOCIAI ) 

CALL PLOT ( XX » YLOC ) A JM 1 ).. UP ) 

IF ( BETWEN { Al . A2 ) .GT.O.O ) GO TO" 1 25 

CALL PLOT (XLOC(ATMl)tYLOC( A JM 1 + A I NC ( A 1 . A2 ) ),DOWN) 

CALL PLOT ( XXtYLOC(AJ) .DOWN ) 

GO TO '205 

IF { EQUAL (A4).GT.T0L) GO TO 145 

IF (RETWENI A1,A2) .GT.O.O) GO TO 1 35 

"IF ( BETWENI A 1 , A3 ) .GT^O.O) GO TO 140 
CALL PLOT {XL0C(AIM1+AINC(A1,A3) ) tYLOC(AJMl) tUP) 

CALL PLOT ( XL0C(AI},YL0C(AJ)j,DOWN) 

"TALL’ PLOT ( X L OCTA I MlTt YLOC ( AJM1+A INC ( Al , A2 ) ) tOOWN ) 
GO TO 205 

IF (BETWEN! A2 , A3 KGT. 0.0) GO T O 205 

I F (BETWE'N ( AT, A3 ) tGT'. 0'.0rG0 TO 205 

CALL PLOT (XL0C(A!Ml+AINC(AltA3) ) tYLOC(AJMl).UP) 

_01FF_= A I NC I A2 . A3 ) _ _ 

CALL PLOT ( XL Of. ( A IM 1+D IFF ) , YLOC ( A J-DI FF ) .DOWN ) 

CALL PLOT (XLOCtAl ) , YLOC ( A J ) t DOWN ) 

GO TO 205 

“IFTBETWENI A2.A3) .GT.O.O') GO TO 205 
CALL PLOT (XL0C(AIM1).YL0C(AJM1+AINC(A1,A2))»UP) 

_0l£F_= AINC<A2,A3) 

CALL PLOT ' ( X LOG ( A'j M 1 +D T F F ) .YLOC < AJ-OIFF ) ,DOWN) 

CALL PLOT (XLOCtAl ) . YLOC ( A J ) , DOWN ) 

GO TO 205 J __ 

I F'Tb'E TW EN ( A 1 , A 3 ) . Lt . 0 . 0 ) GO TO 170 
IF (BETWENU1 ,A2 ) .LE.O.O) GO TO 150 

IF (BETWENt A3,A4) .GT.O.OJ GO TO 2 05 

CALL PLOT ( XL OC (A I T t YLOC ("A JM 1 + A I NC ( A3 ♦ A4 ) ) , UP ) 

0 IFF = AINCt Al ,A4) 

CALL PLOT ( XLOC (AIM1+DIFF ) .YLOC < AJMl+DIFF ), DOWN ) 

— CALr'PLOTn X L 0 C'< AT M 1 +'Af NCTA 2»A 4 ) ) , YLOC ( AJ ) , DOWN ) 

GO TO 205 

i IF (RETWENt A2,A4) .GT.O.O) GO TO 155 

'••CALL PCOT'I XLOC (£ I M I } VYUtTCT A J M' ITaTFIC I Al , KT ) ) .UP ) 

0 1 FF = A INC ( A2 * A3 ) 

CALL PLOT ( XLOC (AIM 1 + D IFF). YLOC ( A J-P IFF) .DOWN ) 
CALL" PLOT ( XLOC ( AIMl+A INC ( ‘A 2 , 'A 4] rrVCnCTTTJ 1TD 0 WM1 
GO TO 205 

IF ( ( ABS( A4-AI ) .GT.AftS! A3-A2 ))) GO TO 150 
"IF “I BETWEN ( A2 ,A3 Y TGT r0;0T"'G"0“Tfr“205 
DIFF = A 1 NC ( A2 t A3 ) 

XX = XLOC ( AIM1+0JFF ) 

" YY =->LOC(AJ-l)lFFf 
GO TO 155 

I IF (RETWENt Al ,A4). GT.O.O) GO T0_?05 

DTFF'"= A INC ( A it A4 ) 

XX = XLOC ( A I M l > 01 F F ) 

YY = YLOC( AJMl+DIFF ) 
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r ( XLOCUI ) ,YI.OC(AJMl + AINCI A 3 ,A 4 ) J ,UP ) 

T I X X f Y Y * DOWN ) 

T f ^ r_^i ) » YI.OC ( AJMi+A f NC ( &1 1 A2 ) ) * DOWN ) 

5 N< Al,A 2 ).GT.O.O) GO TO 185 
EN | AX ,A 4 ) .GT. 0 . 0 > GO TO 175 

r ( XLOC! AIMl + AINCI A l ,A 3 ) ) , YI.OC < A JY 1 T,'UP) 

I NC ( A 1 * A4 ) 

r < XI.OC ( A IM 1 + DIFF ) , YI.OC I A.JM 1 +OIFF ) , DOWN ) 
r ( XLOC ( AIM 1 V , YUOCI AJM 1 +AINCIA 1 «A 2 )T,OOWNT 
) 

r ( XLOC ( A I M 1 ) t YLOC ( A JM 1 + A INC ( A 1 1 A 2 ) ),UP) 

D.GT.TOi.) GO TO 180 — 

(XL 0 C(AIMl + AINC(A 2 ,A 4 ) | ,YLOC<AJ) ,OOWN) 

: ( XI. OC ( A I )_, YI.OC ( AJM1 + AINC ( A3 t A4 ) I, UP > 

[ (XLOC (AIM 1 + A INC ( Al ,A 3 ") I ,Y L 0 C (A JM'TT , D OW'NT 
) 

I X L OC I A I M 1 ♦ A I N C ( A 1 » A 3 ) ) »YL 0 C(AJM 1 ) t DOWN) 

( XLOC ( A I ), YLfJC(AJMl+AINC( A 3 ",A 4 ) ),UP ) 

( XLOC (AIMl+AlNC(A2tA4) ) ,YLOC < A J ) , DOWN } 


CALL PLOT { XLOC ( AIM1+A1NCI Al , A3) ) » Yl 0CIAJM1 ) .Up ) 
_OIFF = AINC ( A2 t A3 ) 

CTixnprrn rxxuc ( a imt+ui fftt^uici atj-uiff j t unwNi — 

CALL PLOT {XLOCIAI ) T YLOC( AJMl + AINCI A3,A4 ) ) , DOWN) 
GO TO 205 

IF < ( ABS ( A4-A U .GT. ABS ( A3-A2 > > ) GO TO 195 
IF (BETWENIA2 »A3) „GT .0.0) GO TO 205 

“DTF F =” A'l WCTT2TA3 - ] 

XX = XLOC (AIM1+0IFF) 

YY = YI.OC ( A J-0 I F F ) 

"G O T O “200 

IF (BETWENI Al , A4 ) .GT.O.O) GO TO 205 
DIFF = A I NC ( A 1 * A4 ) 

"TTST = XL'OCTA IMl + DTFTT : — 

YY = YLOC( AJM1+DIFF ) 

CALL PLOT ( XLOCI AIM1 + A INC( Al , A3 ) ) , YI.OC I AJM1 ) -UP) 

"X'ALC“PTm~rXX7^Y7DnWN5 

CALL PLOT ( XLOC ( A I Ml+A INC < A2 * A4 ) ) » YI.OC ( AJ ) ,00WN ) 

~CWO-OF-|;tiOP’S 

CONTINUE 

"CONTINUE - — 

IF I.NOT.ISW) GO TO 10 
CONTINUE 

“RETURN — 

END 
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Algorithm to Enhance Variations 
Within a Category 
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Let u 4 , S. be the mean and covariance of category ii i ~ 1». . . ,k. Define 
>11. 

^ik “ ( li i ” ^ ^ u i " \)’* Wc "'tmld Ii3<e to bring to zero all points from cate- 
gories other than Ci t , enhance all the variation in category c^, and show this 
variation about a point away from zero. 

The Algorithm 

-1 

Compute .t., i - 1, 2, . . . ,k, i 4 k where t. is the eigenvector of (S^ + B^) 
with largest eigenvalue. 

Normalize t. so that 

9 

l i ~ V “ 1 i=:1 *--*« lc » i5 ^ k * 

This normalization makes it so that 

i 

if x e c k , E |V [x - u. jj= - u.) = 1 

and if x £ c^, x e c ^ for some j, and E [t.' (x - Uj)] = 0 

For any x compute t. f {x - u.). Define f(?f) = min t/(x - u.).. 

i 

i=3 k 

Tiie following theorem and proof show how t. is determined. In this 

theorem the covariance of category Cj is used rather than the combined 

private Communication from R. M. Haralick, Remote Sensing Laboratory, 
University of Kansas, Center for Research, Inc. . - 
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covariance Sj c , within category Cj t , and B^, between categories Cj and c^.- 


Hence the theorem obtains ^ from S." 1 ^ rather than from S i A (Sj c + D^). delu- 
sion of the bctwccn-catcgory vax'iance term B^ c m alecs use of the separation 

of the category sample means. 

N V ft'xlcol ,, 

Theorem: Let f: R R be defined by f(t) = Then f has a 

global maximum of e where t is the eigenvector ofS-^Sg having largest eigen- 
value e. (Sj is the covariance matrix for category c^, i = 1,2). 

Proof: To maximize — ^ C ~] , we will maximize V [t’x | C2’ under the 

V [t l x | Cj] 

constraint V [t’x [ c-^] = k. 

Consider g(t) = V [t’x [ c 2 ] - e {V [t'x [ - k) . 

® E^Jlt'x - E [t’x|c 2 )]| 2 | c 2 - e^||t’x - E [t'x[c 1 ][| 2 |c 1 j-k 

Expanding the norm and letting u 2 = E [x ] c^] and - E [x | c^], we obtain 

g{t> = E [(x - u 2 )'t t’(x - u 2 ) | c 2 ] - e^E [(x - u^’t t’(x - | cj 


- 1 , 


a E [trace (tt 1 (x - u 2 ) (x - u 2 )' j Cg] 

- c^E [ trace (tt' (x - {x - u^’ -kj 


Since tx'acc is a linear operator, wo may interchange die expectation with the trace 



J 


C<t) 


= trace [tf (x - u 2 ) (x - x^V | c 2 lj 

- c /trace E [tt» (x - uj (x - u^’ | c^ - 


G- 


» trace (tt f E [x - \x^ ~ u 2^ ? | c 2^ 

- e^trace tt' E [(x -Uj) (x - Uj)' | c 1 ] - 

- trace (tt’S^ - c (trace tt’S^ k) 

» 

= t’S 2 \- e (t'S 1 t - k) 


Lett* 


t\,S 2 = s M,and Sl * s 



.( 2 ) 


Jj I)’ 




E (tx y-E *i E ?) + ek 

i = 1 3 = 1 V 7 

A necessary condition lor g to be maximized is for |^ = 0, n = 1 N. 

Expanding the summation for g so that we can talte partial derivatives easily, 
g(t r t H ,= gt 1 |t i ( s<f-es 1 X) ) 


i/n 


N 


+ *n E l j (s jn ~ ° S J^ 

j=l J 
• j^n 

+ Y^t. t (s< 2 > - c s^ 1 ) ) + c k 

+ Z-» i • n ' m ni 

i-1 

i/n 


Hence, 


N 


4Vgv 

j/n 


( 2 ) 

in 


c .« ) + ^ h ( s<2) 
i=l 
tin 


in 


m 


0 s'" ) 

Til ■ 


= 2 


N 


E <4> - • W s 

i=l 


+ 2t <sg) - e s£>) 
n ’ nn nn 


t by symmetry of and S^. 


Setting = we have 

N 

N 

S* te(? - e ) t. — 0, n - 1,2, .... N; 

jr[ m ni 1 


or in matrix form 


<S 2 ~ eS^ t= 0. 

Hence, S 2 t= eS^ or et so that it is necessary for t to be an eigen 

vector of S-^Sg having corresponding eigenvalue e. 

To see which eigenvector, we will evaluate f(t) for those cases when t is an 


eigenvector 




f(l) = and S t \t« ct. 

t'Sjt 1 2 
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f (t) = -. t<S 2 t 


VS ^l S 2 t 

0 


• t’S 0 t 


t ’ S l S l ' ls 2 t 


t'SA ‘ 

“ C 2_ - 0 

fs 2 t 


Therefore it should be that eigenvector of with largest eigenvalue e 
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Computer Programs to Implement the 
Algorithm of Appendix G 
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C 

C THIS PROGRAM IS AN IMPI.ENTAT10N OF AN ALGORITHM TO GNCHANC6 CATEGORY 

C VARIATIONS, PERIODIC REFERENCES ARE MADETO A PAPER E NTITLE D ENCHANCIN G 

C CATEGORY VARIATION FROM WHICH THE' ALGORITHM TTCiMES FROM AND WILL RE REFER ED 

C TO AS THE ECV PAPER. 

C 

REAL I AA , I AR , l AC » I AD * I A AS, I ABS , l ACS * J Al)S 

REAL I ARM, IRCM, ICDM, 1ACM, 1 RDM , I ADM 

INTEGER*? CAT EG, CL EN ,COL ,RDW ,MCOU_ ,MR 0WL,CAT ,HOWRED 

I N T EGER* 2 X ( 5 ,12ft, 12 ft) , Ft. AG I 100 ) 

DIMENSION B( 10) ,WORK V 1(4) ,WORKV2<4> 

DIMENSION RU I K ( 4 , 4 ) ■ 

‘DIMENSION FMT ( 15’)” 

DIMENSION S ( 5 , 4 ,4 ) 

DIMENSION SI (4,4) 

DIMENSION SK ( 4,4) * 

DIMENSION SKA (4,4) 

DIMENSION T I (4,4) — 

"DIMENSION' 0(4,2)" 

DIMENSION 0 A { 4 > 

DIMENSION UCAT EG ( 4 ) 

DIMENSION Ol (4) ,t)2(4 r ) 

C 

C THE FIRST DATA CARO CONTAINS THE NUMBER OF T HE CAT EGO RY OP INTEREST, 

c The" numreTR - of~c ategdr i es * The nuiMbFr OF" columns and rows in the input 

C MATRIX, AND A NUMBER FOR THE READ TYPE FOR THE INPUT MATRIX. 

READ (5,200) CATEG ,CLEN,MCOLL ,MROWL , NRE AD , HOWftEO 

WRI TE"T6",2rn CATEG', C'CEN ,MCOLL‘,MROWL /NRE AD , HOWRED 

c 

C READ THE FORMAT FOR READING IN THE INPUT MATRIX X. 

READ" 1 '572T0T"F W'T 
C 

C INITIALIZE THE INPUT MATRIX X. 

D 0 12 r~ K^lTl 28 ‘ 

DC 120 J= 1 ,128 

X(l,J,K)=l 

00 1 2 0' T= 2 V5 

X ( I , J , K ) =0 

120 CONTINUE 

T 2 1" CONT INtJE ' ' 

C 

C READ IN THE MATRIX DATA IN THE FORM X < CAT , ROW , COL I - 

— TF ‘ ( H0WRED;NE-.TT — G0-T0 _ 45 0 

READ ( NREAO ) X 

GO TO 451 

— '45 a "R E A D " I NRFSDTFNT! TTTXT I , J , tCT , T=! , 5T, J=1 ,MkOWLT,"K=T,MCOCO 

C 

c GO THROUGH THE INPUT DATA FOR EACH CATEGORY. 

A 51 DO i OO CAT*!", 'CLEW “ 

NUMRER*0 

I AA=0 

TAR=0 ' 

I AC = 0 

1 AD*0 ; 

J A AS=0‘ " 

I ARS=0 

|ACS»0 

rADS*0 

IARMsO 
IRCM = () 


H-3 


1CDM=0 
I ACM=0 

] HDM=0 _ _ _____ _ 

I ADM=0 “ 

C GO THROUGH THE INPUT MATRIX TO GST EACH DATA POINT IN THE CATEGORY. 

DO 101 COt_=l tMCOIL ___ _ 

DO 102 ROW = 1 , MROWL 

IE (X< l.ROVI T COU .NE.CAT ) GO TO 102 

NUMRPR=NUMBER+1 

I I = X ( 2 * ROd »COi. ) 

I? = X(3 , ROW, COL) 

13 = X(4,ROW,C0U 

14 = X I 5 ,ROW,CUU 
IAA=IAA+Il 

I AR= I A R+ ! 2 __ __ _ ' 

I AC = I AC +13"” 

I AD= l AD+I4 

I A AS= 1 AAS+ 1 1*11 

~ 1 ARS = I A BS + I 2# I 2 
IACS=IACS+I3*I3 

IADS = I A0S+I4»I4 

IABM=IAHM+I 1#12 
I RCM = I RCM + I 2*13 

ICDM= ICDM+I 3*14 

“ “rACM=IAC”M+Ti*i3 
IBDM=IRDM+I2*I4 

IADM = I AOM + I 1*14 

"102" CONTINUE 
101 CONTINUE 

FLAG I CAT ! =0 __ 

C“ THIS TEST IS Tfl PREVENT DIVISION BY ZERO. 

IF (NUMBER. NE.O) GO TO 400 

_F I. AG ( C AT) =1 

GO TO 100 

C FORM THE mean FOR EACH BAND. 

400 AN= I AA/NUMBER 

B N= j AB / NUMB ER 
CN= I AC/NUMBER 

. ON=I AD/NUMBER 

W RITE' (4 ,2 02) AN, BN, CN, ON, CAT, NUMBER 
C FORM THE MEAN MATRIX 

U( 1 ,CAT) =AN 

Urr,CAT)='8N 
U( 3,CAT)=CN 
U( 4,CAT )=0N 

z — form the cdvar'i ance~m a nrrx 

S(CAT,1,1 )=I AAS/NUMHER-AN*AN 
S (CAT, 1,2 )= I ABM/NUMBER- AN*BN 
SCCAT, 1 ,3 ) = I ACM/NUMKER-aN«CN 
5 (CAT , 1,41 = 1 A DM /NUMBER- AN *0N 
S ( CAT , 2 ? 1 ) *S ( CAT , 1 , 2 ) 

S(CAT,3, i )-S<CAT7r,3)'' 

S(CAT ,4,1 ) = S ( C AT ,1,4) 

S(CAT»2»2)=IABS/NUMHER-BN*BN 

'SI CAT';'2'»3 )"• I HCM/N(JrtHER-HN»CN 
S(CAT,2,4)=IB0M/NUMBER-RN*DN 
SIC AT ,3*2)=S1CAT ,2,3) 

S ( CAT ,4 , 2 ) = S ( CAT , 2 ,4 )"* " " . 

SI CAT, 3 ,3 ) = ! ACS/NUMBER-CN*CN 

S(CAT,3,4)=IC0M/NUMBER-CN*DN 
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S(CAT,4,3)=SICAT,3,4) 

S ( CAT ,4, 4 ) = I ADS/NUMRER-DN^DN 

100 CONTINUE 

— on 'i05”j«iT*r 

UCATEG(.n=U< J.CATEG) 

DO 106 X = l,4 

SK < J »K ) =S ( CAT6G, J,K ) 

106 CONTINUE 

103 CONTINUE ____ 

WRITE T6,2 03) CATE G”“ “ 

WRITE 16,204) (UCATEG{ 1 ) , < SK( I , J) ,J=1,4) , 1=1 ,4) 

l)l) 103 CAT=1,CLEN 

— ... j p (CAT.EO.CATEG) GO TO” 103 

I E ( FLAG ( CAT ) . EO . 1 ) GO TO 103 

_ 00 104_J=1,4 ; 

~ ■ U1 ( J )=U( J,CAT ) 

U2 ( J ) =U 1 ( J ) 

DO 107 <=1,4 

’Sit J » <T= S I C A T , J , K ) 

107 CONTINUE 

104 CONTINUE • 

— WRITE~nr,T(T5TXST 

WRITE (6,206) ( U 1 ( I ) , ( S I < 1 1 J ) , J = 1 *4 ) , 1 = 1 ,4 I 

C _ _ 

C SOfiT R~ACT — TMF~M FAKTITAT R I X . 

CALL GMSU6(U1,UCATEG,UA,4,1) 

C 

C lTT”NOW”B'ETrO'ME3 THE TR'ffTTS P U S fc UF UA. 

CALL GMTR A (UA,U1,4, 1 ) 

C 

c FOR’M~T H E pRODIJCT^TTICTh - TT”B"T I K ) IN 'ThE ECV pAp'E R . 

CALL GMPRD ( U A,U1 , BIJ I K ,4 , 1 ,4 ) 

c AOT) THr 3 R~T<ATRnrAW-B”TTATin x WhIChHB AS JUST form err: 

CALL GMAOO ( SK,BUIK , SKA, 4, 4) 

c 

c FORM “TRE INVERSE OFT hFS I “"RAT 111 X . 

CALL MINV (SI,4,DET,W0RKV1,W0RKV2) 

C ^ N Vg<SIK)+H( 1,0). - 

CALL GMPRD ( SI , SKA, T I ,4,4,4) 

C 

C FORM " THE E TCENVAUU E"S AND ' ETGEN VECTTTR'SC 

K=0 

00 113 .1=1,4 

DO 1 13— r=TTX 

K = K + 1 

113 B(K)=TI { I,.JI _ 

CALL EIGEN <rt,TI *4,0) 

C FORM A SCALAR TO NORMALIZE THE EIGENVECTOR. 

- SC ALA1=-T 1(1,1 )*UA ( 1 ) -T I I 7, l ) *UA ( 7 ) -TI T3 VT1 »CJA T3 ) ^TT ( 4 ; IT^UAT^T 

WRITE 1 6,223) SCAt.Al,R(l) 

IF (SCALAl.NE.O) GO TO 401 

--- WRITE (6,217) CAT, (TI ( T ,1 )',T‘=l'»4r 

GO TO 103 

401 00 114 1 = 1 ,4 

C FORM THE NORMALIZED EIGENVECTOR. * 

114 Til 1 ,1 J=T! (I , D/SCAI.A1 

WRITE (6,2<>U CAT , ( T I I I , 1 ) , I = 1 ,4 ) 
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WRITE <7,215) (02<J),J“1*4),CAT 
WRITS (7,215) (Till, 1), 1*1, 4) , CAT 

103 CONTINUE 

CAT IS USED" TO ''tDENT'lFY'THE 'CATEGORY IN DI FFEREFnT LOOPS THROUGH THE 
PROGRAM ANO ALSO AS A SUBSCRIPT FOR ARRAYS. 

CATSG IS THE CATEGORY TO RE ENHANCED. 

CLEN IS THE NUMBER OF CATERGDR I ES • 

COL f. ROW ARE USED AS SUBSCRIPT VARIABLES IN A LOOP TO GET THE DATA 
FROM THE INPUT MATRIX WITH MCOLL A ND MRDW t, BEING THE LIMITS OF CDL AND 

" ROW IN 'TtifiTOOP. ' “ ~ TJ _„ 

HOWREl) IS CODE FOR READING INPUT OATAi IF 1, READ FROM DISK UNEORMATED. 

IAA IS THE SUMMATION OF THE VALUESIN BANO_U 

1 All IS THE"SUMMATI0N'0F THE 'VALUES TN'~RAND 2. 

I AC IS THE SUMMATION OF THE VALUES IN BAND 3. 

I AD IS THE SUMMATION OF THE VALUES JN BAND_4. _ 

"I A AS- 1 AOS ARE THE SUMMAT fOfiS OF* THE "SOUARES OF THE INPUT NUMBERS. 
IABM-IADM ARE THE SUMMATIONS OF THE CROSS PRODUCTS OF THE INPUT NUMBERS . 
11-14 ARE THE INPUT VAI,UES IN E ACH LIGHT BAN D. _ 

'mcoll is the number op COLUMNS In the data matrix, 

NREAP IS USED TO SPECIFY INPUT TYPE. 

NROWL IS THE NUMBER_DF ROWS IN T HE DATA MATRIX. 

WMBEr ITThI’^ 'N llMHER' OF" VALUES A CATEGORY. 


200 FORMAT ( 1 3 , 3X , 13 ,3X , I 3 , 3 X , I 3 , 3X , I 3 , 3X , I 3 ) 

"201”' FORMAT" < ' 0" r ,3 X , • CATEGORY IS ',I2»6X,'TI=' f E14.7,3X,E14.7,3X,E14.7, 
13X , E L4.7 ) 

202 FORMAT (' ' , » A N« ' , Fft . 2 1 1 BNs»,Ffe.2,* CN^tFft.2,' PNs',Pfe.2,l C_ 

“ rA'TEG' 0 "RT~lT -r n 2T* FjUMB6R=STS1 

203 FORMAT I ,3X CATEGORY IS • ♦ 1 2 , 10X , 1 UCATEG ' , 10X , * SK MATRIX') 

204 FORMAT ('0'.27X.F6.2.12X,F6.2,5X,EA.2 t 5X,F4.2,5X,F6.21 


''205'" FORMAT T T:ri T3 X , '"C AT 6GDR Y I S 


► X, 'SI MATRIX' ) 


206 FORMAT I • 0 ' , 27X , F6 . 2 , 1 3 X , F6 . 2 , 5X , F A. 2 , 5X , F4 . 2 , 5X , F<S. 2 I 

210 FORMAT (20A4) 

TlT~FORMAT ( ' Six , r CA 1 TEG = ' ,7TP CLEN= 1 , I 2 , ' MCOLL= 1 , 12, 
12,' NREADa',12,* HOWREDs',12) 

212 FORMAT (' ',3X,'F0R CAT EGOR Y = 1 ,12, 1 THE SCALAR IS ZERO 
1 , 4 ( 3 X ,"E" 14 ,T)'Y 

215 FORMAT ( 614.7, 3X, 614.7, 3X,E14. 7, 3X,E14.7,3X, 12, 10X) 


MROWL= ' , I 


S 10X , ' T 1 = ' 


FORMAT ('O', 20X , 1 SC ALAI a • , F B . 4 , 1 0 X , 'EIGENVALUE 


' , E 14 .7 ) 


C B , WORK V 1 , ANO W0RKV2 ARE VORK AREA VECTORS. 

C HU I K REFERS TO THE 6(TK) IN THE EC V PAPER. 

c — f'lag t s "u s e 0 ”F w " at ’ c o no" iTi'o n’ coot; when a ctte'Gory rr 

c FMT IS READ IN AS THE FORMAT OF THE INPUT MATRIX X. 

C S REFERS TO THE COVARIANCE MATRIX IN THE ECV PAPER. 

c S I "RE F E R'S "TO”TKE"'S rn MAT'RTX'" I'N"'rR^"'Et:V~P'ffPll _ . 

C SK REFERS TO THE S(K) IN THE ECV PAPER. 

C SKA IS THE SUM OF S(K> AND H(IK> IN THE ECV y APER . 

C TT "I S THE" EIGENVECTOR 'W I TH' TH(E C/TR'KE ST' ETG EN V AIDE". 

c U IS A MATRIX USED TO FIND THE MEAN. 

C UA IS THE DIFFERENCE OF U(!> AND IKK). 

C'~ 'UC ATEG""REF ERS"TO THF1) ('! )'" IN~THF"FCV”p-A'P'Ein ’ 

C U1 AND U2 REFERS TO THE 11(1) IN THE ECV PAPER. 

C X IS THE INPUT MATRIX STORAGE AREA. 

- STOP " 

END 

//GO.SYSIN DO * 

//GO.FTlflFO 01 — DD DSN = C O'R'L'A'2 ,‘A S HCfiMD'lTFPlF-T 2,01 SPHEW 

1 2 4 4 5, 

(12,12,12, 12,12) 


USED. 



c 

c 



THIS PROGRAM USES THE EIGENVECTORS AND MEAN VECTORS ♦ FORMER IN AN EARLIER 
PROGRAM, TO FIND MINIMUM VALUES FOR THE INPUT MATRIX. 

INTEGER COUNT, FLAG, HOWRED, VECLEN . 

INTEGER*? COL, DEPTH, ROW " ' 

INTEGER* 2 CAT ( 99 I , CAT 2 I 12ft , 1 2« ) , XD ATA ( 5 , 12R, 12H ) 

DIMENSION’ AMEAN(99,4),DIFF(99,4) ,EVECT<99,4) _ __ 

DIMENSION FMTU 15) ,FMT2(15),FMT3< 1 5 ) , FMT4 (IS) ,'FMT'M 15) 

DIMENSION PROD I 99 | ,PROD2I 12fl,12R> 

C THE FIRST RATA CARD CONTAINS: CODE_ FOR TYPE OF INPUT AND OUTPUT, 

C THE DIMENSIONS OF 'THE INPUT DATA," AND' THE NUMBER OF INPUT RAT A BLOCKS. 

READ ( 5,200) NRE AD , ROW , COL , DEPTH, N WRITE, NUMX , HOWRED 

WRITE (6,402) NREAD, ROW, COL, DEPTH, NWRITE, NUMX, HUWREO 

C READ THE INPUT AND OUTPUT FORMAT S. 

READ (5,201) FMT1 

READ (5,201) -MT2 

~ READ < 5,201 V FMT3 
READ (5,201) -MT4 

_ READ (5,201) FMT5 _ 

WRITE (6,411 ) FMTL 
WRITE (6,412) FMT2 

WRITE (6,413) FMT 3 ___________ 
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C aTFROR V . ' 


WR I TE (5,400) PRODI I ) ,CAT< 1>,K1,K3,CAT2(K1.K3),I 
GO TO 302 

103 IP (PROO( 1) .GT.PR007 < K1 ,K3 )) GO TO 302 _ 

SAVE" THE MINIMUM PR0D2 AND THE RESPECTIVE 
PROD? ( K 1 » <3 ) = PROD ( I ) 

CAT2IK1 ,K3 )=CATII) _ 

302 CONTINUE 

301 CONTINUE 

300 CONTINUE 

WRITE (5,403) "NUM " 

WRITE ( NWR T TE , FMT4 ) 

WRITE (6, POP! 

WRITE (NWRITE,FMT5 ) 

IF (NUM.NE.NUMX) GO 


< ( PR0D2I 12,11), 1 1 = 1 ,CDL ), 1 2 = 1, DEPTH) 

< (C AT2TT27I n"7Tl"^CTTCT'n"2^r71TEPTH"T' 
TO 100 


C 

C' 

c 

C 

c 

c 

c 

C" 

c 

c 

c" 

c 

c 

"XT 

c 

c 

"c 

c 

c 

~C“ 

c 

c 

"C" 

c 

C 

c~ 


" AMEAN'"TS~ ^TH'E' ME AN ' VECTOR'S'. ‘ ’ 4 . 

CAT IS A VECTOR WITH CATEGORY I DENT I E IC AT I ON PROM THE FIRST PRUGRAM. 

CAT2 IS USED TO STORE THE CATEGORY AS SOCI A TED WITH THE MINI MUM PROD. 

COL IS A DIMENSION OF ’THE TNPilT MATRIX XDA T A . 

COUNT IS THE NUMBER OF EIGENVECTORS AND MEAN VECTORS. 

DEPTH IS a DIMENSION OP THE INPUT MATR IX XO AT A . . 

DIF F" TS " A VECT OR “0 F " T H E "f N pTi T” D"ATA‘"MTNiT5 Th t iTE S . 

EVECT IS the eigenvectors. 

flag is used to stop a read LOOP. . 

"'FMT 1 IS T H E FORM AT FOR R E A D I NG AMEAN . 

FMT2 IS THE FORMAT FOR READING EVECT, CAT, AND THE FLAG. 

FMT3 IS THE FORMAT_FOR READIN G T H E INPUT PATA (XDATA). 

fmT4~Ts TheoutpuT FORMAT for PROD 2. 

FMT5 IS THE OUTPUT FORMAT FOR CAT2. 

HOWRED IS CODE FOR READING JJVPIJT DATA , IF It READ FROM DISK FORMAT 40A? 
~N rB a D — rS" " a“ C’d OF F OR - T H E — T"y"p E — 0 F HiPt tTC 
NUM IS USED TO TEST FOR THE LAST DATA BLOCK. 

MO MX IS THE NUMBER OF INPUT DATA BLOCKS. 

~ NWR I TE I S - CODE FOR~"ThE TYPE" OF~OTJT pTH . 

PROD is the scalar product of THE EVECT vector and diff VECTOR. 

PR002 IS USED TO STORE THE MINIMUM PRODUCT I PROD) , 

~ ft OvT I S~A"D I MENS 1'0 N~OT" T'H 6 I N P I IT M A TRTX - XO'ATfiC 
VECLEN IS USED AS A LIMIT FOR LOOPS AND IS THE NUMBER OF LIGHT BANDS. 
XDATA IS THE INPUT DATA. 


tatthti"?', ' , 1 , I 2 * ' )'=" 


200 FORMAT ( 1 3 , 3X , I 3 , 3X , I 3 , 3X , 1 3 , 3X , 1 3 , 3X , I 3 , 3X , I 3 ) 

201 FORMAT (20A4) _ 

”400 FORMA TT r “ *”,'3X V ' P'ROD'A* ,'F6TZ , 1 U A I * ■ , X 2 , 1 

1.12,' I * • » I 2 ) 

401 FORMAT ( ' 0 ' » 10X * ' THFRE ARE TWO EQUAL MINIMUM VALUES.') 

-402 FORMAT-I • 1 ' , 3X V ’ NRE A0= '7T7V3XV •RTIW=-r7rT73XT , rm_ '=’ V 1 21 3X , « OEM 1 H* 1 , I 
L2,3X, 'NWRITE=' , 12 ,3X, 'NIJMX=* , 12, ' HOWRED* 1 , 12) 


403 FORMAT ( 'O' , AX , 'MINIMUM F(X) MATRIX 

404 FORMAT "( r 'VAX, » CATEGORY "MATR TX 'T' 

411 FORMAT (' «,' FMT1 ' » 20 A4 ) 

41? FORMAT (' ',' FMT? ',?0A4) 

413 FORMAT T T * V FMT 3 'T70 A4 F~ 

414 FORMAT (' FMT4 '.20A4) 

415 FORMAT (' ',' FMT5 ’,20A4) 

STOP” 

END 


FROM DATA SET NUMBER' ,13) 


//GO. SYS IN DO* 
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//GO. P T 18F001 00 0SN=C0BL&2. ASHLAND. S6P19-72, 01 SP=SHR 

6 5 4 4 6 1 

<614.7,3X,614.7.3X,614.7,3X,614.7 T 15X ) _ , 

«El4.7,3X t )rl4.7,3X t ei4.7 t 3X,fil4.7t3X,I2,2X,l3,SX> 

n?»i?»i2»i?»i2> 

{ • • *tL4.7,3X ,614.7 t 3X,614.7,3X,E14.7 ) 

(• * t I ? * 3X» 12»3X»I2*3X» 12) 

0.1 770000E+02 0.22500006+02 0. 1 520000E + 02 0.30500006+02 

0. 66877199-01 0.17276246-01 -0.55137576-02 .^•^2^4_936 r pj 2 — 999 

22420 731 
220311828 

_ 217191934 _ „ _ — 

214231227 

219231730 

2 9141224 __ _____ ; 

136301425 ' “ 

23223 931 

213271133 

129281026 

127231533 

13 224 934 

"213222530" 

216232237 

132 2514 29 __ _ - 

‘"'13423 934 
/* 


